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Geometric  Aspects  of  the  Linear  Complementarity  Problem 


by  Richard  E.  Stone 

ABSTRACT 


A  large  part  of  the  study  of  the  Linear  Complementarity  Problem  (LCP) 
has  been  concerned  with  matrix  classes.  A  classic  result  of  Samelson,  Thrall, 
and  Wesler  is  that  the  real  square  matrices  with  positive  principal  minors 
(P-matrices)  are  exactly  those  matrices  M  for  which  the  LCP  [q,  M)  has  a 
unique  solution  for  all  real  vectors  q .  Taking  this  geometrical  characteriza¬ 
tion  of  the  P-matrices  and  weakening,  in  an  appropriate  manner,  some  of  the 
conditions,  we  obtain,  and  study  other  useful  and  broad  matrix  classes  thus 
enhancing  our  understanding  of  the  LCP. 

In  Chapter  2,  we  consider  a  generalization  of  the  P-matrices  by  defining 
the  class  U  as  all  real  square  matrices  M  where,  if  for  all  vectors  x  within 
some  open  ball  around  the  vector  q  the  LCP  [x,  M)  has  a  solution,  then 
(q,  M)  has  a  unique  solution.  We  develop  a  characterization  of  U  along  with 
more  specialized  conditions  on  a  matrix  for  sufficiency  or  necessity  of  being 
inU. 

Chapter  3  is  concerned  with  the  introduction  and  characterization  of  the 
class  INS.  The  class  INS  is  a  generalization  of  U  gotten  by  requiring  that  the 
appropriate  LCP’s  ( q ,  M)  have  exactly  k  solutions,  for  some  positive  integer 
k  depending  only  on  M .  Hence,  V  is  exactly  those  matrices  belonging  to 
INS  with  k  equal  to  one. 


Chapter  4  continues  the  study  of  the  matrices  in  INS.  The  range  of 
values  for  k ,  the  set  of  q  where  (q,  M)  does  not  have  k  solutions,  and  the 
multiple  partitioning  structure  of  the  complementary  cones  associated  with 
the  problem  are  central  topics  discussed. 

Chapter  5  discusses  these  new  classes  in  light  of  known  LCP  theory,  and 
reviews  its  better  known  matrix  classes. 

Chapter  6  considers  some  problems  which  remain  open. 
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CHAPTER  1. 

BACKGROUND  TO  THE  LINEAR  COMPLEMENTARITY  PROBLEM 


1.1  Introduction 

The  central  topic  with  which  this  work  is  concerned  is  the  linear  com¬ 
plementarity  problem  (LCP).  The  LCP  is  a  nonlinear  system  of  inequalities 
where  we  are  given  as  data  an  n  X  n  real  matrix  M ,  a  real  n -vector  q ,  and 
are  asked  to  find  a  real  n -vector  z  such  that 

z>0,  (1.1) 

Mz  +  q>  0,  (1.2) 

*t(M*  +  $)  =  0.  (1.3) 

Although  we  shall  not  do  so  here,  one  can  consider  the  more  general  com¬ 
plementarity  problem:  given  a  closed  convex  cone  K  C  Kn ,  with  positive 
polar  cone  K*  =  {y  €  :  yTx  >  0,  for  all  x  €  K } ,  and  a  function 

F  :  K  -*  Rn ,  find  z  6  Rn  such  that 


*€K, 

(1.4) 

F{')  €  K* , 

(1.5) 

,*>(*)  «0. 

(1.6) 
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These  problems  may  be  thought  of  as  the  natural  formulations  to  use  in 
situations  where  an  equilibrium  point  is  being  sought.  They  arise  in  quite  a 
number  of  fields  including  engineering,  economics,  optimization,  game  theory 
and  control  theory.  For  more  on  these  applications  see,  for  example,  Lemke 
and  Howson  (1964),  Cottle  and  Dantzig  (1968),  Cohen  (1975),  Koehler  (1979), 
and  Cottle,  Giannessi  and  Lions  (1980). 

As  the  previous  references  would  suggest,  the  LCP  has  been  extensively 
studied.  Most  of  this  research  has  emphasized  the  algebraic  nature  of  the 
problem.  In  the  present  work  we  study  the  LCP  from  a  geometric  view¬ 
point.  Other  authors  have  also  taken  this  direction,  see  Saigal  (1970b,  1972a, 
1972b),  Murty  (1972),  Eaves  (1979),  Kelly  and  Watson  (1979),  Garcia  and 
Gould  (1980),  Howe  (1980),  Cottle,  von  Randow,  and  Stone  (1981),  and 
Doverspike  and  Lemke  (1981).  This  work  studies  the  characterization  and 
general  properties  of  matrices  M  for  which  (q,  M)  has  the  same  number  of 
solutions  “globally,”  and,  as  a  special  case,  has  a  unique  solution  “globally.” 
(Here  “globally”  is  used  to  mean  “for  all  q  €  9Rn  for  which  (q,  M)  has  a 
solution,  except  possibly  for  a  set  of  measure  zero.”  This  will  be  explained 
in  more  detail  later.)  Other  works  often  are  concerned  with  exhibiting  algo¬ 
rithms  that  “process”  the  LCP  for  a  specified  matrix  class,  and  then,  possibly, 
using  the  algorithm  to  show  various  properties  of  that  class.  In  this  work  we 
are  concerned  with  existence  proofs  and  properties  of  matrix  classes  rather 
than  with  algorithms.  The  questions  studied  do  not  seem  to  lend  themselves 
to  algorithmic  techniques. 

In  Chapter  2  we  will  study  LCP’s  which  have  either  no  solutions  or 
unique  solutions  at  almost  every  point.  We  will  derive  necessary  and  sufficient 
conditions  for  a  matrix  M  to  have  the  property  that  if  for  some  qo  €  SRn 
and  some  e  >  0  the  LCP  (1.1)— (1.3)  has  a  solution  for  all  S"  within 
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a  distance  of  e  from  q<j ,  then  the  LCP  has  precisely  one  solution  for  qo . 
Further  results  on  matrices  where  the  related  LCP  has  this  “global”  unique¬ 
ness  property  will  be  derived.  A  few  papers  show  that  some  known  matrix 
classes  are  of  this  type.  We  will  examine  these  papers  more  closely  in  Chap¬ 
ter  5.  In  another  direction,  the  question  of  local  uniqueness  in  an  LCP  was 
studied  by  Mangasarian  (1980).  That  paper  exhibits  necessary  and  sufficient 
conditions  for  a  solution,  z,  to  a  given  LCP  to  be  within  an  open  ball  in 
5Rn  that  contains  no  other  solutions  to  the  LCP.  Aside  from  keeping  an  al¬ 
gebraic  outlook,  these  results  are  in  a  different  vein  from  the  questions  we  are 
presently  considering  and  do  not  appear  to  be  helpful  to  the  current  study. 

In  Chapters  3  and  4  we  relax  the  condition  of  global  uniqueness.  In 
essence  we  replace  the  italicized  word  “ one  "  in  the  previous  paragraph  with 
k ,  where  k  is  some  fixed  positive  integer.  We  will  derive  characterizations 
for  these  matrices  and  related  results  concerning  the  geometric  structure  of 
LCP’s  with  this  property.  There  are  a  few  papers  that  deal  with  the  property 
of  an  almost  globally  invariant  number  of  solutions,  see  Murty  (1972),  Saigal 
(1972b),  Kojima  and  Saigal  (1979),  and  Mohan  (1978,  1980).  These  papers 
deal  with  special  matrix  classes  for  which  a  specialized  result  is  sought. 
They  do  not  attack  the  problem  in  full  generality,  and  some  do  not  look 
for  underlying  geometric  structure.  Saigal  (1972b)  contains  some  errors  - 
inherited  by  Mohan  (1978)  -  which  will  be  discussed  in  Chapter  5.  These 
papers,  along  with  some  others,  e.g.,  Saigal  (1972a),  discuss  the  property 
of  an  almost  globally  invariant  parity  in  the  number  of  solutions.  That 
is,  the  number  of  solutions  to  (1.1)— (1.3)  for  a  particular  matrix,  M ,  will 
be  either  odd  or  even,  not  both,  for  almost  all  q.  This  is  a  much  weaker 
property  than  that  of  an  invariant  number  of  solutions,  and  will  not  be  given 
much  consideration  here.  The  interested  reader  should  see  S jugal  (1972a) 
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for  a  complete  geometric  characterization  of  LCP’s  with  this  constant  parity 
property. 

Chapters  5  and  6  discuss  other  matrix  classes,  related  LCP  theory  and 
some  open  questions.  It  is  typical  for  dissertations  in  this  field  to  begin  with 
one  or  two  chapters  reviewing  the  known  classes  of  matrices  and  the  history 
of  the  area.  In  this  work  it  seemed  better  to  leave  this  to  a  later  chapter.  It 
is  Chapter  5  that  contains  such  a  summary. 

The  next  section  of  this  chapter  will  go  over  preliminary  results  that  are 
needed  throughout  this  work.  The  last  section  of  this  chapter  is  a  glossary  of 
the  notation  that  is  used.  It  is  suggested  that  the  reader  first  look  over  this 
last  section  to  see  the  basic  style  of  notation  used.  It  should  be  pointed  out 
that  throughout  this  work  the  word  interior  is  used  to  mean  relative  interior. 


1.2  Background  Material 


As  was  stated  before,  the  Linear  Complementarity  Problem  is:  Given 
M  €  5RnXn  and  q  €  ,  find  z  €  3ln  such  that 

*  >  0,  (1.1) 

Mz  +  q>  0 ,  (1.2) 

zt(Mz  -f-  q)  =  0 .  (1.3) 

The  LCP  with  M  and  q  as  the  data  will  be  denoted  as:  ( q ,  M) .  For  our 
purposes  it  will  be  useful  to  define  tv  =  Mz-\-q  .  Thus  we  can  express  (q,  M) 
as  the  problem,  given  M  6  SRnXn  and  q  G  5ln,  of  finding  z,w  6  Sn  such 
that 
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Iw  —  Mz  =  q , 
z,  w  >  0, 
zTw  —  0, 


(1.7) 

(1.8) 
(1.9) 


•where  I  is  the  n  X  n  identity  matrix.  This  formulation  of  the  problem  makes 
it  clear  that  we  are  just  trying  to  find  a  nonnegative  linear  combination  of 
the  column  vectors  of  I  and  — M  that  equals  q,  where  we  may  not  “use” 
both  I.i  and  — M.%  for  any  t  €  H .  This  idea  suggests  making 


Definition  1.1  For  M  E  9?nXn  and  a  E  (n)  define  Cm  (a)  E  3?nXn 


as 


CM[a)i 


,  if  *  a 

,  if  iE  a 


(1.10) 


where  the  subscript  M  will  be  dropped  when  it  is  clear  to  which  M  we  are 
referring.  The  are  called  the  complementary  matrices  associated  with 

M .  There  are  2n  such  matrices,  not  necessary  distinct. 


Associated  with  each  complementary  matrix  is  the  finite  convex  cone 


posC^a)  =  {y  E  »B  :y  =  CM(a)x,  x>0 }. 

The  cone  pos  is  called  a  complementary  cone  of  the  matrix  M ,  and 

the  subscript  M  is  dropped  when  it  is  clear  which  M  is  meant.  There  are  2n 
such  cones,  not  necessarily  geometrically  distinct.  Notice  that  two  distinct 
complementary  matrices  may  be  associated  with  complementary  cones  that 
are  geometrically  identical.  For  example,  the  matrix 

:i  (iu) 

will  have  pos  C({  1 })  geometrically  equal  to  pos  (7(5) ,  even  though  C7({  1 }) 
and  (7(3)  are  distinct  matrices. 
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If  we  have  /3  6  (n)  such  that 


dim[pos  C\i(a).p\  =  r 

then  pos  Cm{<*)p  is  referred  to  as  a  r- dimensional  facet  of  the  complementary 
cone  pos  Cm(oc)  .  Furthermore,  if  \/3\  —  n —  1  then  pos  Cm (a).p  is  referred 
to  as  a  face  of  the  complementary  cone  pos 

Let  8ol(<7,  M)  be  the  set  of  ordered  pairs,  (u>,  z) ,  of  solutions  to  the  LCP 
(9,  M) .  If  (w,  z )  6  sol(g,  M)  then,  letting  x  =  w  -\~  z  >  0  and  o  =  supp  z , 
we  have  C(a)x  =  q.  Conversely,  if  we  find  for  some  a  €  (n)  that  there 
is  an  x  >  0  with  C(ci)x  —  q  then  with  za  =  xa,  z&  =  0,  wa  —  0 
and  w&  =  x&,  we  have  (w,  z )  6  sol (9,  M) .  In  this  way,  each  solution, 
(t u,  z)  €  sol(g,  M) ,  will  be  associated  with  at  least  one  complementary  cone 
of  M .  Also,  in  this  way,  each  point  in  a  complementary  cone  of  M  will  be 
associated  with  at  least  one  solution.  We  can  now  state - 

Definition  1.2  For  M  €  SHnXn  let 

K(M)  =  U  pos  CM(a) . 

a€(n) 

We  then  see  from  the  previous  discussion  that 

K(M)  =  {g€S?n:sol(?,M)^0}. 

In  Figure  1.1  we  show  the  complementary  cones  for  the  matrix  in  (1.11). 
In  Figures  1.2  and  1.3  we  show  the  complementary  cones,  respectively,  for 
the  matrices  (1.12)  and  (1.13),  where 


■  — 1  — 1  ■ 

■  0  1  ■ 

.  -1  1 . 

.  1  1 . 

(1.12)  (1.13) 
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In  these  diagrams  the  column  vector  /.»  is  indicated  by  an  i  and  the  column 
vector  — M.i  is  indicated  by  an 

Each  solution  of  (q,  M)  must  be  associated  with  at  least  one  complemen¬ 
tary.  cone  containing  q ,  and  each  complementary  cone  containing  q  must  be 
associated  with  at  least  one  solution  of  ( q ,  M) .  However,  the  exact  relation¬ 
ship  between  complementary  cones  and  solutions  is  often  not  simple.  For 
example,  consider  the  problem  with  M  given  by  (1.12)  and  q  —  (1, 1)T . 
Then  q  is  contained  in  three  complementary  cones  pos  C(0) ,  pos  C({  1 })  and 
pos  C(2) .  However,  )sol(g,  M)|  =  2;  where  the  solution  («/,  z)  —  (1, 1, 0, 0) 
is  associated  with  pos  C(0) ,  and  the  solution  (w,  z)  =  (0, 0, 1, 0)  is  associated 
with  the  other  two  cones.  With  M  given  by  (1.13)  and  q  =  (0, — 1)T ,  we 
find  q  is  contained  in  the  complementary  cones:  pos  C( 2)  associated  with 
the  solution  (w,  z)  —  (0,0, 1,0);  posC({2})  associated  with  the  solution 
[w,  z)  =  (1, 0, 0, 1) ;  and  pos  C({  1 })  associated  with  the  infinitely  many  solu¬ 
tions  ( w ,  z )  =  (0, 0,1  -f-  0, 0) ,  where  0  ranges  over  all  nonnegative  reals.  In 
the  first  case  we  have  have  more  complementary  cones  containing  q  than 
solutions  to  ( q ,  M ) ;  in  the  second  case  there  are  more  solutions  to  (q,  M) 
than  there  are  complementary  cones  containing  q . 

To  help  in  our  discussion,  we  make  the 

Definition  1.3  For  M  €  3?nXn,  we  say  the  complementary  cone 
pos  Cm  (a)  is  full  or  nondegenerate  if  and  only  if  det  Cm  (a)  ^  0 ;  otherwise 
we  say  the  cone  is  degenerate.  Notice  detC^ffa)  —  ( — l)l“*  det  Maa  .  More 
over  a  complementary  cone  is  full  if  and  only  if  it  has  positive  n -dimensional 
volume,  and  a  complementary  cone  is  full  if  and  only  if  it  is  not  contained  in 
an  (n— 1) -dimensional  hyperplane.  In  addition  to  the  above,  we  say  M  itself 
is  nondegenerate  if  for  all  a  6  (n)  the  cone  posCjvf(a)  is  nondegenerate,  i.e., 


7 


all  the  principal  minors  of  M  are  nonzero. 

Definition  1.4  For  M  E  StnXn ,  we  say  the  degenerate  complementary 
cone  posC^O*)  is  strongly  degenerate  if  and  only  if  there  exists  a  2  6  Sn 
such  that  0  z  >  0  and  Cm  (a)  z  =  0 ,  i.e.,  if  and  only  if  for  q  =  0 ,  which 
is  in  every  complementary  cone,  we  find  ( q ,  M )  has  a  non-trivial  solution 
(w,  2)  7^  0  associated  with  pos  Cm{o)  ■  Otherwise  we  say  the  cone  is  weakly 
degenerate.  We  say  M  is  weakly  degenerate  if  not  all  of  its  complementary 
cones  are  nondegenerate,  but  none  of  the  complementary  cones  are  strongly 
degenerate.  That  is,  M  is  weakly  degenerate  if  it  has  a  zero  principal  minor 
and  sol(0,  M)  =  { (0, 0) } . 

Definition  1.5  For  M  6  5Kn><n,  we  say  (u/,  z)  E  sol (q,M)  is  a 
nondegenerate  solution  if  and  only  if  w  -f-  z  >  0 .  Otherwise  the  solution  is 
said  to  be  degenerate.  We  say  a  point  q  E  is  nondegenerate  with  respect 
to  M  if  all  solutions  to  (g,  M)  are  nondegenerate.  Otherwise  q  is  said  to 
be  degenerate. 

Consider  the  matrix  M  as  given  by  (1.11).  M  is  a  degenerate  matrix 
as  all  of  its  complementary  cones  are  strongly  degenerate,  except  for  the 
nondegenerate  complementary  cone  pos  C(0) .  The  matrix  M  as  given  by 
(1.13)  is  degenerate  as  it  contains  the  weakly  degenerate  complementary 
cone  pos  C({  1 }) .  Finally,  let  M  be  the  nondegenerate  matrix  given  by 

(1.12) .  If  q  =  (1,1)T,  then  q  is  degenerate  as  (w,  z)  —  (0,0, 1,0)  is  a 
degenerate  solution  to  ( q ,  M ) .  If  q  =  (3, 1)T ,  then  q  is  nondegenerate 
as  ( w 1,zi)  —  (3, 1,0, 0)  and  (w2,z2)  =  (0,0,  2,1)  are  the  only  solutions 
to  (q,  M)  and  both  are  nondegenerate.  The  reader  should  refer  to  Figures 
1.1,  1.2,  and  1.3,  respectively,  when  considering  the  matrices  given  by  (1.11), 

(1.12) ,  and  (1.13). 
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Now,  suppose  q  is  contained  in  the  interior  of  the  degenerate  complemen¬ 
tary  cone  pos  Cjvr(a) .  Thus  there  is  some  0  <  x  G  5Rn  such  that  C(a)  x  =  q . 
As  this  is  a  degenerate  cone,  there  exists  0  ^  y  €  5Rn  such  that  C(a)  y  =  0 . 
Thus  we  may  select  some  real  number  X  ^  0  such  that  0  x  -f-  Xy  >  0 . 
Hence,  if  we  let  za  =  (x  +  0y)a  ,  z&  =  0 ,  wa  =  0 ,  and  w&  =  (x  -f-  0y)&  , 
then  (iu,  z)  is  a  solution  to  ( q ,  M)  for  all  0  such  that  |0|  <  |X| .  Hence,  if  M 
is  degenerate  then  we  have  some  q  €  9 i  with  jsol(g,  M)\  —  oo .  Notice  also 
that  (q,  M)  has  a  degenerate  solution  when  we  let  0  =  X .  In  fact,  this  holds 
even  if  we  have  q  on  the  boundary  of  the  degenerate  complementary  cone. 
However,  now  we  might  have  X  =  0,  and  so  we  might  not  have  infinitely 
many  solutions,  but  we  will  still  have  a  degenerate  solution. 

Suppose  q  is  contained  in  the  nondegenerate  complementary  cone 
P03<7m(“)-  Thus  there  is  some  0  <  i  6  5?n  with  C(a)x  —  q.  But  now 
C(a)—1  exists  and  x  =  (^(a)-1?.  So  with  za  —  xai  z&  =  0,  wa  —  0, 
and  w&  —  x* ,  we  have  ([wa,  We,\,[za,  z&\)  is  the  only  solution  to  ( q,M ) 
associated  with  the  complementary  cone  pos  Cat  (a).  In  fact,  if  the  solution 
(u>,  z)  is  nondegenerate,  i.e.,  if  za  >  0  and  w&  >  0 ,  then  this  solution  is 
associated  with  no  other  complementary  cone.  For  if  it  were  associated  with 
pos  Cm{P),  P  €  (n),  then  we  would  need  2^=0  and  wp  =  0  which,  with 
the  previous,  would  imply  a  ==  /?.  We  now  have 


Proposition  1.6  Given  M  e  SRnXn: 


(i)  ( q ,  M)  has  finitely  many  solutions  for  all  q  6  5Rn  if  and  only  if  M  is 
nondegenerate; 

(ii)  if  q  6  9in  is  in  the  interior  of  a  degenerate  complementary  cone  then 
(q,  M)  has  infinitely  many  solutions; 

(iii)  each  degenerate  complementary  cone  is  associated  with  exactly  one 
solution  of  (q,  M)  for  each  q  (E  that  it  contains  (and,  of  course,  it 
is  associated  with  no  solutions  for  the  q  it  doesn’t  contain); 

(iv)  if  q  €  9?n  is  nondegenerate  then  there  is  a  bijective  correspondence 
between  solutions  of  (q,  M)  and  complementary  cones  containing  q . 

□ 


The  concept  of  complementary  cones  is  first  seen  in  Samelson,  Thrall 
and  Wesler  (1958),  and  was  later  given  a  comprehensive  treatment  by  Murty 
(1972).  Proposition  1.6  is  an  expansion  of  theorems  proved  in  Murty  (1972). 
Before  moving  on  to  discuss  other  areas  of  LCP  background  material,  it  is 
important  to  bring  up  the  following 

Definition  1.7  For  M  €  3?nXn,  we  say  the  two  complementary  cones 
posC(a)  and  posC(/3),  with  a,/9  6  (n),  are  adjacent  if  and  only  if 
|a  a  j3|  =  1 .  That  is,  two  distinct  complementary  cones  are  adjacent  if 
they  share  a  common  face.  If  a  a  /3  =  {i},  then  that  common  face  is 
pos  C(a).{  =  pos  . 

Definition  1.8  For  M  6  5RnXn,  we  say  the  common  face 
pos  Cm(o).%  between  the  complementary  cones  pos  Cm  (a)  and  pos  Cm  ((3) , 
where  a  a  /3  =  {»},  is  proper  if  and  only  if  (detCAf(a))(det(7M(/?))  <  0. 

As  det  CM(a)  =  (-l),a|  det  Maa ,  we  have 
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pos  \  is  proper  if  and  only  if  (det  Maa)(det  Mpp)  >  0. 

Geometrically,  pos  is  proper  if  and  only  if  it  is  (n—  1) -dimensional 

and  the  vectors  I.i  and  — M.,  lie  on  strictly  opposite  sides  of  span  Cm (a).j  . 

Definition  1.9  For  M  G  9?nXn,  we  say  the  common  face 
posC,A/(a).j  between  the  complementary  cones  pos  Cm  (a)  and  pos  Cm(P)  , 
where  a  A/3  =  {»},  is  reflecting  if  and  only  if  (det  C*M(«))(det  >  0. 

Similar  to  the  above  we  have 

pos  Cjif(a)  i  is  reflecting  if  and  only  if  (det  Maa)(det  Mpp)  <  0. 

Geometrically,  pos  (a).;  is  reflecting  if  and  only  if  it  is  (n — 1) -dimensional 
and  the  vectors  I.i  and  — lie  on  the  same  side  of  span  CV(a).| . 

Definition  1.10  For  M  €  9?nXn ,  we  say  the  common  face 
posC^f  (<*)•?  between  the  complementary  cones  pos  Cm  [a)  and  pos  Cu(fl), 
where  a  a/3  =  { i } ,  is  degenerate  if  and  only  if  (det  CM  (a))(det  Cm(P))  =  0 . 

As  above,  it  can  be  shown  that  posCm(c«)-»  is  degenerate  if  and  only  if 
(det  Maa)(det  Mpp)  =  0 ,  if  and  only  if  pos  Cjv^a).*  is  a  face  of  a  degenerate 
complementary  cone. 

For  examples  of  the  preceding  definitions  see  to  Figure  1.3,  which  shows 
K(M)  for  the  matrix  (1.13).  Here  pos  C(0).i  is  proper,  posC($).2  is  reflect¬ 
ing,  and  pos  (7(5).!  and  pos  C(0).2  are  degenerate. 

We  now  move  on  to  consider  the  algebraic  concept  of  principal  transforms 
of  the  matrix  M .  For  a  more  detailed  discussion  see  Tucker  (1960,  1963), 
Cottle  and  Dantzig  (1968),  and  also  Cottle  (1974).  Suppose  we  are  given  a 
matrix  M  €  3?mXn  which  is  not  necessarily  square  and  can  be  permuted  to 
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look  like 


M  = 


M&0  \ 

.  I 

6  SRmXn 


Moreover,  suppose  a  €  (m),  /?€(«),  ja|  =  \/3\,  and  det  Map  y^O.  We  then 
say  the  matrix 


-M-p1  Ma4 

M&0 

M&$  ~  M&p  Mafl 

(1.14) 


is  a  pivotal  transform  of  M .  We  also  say  ~M  is  gotten  from  M  by  block 
pivoting  on  Ma0 .  If  a  =  ,  we  then  say  M  is  a  principal  transform  of 

M .  Notice  from  (1.14)  that  if  a  C  y  €  (»)  then  the  principal  transform  of 
Min  resulting  from  a  block  pivot  on  Maa  is  just  (M)~n  •  In  other  words, 
the  principal  transform  of  a  submatrix  will  be  the  submatrix  of  a  principal 
transform.  (The  converse  is  not  necessarily  true.)  The  following  two  theorems 
are  straightforward  algebraic  consequences  of  the  definition  of  A?.  They  can 
be  found,  for  example,  in  Cottle  (1974)  and  Parsons  (1970). 


Theorem  1.11  Given  M  6  5Rrn><n  with  M  6  9?mXn  being  the  trans¬ 
form  of  M  obtained  by  block  pivoting  on  Map ,  then  for  all  z  (=  5Rn  and 
y  £  we  have 


if  and  only  if 


Map 

tt&fi 

ya 


H 


Xfi 


y& 


□ 

Theorem  1.12  (Tucker  (I960))  Given  M  6  5fcn*n  and  a  €  (n).  If 
€  5RnXn  is  the  principal  transform  of  M  obtained  by  block  pivoting  on 
Maa  ,  then  for  all  /3  £  (n) 


det  Mpp 


det  A fa&P,a&0 
det  Maa 


□ 

We  will  now  obtain  a  few  facts  concerning  principal  transforms  and  their 
relation  to  the  LCP. 

Theorem  1.13  Given  M  €  SRnXn  and  q  £  5Rn,  consider  the  matrix 
[  M  |  q  ]  6  5RnX(n+1)  and  let  [  M  |  q  )  £  sjjnx(n+i)  be  its  principal 
transform  obtained  by  blocking  pivoting  on  Maa  for  some  a  £  (n).  Then 
|sol(g,  M)\  —  |sol(g,  M)\ . 

Proof.  From  Theorem  1.11  we  have  for  any  w,z  £  3?n  that 


Maa 

Map 

9a 

Mpa 

M&& 

9* 

to<* 


(1.15) 


1 


f 
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(I 


i 


! 1 
hi 


if  and  only  if 


M  act 

A fa& 

q* 

Za 

*  a 

z& 

M&a 

A 

9a 

1 

(1.16) 


Hence,  if  ((toa,u>&I>  [za,  *&\)  €  soI(g,  M)  then  ([za,  tn&J,  [tua,  z&j)  €  sol(q,  M) , 
and  vice  versa.  This  gives  us  a  bijective  correspondence  between  solutions  to 
(q,  M)  and  solutions  to  (q,  M) .  Thus,  the  number  of  solutions  must  be  the 
same  for  the  two  LCP’s. 

□ 


Theorem  1.14  Given  M  6  3£nXn,  q  6  and  a  6  (n),  let 
(  M  |  q  1  €  3fnX(n+1^  be  the  principal  transform  of  [  M  |  q  ]  €  SRnXn 
obtained  by  block  pivoting  on  Maa  •  Then  q  €  int  K(M)  if  and  only  if 
q  €  int  K (M)  . 

Proof.  For  any  z,  w,  x  €  SRn ,  Theorem  1.11  implies  that 


r 

A faa 

Ma& 

Qa 

Iaa 

0 

M&a 

M&& 

q& 

0 

I&& 

Za 

z& 

1 

Xa 

x& 


if  and  only  if 
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wa 


u>& 


(1.17) 


f 


span  5Rn.  Suppose  g  €  K(M)  .  We  then  have  an  t  >  0  such  that  x  £  B(q,  e) 
implies  sol(g -f  x,  M)  0 .  Let  x  —  (— M~*  xa,x&)T  £  3in.  Thus,  by  (1.17) 
and  (1.18)  we  see  that  sol(q-\-z,  implies  sol(g-f-i,^7)  0 .  As  (1.19) 

spans  5Rn ,  the  set  of  S  corresponding  to  all  x  6  B(q,  e)  contains  an  open  ball 
around  q.  Thus  q  £  int  K(M)  .  This  proves  one  direction  of  the  theorem. 
The  other  direction  is  proved  by  the  same  argument. 


Theorem  1.15  Given  M  £  5ftnXn,  q  6  3in  and  a  £  (fx),  let 
(  A?  |  q  ]  £  $ftnX(**+1)  be  the  principal  transform  of  [  M  |  q  ]  6 
obtained  by  block  pivoting  on  Maei .  Then  q  £  intpos  Cm(0)  if  and  only  if 
q  £  int  pos  Cjf(a  a  /3) . 

Proof.  Suppose  q  £  intpos Cm(0).  Then  there  is  an  e  >  0  such  that 
x  £  B(q,e)  implies  q  -j-  x  £  pos  Cm(P),  the  latter  thing  implies  there  is 
a  (to,  z)  £  sol(g  -}-  x,  M)  such  that  top  =  0  and  —  0 .  As  before, 
let  2P  =  ( — xa,xt,)T  6  9in  •  By  (1.17)  and  (1.18)  we  see  that  (57, 1)  — 
{[za,  toa],  [to0,  a^J)  6  sol(5-|-S,K7).  Hence,  with  T  =  oAj3,we  have  57^  =  0 
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i 


and  Z*,  =  0 .  This  means  q  -\-  x  £  pos  C-^i) .  Thus  for  the  set  of  x 
corresponding  to  all  x  £  B(q,e),  which  as  before  will  contain  an  open  ball 
around  q ,  we  have  q  x  £  pos  Cjf(a  a  0) .  Hence,  q  £  int  pos  C(a  a  0) . 
The  other  direction  of  the  theorem  is  proved  by  the  same  argument. 

□ 


The  preceding  theorems  show  that,  from  the  standpoint  of  combinatorial 
topology,  the  structure  of  K(M)  is  identical  to  the  structure  of  K(A?) .  The 
positive  orthant  in  K(M)  k  identified  with  pos  C(a)  in  K (M) .  Pivoting  on 
Maa  is,  in  essence,  “swapping"  the  vectors  I.a  with  the  vectors  — M.a . 

As  our  last  topic,  *>'»'  turn  to  look  at  some  classes  of  matrices  that  we 
will  need.  We  will  be  discussing  many  more  matrix  classes  in  Chapter  5,  but 
for  now  we  will  mention  only  the  classes  P,  Po ,  Q,  Qo ,  and  Eo  . 

We  say  a  matrix  M  £  tftnXn  is  in  P  (Po)  if  and  only  if  all  its  principal 
minors  are  positive  (nonnegative).  It  is  clear  that  membership  in  P  or  P0 
is  an  inherited  property,  i.e.,  if  a  matrix  is  in  P  (Po)  then  all  its  principal 
submatrices  are  in  P  (Po ).  We  also  see,  from  Theorem  1.12,  that  if  a  matrix 
is  in  P  (Po)  then  all  its  principal  transforms  are  in  P  (Po)*  (This  was 
first  proved  in  Tucker  (1963).)  The  main  theorem  concerning  the  geometric 
structure  of  P-matrices  comes  from  Samelson,  Thrall  and  Wesler  (1958)  and 
states 


Theorem  1.16  For  M£fRnXn, 
for  all  q  £  9Jn . 


M  £  P  if  and  only  if  |soI(g,  M)\  —  1 

□ 


Another  pur  of  matrix  classes  that  are  defined  by  the  LCP  are  Q  and 
Q0.  A  matrix  M  £tRnXn  is  said  to  be  in  Q  if  and  only  if  sol(g,  M)  0  for 
all  q  £  5tn ,  i.e.,  K(M)  =  3?n .  It  is  clear  that  P  C  Q.  However,  the  aero 


16 


matrix  in  any  dimension  is  in  Po  but  not  in  Q.  Also,  the  matrix 


'  1  2 

.2  1  . 


(1.20) 


is  not  in  Pq  ,  so  not  in  P,  but  it  is  in  Q,  as  can  be  seen  in  Figure  1.4. 


The  definition  of  Qo  requires  the  concept  of  being  “feasible”  with  respect 
to  a  LCP.  We  say  z  €  9in  is  feasible  with  respect  to  (q,  M)  if  and  only  if  z 
satisfies 

z  >  0,  (1.1) 

Mz  +  q>  0.  (1.2) 

We  now  can  define  M  6  9?nXn  to  be  in  Qo  if  and  only  if  for  all  q  €  3tn ,  for 
which  there  is  a  z  €  SR  which  is  feasible  for  (g,  M) ,  we  have  sol(g,  M)  0 . 
Clearly  Q  C  Q0  .  Notice  that  the  negative  of  the  identity  matrix  in  any 
dimension  is  not  in  Po ,  but  is  in  Qo  as  then  (g,  — I)  has  a  feasible  z  if  and 
only  if  g  >  0 ,  in  which  case  (g,  0)  €  sol(g,  M) .  Also  the  matrix 


’  0  1  ' 

0  1  . 


(1.2X) 


is  in  Po  but  not  in  Qo,  for  when  g  =  ( — 1,0) T  we  find  sol(g, 
although  z  =  (0,  l)T  is  feasible.  From  Eaves  (1971),  we  have  the  following 
geometric  result  pertaining  to  Qo -matrices. 

Theorem  1.17  Given  M  €  5RnXn,  Me  Qo  if  and  only  if  K(M)  is  a 
convex  set  in  Sftn  . 


Before  leaving  this  section,  we  mention  the  matrix  class  Eo  •  A  matrix 
M  6  5ftnXn  is  said  to  be  semi-monotone,  denoted  M  E  Eq  ,  if  and  only  if 
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for  all  x  €  3in ,  where  0  7^  x  >  0 ,  there  is  an  index  k  €  n  for  which  X*  >  0 
and  (Mi)it  >  0 .  (This  class  was  introduced  in  Eaves  (1969).)  Consider  the 
matrices 


0  1 
.  1  0  . 


—1  2 

— 1  1. 


(1.22) 


(1.23) 


Notice  that  the  matrix  (1.22)  is  in  Eo  but  not  in  Pq  .  It  isn’t  in  Qo  since 
z  —  (1,  Of  is  feasible  for  q  =  (1,  — l)r ,  yet  there  is  no  solution  to  the  LCP 
with  this  q  and  matrix  (1.22).  Also,  matrix  (1.23)  is  in  Q,  as  can  be  seen  in 
Figure  1.5,  but  is  not  in  Eo .  It  is  fairly  obvious  that 


M  >  0  =»■  M  €  E0  .  (1.24) 


It  is  also  fairly  obvious  that  being  in  Eo  is  an  inherited  property,  i.e.,  if  a 
matrix  is  in  Eo  then  so  are  all  its  principal  submatrices.  For  if  the  vector 
x  €  5Rn  with  xa  >  0  is  such  that  ( Maa  xa)k  <  0  for  all  A:  €  a  where 
I*  >  0 ,  then  letting  x&  =  0  we  have  0  x  >  0  with  (Mx)*  <  0  for  all 
k  €  (n)  where  x*  >  0.  A  is  less  obvious  fact  is  the  following  (see  Fiedler 
and  Ptak  (1966),  Lemke  (1970)  and  Eaves  (1971)). 

Theorem  1.18  Po  C  Eo  . 

□ 


1.3  Notation 

For  easy  reference,  this  section  lists  the  notation  that  will  be  used  in  this 
work  and  specifically  the  notation  which  is  not  standard. 
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Item 


Explanation 


n  The  set  { 1, 2, 3, . . . ,  n  } . 

a,  /?,  'i ,  etc.  Index  sets.  Example:  the  ordered  A; -tuple  a  — 

(a:i,...lOffc)  with  1  <  ai  <  •••  <  a*  <  n. 

(n)  The  collection  of  all  index  sets  formed  from  n  (in¬ 

cluding  the  empty  set,  0  ). 

&  The  index  set  “complementary”  to  a  (relative  to 

?j).  a  is  obtained  from  (1, 2, ...,n)  by  deleting 
the  components  in  a . 

t  q  for  ot  =  { t } . 

sgmxn  The  class  of  all  real  m  X  n  matrices. 

2+  The  class  of  all  positive  integers. 

Ma(}  The  submatrix  of  M  with  rows  indexed  by  a  and 

columns  indexed  by  P .  If  a  =  P  we  say  then  call 
Maa  a  principal  submatrix  of  M .  The  determinant 
of  a  principal  submatrix  of  M  is  called  a  principal 
minor  of  M .  By  convention  det  M%%  —  1 . 

AC*  (Map)"1- 

The  itk  row  of  M . 

M.j  The  jth  column  of  M . 

Ma.  The  rows  of  M  indexed  by  a. 

M.0  The  columns  of  M  indexed  by  /? . 

t 

1 
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The  entries  of  the  vector  z  indexed  by  a . 


I 


Za 

CM{a)  or  C{a) 


span  C 
aff  X 


dimX 


posC 


K(M) 


A  complementary  matrix  relative  to  M  and  the 
index  set  a .  If  C  =  C(a) ,  then 


if  j  €  a 
if  j  £  a 


The  subscript  M  is  normally  dropped  when  it  is 
clear  from  the  context. 

The  column  space  of  the  matrix  C . 

The  affine  hull  of  the  set  X .  This  is  the  set 
{x  4"  0(y  —  x)  :x ,  y  €  X,  and  0  6  5R  }  . 

The  dimension  of  the  affine  hull  of  the  set  X  .  This 
is  the  minimum  number  of  columns  needed  in  a 
matrix  C  so  that,  given  some  q  €  X,  we  have 
affX  =  {q  z  :  z  €  span  C  } . 

The  set  { Cx  :  i  >  0 }  where  C  is  a  matrix 
(not  necessarily  square).  If  C  is  a  complemen¬ 
tary  matrix  relative  to  M  and  some  a  6  (n), 
then  pos  C  is  called  a  complementary  cone.  A  com¬ 
plementary  cone  is  said  to  be  full  or  nondegenerate 
if  det  C  0 .  Otherwise,  it  is  called  degenerate. 

The  set 

(J  posCWte). 

<*€{*) 
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The  set 


<(M) 

B(q,  e) 

int  X 

dX 

X 

P 

Po 

Q 

Qo 

Eo 

1*1 

(q,M) 


U  P°s  CM(a).t 

«e(5T) 

iSn 

The  open  ball  centered  at  q  with  radius  e .  This  is 
the  set  {  x  G  9?n  :  ||a:  —  g||  <  e  }  . 

The  relative  interior  of  the  set  X  with  respect  to 
aff  X .  This  is  the  set  of  all  q  G  X  such  that  there 
is  some  e  >  0  such  that  aff  X  n  B(q,  e)  C  X . 

The  relative  boundary  of  the  set  X  with  respect  to 
affX.  Thus  dX  =  X\intX. 

The  closure  of  the  set  X .  This  is  the  set  of  all 
z  G  5Rn  such  that  for  all  e  >  0  there  is  a  q  E  X 
where  q  6  B(z,  e) . 

U„{  M  6  5RnXn  :  det  Maa  >  0,  for  all  a  6  (n)  } . 

Un{  M  6  SRn><n  :  det  Maa  >  0,  for  all  a  6  (n) }  . 

Un{  M  G  3 ?nxn  :  K(M)  =  3Rn  }  . 

Un{  M  G  9?nXn  :  K(M)  =  pos[  1 1  —  M  J  } . 

Un{  AT  G  3S"Xn  :  O^i  >  0  =» 

3fc  xk  >  0  &  (Mi)fc  >0}. 
Matrices  in  this  class  are  said  to  be  semi-monotone. 

The  cardinality  of  the  set  X . 

The  LCP  given  by  (1.1),  (1.2)  and  (1.3). 

21 


sol(g,  M ) 

The  set  of  all  solutions  of  the  LCP  (q,  M) . 

a  a  ft 

The  symmetric  difference  of  a  and  /3 .  This  is  the 

set  (a  \  f3 )  u(/3  \  a) . 

supp  z 

The  support  of  the  vector  z .  This  is  the  set 

{j  :  Xj-^0}. 

One  last  point  before  ending  this  list:  we  say  a  set  X  is  star-shaped  at 
q  if  and  only  if  for  every  z  £  X  we  have 

{  Xqr  -f-  (l  —  \)z  :0<\<1}CX. 

This  says  that  the  line  segment  between  q  and  z  is  contained  in  X. 
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► 

2,  1* 


Figure  1.2 
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2 


Figure  1.5 


CHAPTER  2. 

THE  CLASS  OF  U-MATRICES 


2.1  Preliminary  Definitions  and  Results 

In  Chapter  1  we  exhibited  several  matrix  classes  that  are  related  to 
the  LCP.  It  is  often  the  case  that  one  useful  class  of  matrices  leads  to 
the  consideration  of  other  interesting  matrix  classes,  gotten  by  weakening  or 
strengthening  the  conditions  that  define  the  original  class.  For  example,  the 
class  P  suggests  considering  the  more  general  class  Po .  The  class  P,  viewed 
as  the  class  of  all  matrices  M  for  which  (q,  M )  has  a  unique  solution  for  all 
q ,  suggests  defining  the  class  Q  by  dropping  the  uniqueness  requirement  and 
just  requiring  that  for  each  q  a  solution  to  (q,  M)  must  exist.  The  class  Q, 
in  turn,  gives  rise  to  the  class  Qo  when  we  relax  the  definition  to  require  only 
that  (q,  M)  have  a  solution  whenever  (1.1)  and  (1.2)  alone  are  satisfiable. 

We  presently  wish  to  understand  the  basic  geometric  structure  which 
gives  rise  to  unique  solutions  to  ( q ,  M) .  With  this  in  mind,  we  consider  the 
following  class  of  matrices 
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(J  {  M  6  »nxn  :  |sol(g,  M)\  =  1,  for  all  g  €  K(M)  }  . 

n 

This  matrix  class  is  obtained  from  P  by  relaxing  the  requirement  that  ( q ,  M ) 
have  a  solution  for  all  g,  as  Q  was  obtained  by  dropping  the  uniqueness 
requirement.  However,  as  we  will  see  later,  this  “new”  matrix  class  consists  of 
nothing  but  P.  While  this  is  an  insightful  result  by  itself,  a  subtler  weakening 
of  the  definition  of  P-matrices  is  needed  to  get  an  appropriate  matrix  class 
for  our  analysis.  We  find  that  the  appropriate  class  to  study  is  embodied  in 
the  following  definition. 

Definition  2.1.  A  matrix  A  will  be  said  to  be  a  U-matrix,  A  6  U,  if 
and  only  if 

A  €  |J  {  M  €  »nXn  :  |sol(<7,  M) |  =  1,  for  all  q  €  int  K(M)  }  . 

n 

The  goal  of  the  next  section  will  be  to  develop  a  characterization  for  the 
class  U.  Before  embarking  on  this  task,  we  give  two  examples  which  verify 
that  U  consists  of  more  than  just  P;  we  also  discuss  some  needed  definitions 
and  results. 

Example  2.2  Let 


In  this  case, 

K (M)  ={g€3i2  :gi  +  ?2  >  0} 

as  shown  in  Figure  2.1.  Note  here  that  (q,  M)  has  a  unique  solution  for  all 
q  satisfying  gi  +  ga  >  0  including  those  of  the  form 

gi  r  0  I 

<li  >  0  and  ga  >  0 

.  0  J  ga 
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for  which  the  solution  to  ( q ,  M)  is  degenerate. 

Example  2.3  Let 

In  this  case, 

as  shown  in  Figure  2.2.  Here  the  problem  has  a  unique  solution  for  all  q 
satisfying  q  >  0  or  q  <  0 . 

Notice  that  K (M)  is  convex  in  Example  2.2  and  nonconvex  in 
Example  2.3.  In  both  cases,  |sol(9,  M)|  =  oo  for  all  9  G  dK(M)  . 

Perhaps  the  most  important  fact  underlying  the  study  of  uniqueness  is 
expressed  by 

Lemma  2.4  If  M  G  9Rn><n ,  the  following  are  equivalent: 


0) 

Mg  Eq  (that  is,  M  is  semi- monotone); 

(ii) 

(9,  M)  has  a  unique  solution  for  all 

9  >  0; 

(iii) 

for  all  a  G  (n) ,  the  system 

Mo,a®Q  <  0, 

®«  >  0 

M  = 


0  1 
1  0 


K(M)  =  X2+  U  »L 


has  no  solution. 


□ 


The  equivalence  of  (i)  and  (ii)  was  shown  by  Eaves  (1971).  The  equiv¬ 
alence  of  (i)  and  (iii)  was  shown  by  Lemke  (1970). 

Since  int  C  int  K(M)  for  any  M  €  RnXn,  it  follows  immediately 
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U  C  E0  . 


(2.1) 


from  the  definitions  that 


However,  the  inclusion  is  proper  as  shown  by 


M  = 


1 

.  2 


2  ' 
1  . 


€  Eq 


for  which  K(M)  =  5ft2.  In  this  instance  ME  Eo  as  M  >  0.  ME  Q, 
as  seen  by  Figure  2.3,  so  every  point  q  E  9?2  is  interior  to  K(M)  .  But 
some  problems  ( q ,  M)  do  not  have  unique  solutions,  for  otherwise  M  would 
belong  to  P  which  it  does  not. 

Let  M  E  9inXn  and  q  6  3in  be  given.  If  the  matrix  [  M  \  q  ] 
is  a  principal  transform  of  [  M  |  q  ]  then,  by  Theorems  1.14  and  1.13, 
respectively,  we  know  that  q  E  int  K(M)  if  and  only  if  q  E  int  K(M)  ,  and 
that  |sol(q,  M)\  —  |soI(q,  M)| .  From  this  we  find 


M  EV  ME  U. 


This  leads  us  to  the  following  definition. 

Definition  2.5  If  M  E  SRnXT*r  we  say  M  is  fully  semi-monotone  if 
and  only  if  every  principal  transform  of  M  is  semi-monotone.  We  denote 
the  class  of  such  matrices  by  Eg . 

We  remark  that  Eg  C  Eo  as  the  “empty  pivot”  is  always  legitimate: 
M  is  always  a  principal  transform  of  itself.  Notice  that  being  in  Eg  is  an 
inherited  property  of  matrices.  For,  from  Chapter  1,  we  know  that  being  in 
Eo  is  an  inherited  property,  and  also  that  a  principal  transform  of  a  principal 
submatrix  will  be  a  principal  submatrix  of  a  principal  transform. 
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The  matrices  used  in  Examples  2.2  and  2.3  show  that  Eq  is  a  nonempty 
class.  As  a  matter  of  fact,  Eq  contains  Po .  This  follows  as  any  principal 
transform  of  a  Po -matrix  belongs  to  Po,  and  as  Po  C  Eo  .  The  matrix 
used  in  Example  2.3  shows  that  Po  C  Ej  is  a  proper  inclusion. 

Our  remarks  above  the  definition  imply  that 

U  C  E£  (2.2) 


which  strengthens  (2.1).  But,  again,  the  inclusion  is  proper.  Indeed, 


G  Eq  * 


but  with  q  —  (1, 0)  the  problem  (q,  M)  has  the  solutions 


(u/\*1)  =  (1,0, 0,0) 

(w2,z2)  =  (0,0, 0, 1)  . 

The  corresponding  cone  K(M) ,  shown  in  Figure  2.3,  is  quite  revealing.  Notice 
that  int  K(M)  contains  the  interior  of  the  degenerate  complementary  cone 
pos  C({  2  }) . 


2.2  Characterisation  of  U-matrices 

We  have  seen  in  the  last  section  that  U  C  Eq  .  The  task  now  is  to  find 
precise  conditions  under  which  a  matrix  in  Eq  will  also  be  in  U.  It  turns 
out  to  be  easier  to  state  exact  conditions  for  when  an  Eq -matrix  is  not  in 
U .  The  main  result  of  this  section  is: 
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Theorem  2.6  Let  M  £  SR”*".  Then  M  £  U  if  and  only  if  either 
M  £  Eq  or  there  exist  a,(3  £  (n)  and  t,  j  £  n  such  that 

(i)  otj&P,  ij&j, 

(ii)  (detM^aXdetM^)  7^0  and  there  exists  a  nonzero  vector 
v  £  3?n  such  that  vTC(a).i  =  vTC(/3).j  =  0, 

(iii)  there  exists  x  £  5Rn— 1  with  x  >  0  and  C(a).j  x  £  posC(j3).j  . 


Taken  together,  conditions  (i),  (ii),  and  (iii)  say  that  there  are  two  full 
complementary  cones  which  have  an  (n  —  1) -dimensional  intersection  on  two 
differently-labelled  faces. 

To  prove  Theorem  2.6,  we  first  prove  two  lemmas. 

Lemma  2.7  Let  M  £  StnXn  .  M  £  Eq  if  and  only  if  for  all  a,P£  (n) 
with  det  Maa  ^  0  and  a  f3  we  have 

int  pos  C(a)  D  pos  C{0)  =  0 . 


Proof.  Let  [  M  |  q  ]  be  the  principal  transform  of  [  M  |  q  ]  gotten  by  block 
pivoting  on  Maa .  We  know,  by  Proposition  1.15,  that  q  £  intposCM(a) 
if  and  only  if  q  £  int  pos  C^(0)  [if  and  only  if  q  >  0].  If  we  assume  that 
M  £  Eq  ,  then  A?  £  Eq  .  Letting  C  =  and  using  Proposition  1.13  with 
Lemma  2.4  we  conclude  that 

q  £  int  pos  C7(a)  =»  |sol(qr,  M)\  =  1 .  (2.3) 

For  q  £  int  pos  C(a),  we  have  C(a)~1  q  =  *  >  0  giving  the  solution 
(tu,  z)  €  *ol(g,  M ) ,  where  za  =  xa  >  0  and  u/4  =  >  0 .  If  q  £  pos  C(/3) , 

31 


1- 


then  there  is  a  solution  (w,  z)  €  sol(g,  M)  with  wp  =  0 ,  and  as  a^j3,  we 
have  ( w ,  z)  z)  contradicting  (2.3). 

Conversely,  if  we  knew  that  int  pos  C(a)  intersected  no  other  com¬ 
plementary  cones,  thm,  as  above,  i  =  C(a)~~ 1  q  would  give  us  a  solution 
to  ( q ,  M ) ,  and  it  would  be  the  only  solution.  Thus  (2.3)  is  valid;  again  by 
Proposition  1.13  and  Lemma  2.4,  we  have  M  €  Eo  .  Since  this  holds  for  all 
a  €  (n)  for  which  det  Maa  ^  0 ,  we  have  M  €  Eq  . 

□ 

The  preceding  lemma  says  that  when  M  6  Eq  ,  no  point  in  the  interior 
of  a  full  complementary  cone  lies  in  any  other  complementary  cone. 

Lemma  2.8  Let  L  be  an  n -dimensional  linear  subspace  of  9JP,  and  let 
A  and  B  belong  to  W*™  where  m  >  n.  If,  for  i,j  €  m, 

(i)  span  A  =  span  B  —  L , 

(ii)  int  pos  A.%  fl  int  pos  B.j  ^  0 , 

(iii)  span  A. j  =  span  B.j  L , 

(iv)  A.%  and  B.j  lie  on  the  same  side  of  span  B.j  (relative  to  L ), 
then 

int  pos  A  fl  iiit  pos  5^0. 

Proof.  From  (ii),  we  have  the  existence  of  a  positive  vector  x  £  3R2m—a  such 
that 

[Aj,  —B.j]  z  =  0,  x  >  0 .  (2.4) 

If  the  conclusion  were  false,  there  would  be  no  vector  x  such  that 

=  0,  S>0. 
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Then,  by  Stiemke’s  alternative  theorem  (see  Dantzig  (1963))  there  would  exist 
a  vector  u  such  that 

O^u t[A,-B\  >  0. 

(Without  loss  of  generality,  we  may  assume  u  6  L  .)  But,  by  the  same  alter¬ 
native  theorem,  the  existence  of  a  solution  to  (2.4)  implies  the  nonexistence 
of  a  solution  to 

0^uT[A.it- B.j]  >  0. 

From  this  we  deduce  that 

uT  A.i  =  uT  B.j  =  0. 

Thus  u  is  orthogonal  to  the  span  of  A.j  (which  equals  the  span  of  B.j  ).  Yet 
uTA.i  >  0  >  uTB.j  .  Thus  A.i  and  B.j  lie  on  opposite  sides  of  span  A3  , 
since  by  (iii)  neither  can  lie  in  span  B.j  ,  a  contradiction. 

□ 

We  remark  that  this  lemma  could  be  made  stronger;  e.g.,  we  could  allow 
A  £  3?pXr ,  A  G  9?pXa ,  i  £  f ,  j  £  3  and  replace  (i)  and  (iii)  with 

(i')  span  A  C  span  A  =  L 

(iii' )  span  A  7^  span  A.j  C  span  B.j  7^  L . 

However,  stronger  results  are  not  needed  in  what  follows. 

Proof  of  Theorem  2.6:  Sufficiency.  As  we  have  already  noted,  UC  Ej, 
so  M  £  Eq  implies  M  jj  U.  Suppose  then  that  M  £  Eq  and  the  three 
conditions  of  Theorem  2.6  hold.  Let  a,  (3 ,  i,  j ,  v ,  and  1  be  as  described 
therein.  Define 

H  =  {q  :  vTq  =  0} . 


•Jt 
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Then  by  (ii) 

pos  C(a).i  U  pos  C(/3).j  C  H . 

Clearly  H  is  (n  —  1) -dimensional.  By  (ii),  pos<7(a).i  and  pos  C(/?)j  are 
also  (n  —  1)  -dimensional.  Condition  (iii)  implies 

int  pos  C(a). 5  fl  pos  C(/3)  j  ^  0 . 

In  fact,  the  stronger  assertion 

int  pos  C(a).i  fl  int  pos  C(/3).j  7^  0  (2-5) 

must  also  hold.  To  see  this,  let  q  =  C(o!).{  x  .  As  q  is  interior  to  posC(a).}, 
the  dimension  statements  above  imply  that  for  some  e  >  0 ,  all  points  in  H 
within  a  distance  e  from  q  belong  to  int  pos  C(a). i  .  But  clearly  pos  , 

which  lies  in  H ,  contains  interior  points  within  e  of  q  ;  hence  (2.5)  is  valid. 

Certainly  C(a).i  and  C(P).j  do  not  lie  in  H .  If  they  lie  on  the  same 
side  of  H ,  then  as  pos  C(a)  and  posC(fi)  are  full  cones,  Lemma  2.8  implies 
that  int  pos  C(a)n  int  pos  C(/9)t>^0,  contradicting  Lemma  2.7.  So  C(a).,- 
and  C(P).j  lie  on  opposite  sides  of  H .  Hence 

int  pos  C(a).i  D  int  pos  C(/?).  j  C  int{  pos  C{a)  U  pos  C(ft) }  C  int  K(M)  . 

Let  *7  =  a  A  {  i  } .  Then 

C( 7).*  =  C(a).t  and  C(7).,  =  C(P).<  ■ 

Since  iy^j ,  we  have  pos  C(7)  C  span  C(a).t  ;  this  implies  pos  C('y)  is  a 
degenerate  cone.  However, 

int  pos  C(a).j  flint  pos  C(/3).j  C  posC(7). 
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As  this  intersection  is  nonempty,  there  exist  points  of  pos  C(i)  in  int  K(M)  , 
hence  there  are  points  q  E  int  pos  C(i)  n  int  K(M)  ,  and  so  (q,M)  will  have 
more  than  one  solution.  That  is,  M  £  U. 

Necessity.  We  assume  M  (2  U.  Then  |sol(g,  M)\  >  1  for  some  q  belonging 
to  int  K (M)  .  Considering  what  must  be  proved,  we  assume  M  6  Eq  and 
show  that  the  three  conditions  are  satisfied.  There  are  two  cases. 

Case  1:  q  is  in  the  intersection  of  two  full  complementary  cones.  Assume 
for  the  moment  that  one  of  the  cones  is  pos  0(0) ,  i.e.,  the  nonnegative 
orthant.  Let  pos  C(n)  be  the  other  cone  where  Hy£0 .  Then  there  exists 
a  unique  vector  x  >  0  such  that 

C(n )  x  =  q  >  0 . 

If  Xp  —  0,  i.e.,  the  solution  does  not  use  any  columns  from  — M  but  only 
columns  from  I,  then  by  the  uniqueness  of  x,  we  have  x  —  q,  and  the 
solutions  that  arise  from  0(0)  and  0(n)  are  the  same.  If  x^^O,  we  may 
assume  >  0 .  (If  it  is  not,  we  may  replace  n  by  <j  =  supp  x .  Then 
C(a )  x  —  0(fj.)  x  .  If  pos  C(a)  is  degenerate,  the  argument  of  Case  2  applies.) 
Thus,  as  C(fi)  =  q ,  we  have 

—  —  qp  ^  0,  ^  0 . 

But  detAfMM7>^0,  and  €  Ep  as  M  €  Eq  .  Therefore  having 

— AfMftxM  >  0  with  xM  >  0  says,  with  respect  to  the  LCP  (q^,  M^^) ,  that 
an  interior  point  of  a  full  complementary  cone  is  contained  in  3?!^  ,  another 
complementary  cone.  This  contradicts  Lemma  2.7. 

For  two  full  complementary  cones,  say  pos  0(\)  and  pos  O(n) ,  the  ar¬ 
gument  just  given  can  be  made  to  apply  by  performing  a  principal  pivot 
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on  M\\ .  (Let  the  resulting  matrix  be  M  and  use  the  cones  posC~^(0), 
pos  and  the  correspondence  between  the  cone  structures  of  K(M) 

and  K(M) .)  Either  way,  Case  1  cannot  occur. 

Case  2:  q  belongs  to  a  degenerate  cone.  We  now  assume  det  M ^  =  0 

and 

q  6  pos  C[fj)  Hint  K(M)  .  (2.6) 

Let 

dim  pos  C(n)  =  a,  0  <  a  <  n . 

Note  that  if  a  =  0,  then  C(/x)  =  — M  =  0.  But  then  M  belongs  to  U. 
From  (2.6)  we  have 

dim[pos  C{p)  D  int  K(M)  ]  —  a . 


Thus 


lim  |(J( 


dim  <  M(pos  C{n)  f)  int  K(M)  (1  pos  C7(X)|  :  det  C(X)  7^ 


°1 


—  a 


as  int  K(M)  is  contained  in  the  union  of  the  full  complementary  cones.  Since 
the  union  is  finite,  there  exists  a  /?  €  (n)  with  det  C{(3)  7^  0  and 


dim(pos  C(n)  fl  int  K(M)  n  pos  C(0)]  =  a  . 


Lemma  2.7  says  pos  C(n)  D  int  pos  C[f3)  7^  0 ,  so 

pos  C(/x)  D  int  K(M)  Pi  pos  C(/3)  C  d  pos  C(f3) . 
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Since  pos  C{n)D  pos  C(/3)  is  a  convex  cone  and  intpos  C(/3)  C  int  K(Af)  ,  it 
follows  that 


pos  C(n)  fl  int  K(M)  n  pos  C(/3)  C  pos  C(/3)j 


for  some  j .  As 

dim  pos  C(n)  =  s  =  dim[pos  C[n)  n  pos  , 


we  have  posC(n)  C  span  0(0).$  .  The  (n  —  1) -dimensional  subspace 

H  =  span 

is  the  common  boundary  of  the  two  closed  half-spaces  Jf  +  and  H~ .  Let 
H+  contain  C(/3).j  .  Now 


whence 


Suppose 


int  K(M)  fl  pos  C(/3).}  ^  0 , 


dimfint  K(M)  fl  pos  C(/3).j]  =  n  —  1 . 


detC(X)^0 

pos  C(X)  n  int  H  ~  0 


=>  dim[pos  C7(X)  n  int  K(M)  D  pos  C(/9).j J 


<  n  —  1. 


Then  there  exists  q  6  int  K(M)  D  pos  C(/3)j  contained  in  no  full  cone  that 
intersects  H~ .  Thus,  there  exists  a  number  cq  >  0  such  that  for  all 
e  €  (0,  co) 

B(e,  q) Dint H~  flint  K(M)  =0 
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since  int  K(M)  is  in  the  union  of  the  full  complementary  cones.  But  as 
q  €  H ,  we  have  B(e,  q)  D  int  H~  ^0.  Thus  int  K(M)  does  not  contain  an 
open  ball  around  q  €  int  K(M)  ,  a  contradiction.  This  implies  there  exists 
a  6  (n)  with  det  C(a )  ^  0 ,  pos  C(a)  H  int  H~  ^  0 ,  and 

dim[pos  C(a)  D  int  K (M)  ft  pos  C(/3).j]  =  n  —  1 . 

Since  pos  C(P)  C  H~^~ ,  it  is  clear  that  a^p.  Again 

int  pos  C(a )  f|  pos  C7(/3).j  =  0 


by  Lemma  2.7,  so 


pos  C(ac)n  pos  C(fl). }  C  d  pos  C(a). 

As  before  (with  posC(/i)),  we  must  have  pos  C(a)  n  pos  C(P).3  lying  in  an 
(n —  l)-face,  say  posC(a).j,  of  pos  (7(a).  But 

dim(pos  (7(a). j  n  pos  C(P).  3]  =  n  —  1  (2.7) 


and 


dim(pos  C(a).j]  =  dim[pos  C{P).}]  =  n  —  1 ,  (2.8) 


so 

pos  C(a). j  U  pos  C(p).j  C  H . 
Pick  orthogonal  to  H.  Then 

vTC(a).  j  =  vTC(P ).j  =  0. 


38 


Notice  that  C(a).i  €  int//  (for  otherwise,  posC(o) flint//  =  0).  Thus, 
as  C((3).j  6  int  H+ ,  and  C(n).j  6  H ,  we  have  if^j-  In  light  of  (2.7)  and 
(2.8),  there  exists  a  vector  x  €  SRn— 1  such  that 

C7(a).|  x  e  pos  C(/?).3  ,  x  >  0 . 

This  completes  the  proof. 

□ 

Notice,  from  the  proof  of  sufficiency,  that  all  degenerate  cones  of  a 
U-matrix  must  be  in  dK(M)  . 


2.3  Variations  on  the  Characterisation  and  Further  Results  on  U-matrices 

In  the  previous  section  a  set  of  necessary  and  sufficient  conditions  was 
given  for  a  matrix  not  to  be  in  U.  These  conditions  describe  U  as  a  subclass  of 
Eq  by  stating  exactly  what  “goes  wrong”  with  an  E{>  -matrix  when  it  is  not 
in  U.  It  is  of  interest  to  look  at  other  (sufficient)  conditions  on  an  Eg  -matrix 
that  would  “force”  it  out  of  U  vis-4- vis  (necessary)  conditions  that  would 
have  to  hold  were  the  matrix  not  in  U.  This  will  give  us  a  better  idea  of  the 
structure  of  U-matrices,  especially  by  looking  at  why  other  conditions  are  not 
both  necessary  and  sufficient.  With  this  in  mind,  we  have 
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Theorem  2.9  If  M  E  E£  n3fnXn  and  there  exist  a,f3 ,7  E  ( n )  such 

that 

(i)  a4^  =  {i}^{i}  =  aA1, 

(ii)  det  Maa  =  0  and  (det  Mppftdet  M^)  >  0, 

(iii)  C(P).i  and  C(i).  j  are  on  opposite  sides  of  span  C{a) 

[that  is,  with  1  =  <?(/?).»  ,  y  =  ,  and  >1  =  C(a).\  , 

the  inequality  x T{I  —  A(AT A)-1  AT)y  <  0  holds], 

then  M  £  U . 

The  basic  idea  here  is  that  if  in  K(M)  we  have  two  nondegenerate  cones 
“sandwiching  in”  a  degenerate  cone,  then  the  matrix  cannot  be  in  U. 

Proof.  By  (i)  we  have  that  C(/3).?  =  C(a). j  and  C( 7). 3  =  C(a). 3 . 
Since  det Maa  =  0 ,  we  then  have  a  vector  u^O  such  that 

vTC(P).i  =  vTC{ 7).j  =  0.  By  Theorem  2.6  it  remains  to  show  that 
int  pos  C(j3).i  fl  int  pos  Cfa). 3  0 .  Suppose  not.  Since  pos  C{/3). j  and 
posC^.j  lie  in  the  same  (n  —  1) -dimensional  subspace,  and  since 
pos  C(/3)t-  =  pos  C(t f).t-  ,  it  follows  from  Lemma  2.8  that  and  C(qf).< 

lie  on  opposite  sides  of  span  C(/3).,-  .  (Notice  that  C(a).t-  =  (7(/3).t*  = 
C( 7),,-  .)  Thus  there  exists  a  positive  number,  0,  a  nonzero  vector 
v  E  span  C(a),  and  vectors  z,z  E  5Rn—a  for  which 

<?(/?)•;=  +  » 

^(7) •»  =  s  —  Ov. 

From  (iii),  there  exists  a  positive  number,  r ,  a  nonzero  vector  u;  E  ,  and 
vectors  x,  x  €  SRn— 1  for  which 

C(/3).i  =  C(/3).jx  +  u; 

^(7)-;  =  ^(7)5  Z—TW. 
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Thus, 


det  C(p)  =  det[  C(p)r3  \  C(P).~z  +  v  \  C(p).-,x  +  «,  ] 

where  the  matrix  is  represented  with  column  i  on  the  right,  column  j  in  the 
middle,  and  all  other  columns  on  the  left.  Hence 

det  C{p)  =  det[  C{P),  -  |  C(p)^  z  +  v\  w] 

■■=  det[  C{(31  |  v  |  w  ] 

=  ~  det[  C{p). »  |  w  ]  v  ] 

==~h  CW+3  I  ~TW  I  ~6V  1 

=  det(  3  I  — ' ™  I  C(7).C]  *~6v\ 

=  ~h  det^  C^  y>  I  *  *~TW  I  c(i)-$  z  —  9v  J 

=  -^detC(7). 

This  contradicts  (ii),  so  our  supposition  was  false  and  the  theorem  follows. 

□ 

Recall  that  Po  C  Eq  .  If  we  apply  Theorem  2.9  to  a  Po -matrix, 
the  inequality  “  >  ”  in  (ii)  can  be  replaced  by  the  symbol  “  7^  ”  and  the 
condition  more  closely  resembles  that  in  Theorem  2.6.  Notice,  in  the  proof  of 
Theorem  2.6,  that  if  we  knew  M  €  Po  and 

C(a).}  £  span  C(p).,~  , 
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we  could  define 

a  —  /3  a  {  j  } 

?=  /9 

t  =  {*>.?} 

and  then  a,  0,  and  7  would  satisfy  the  hypotheses  of  Theorem  2.9.  All 
we’ve  done  is  verify  that  C(/3).t»  and  C(a).j  together  constitute  a  linearly 
independent  set  of  columns  so  that,  when  C(a).,  is  adjoined,  a  nondeg(  a.te 
complementary  cone  with  the  desired  position  is  formed.  Thus,  we  need  to 
have  C(a).j  €  spanC(^).,^  and  C{0)\  E  span  C(a).t~  to  spoil  this  reasoning. 
It  seems  plausible  that  the  conditions  of  Theorem  2.9  are  necessary  as  well 
as  sufficient  for  ME  Po  \  U .  However,  this  is  not  the  case. 

Example  2.10  Let 


0  ' 
— 1 
0  . 


It  is  easily  checked  that  M  €  Po  •  The  only  full  complementary  cones 
corresponding  to  it  are  posC(0)  and  posC(3)  -  i.e.,  pos/  and  pos — M . 
The  hypotheses  of  Theorem  2.9  cannot  be  satisfied  by  this  matrix  since  the 
index  sets  0  and  7  can  differ  by  only  two  elements.  But  M  £  U  as,  for 
q  =  (1, 1,  Of  €  int  K(M) ,  the  problem  (q,  M)  has  the  solutions 


(1,1, 0,0, 0,0) 

(w2,z2)=  (0,0, 0,0, 1,1). 


This  example  can  be  used  to  disprove  the  necessity  of  the  conditions  in 
Theorem  2.9  because  we  can  construct  the  two  full  cones,  that  “sandwich 
in”  the  degenerate  cone,  so  that  their  index  sets  differ  by  more  than  two 
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elements.  (Clearly  they  must  differ  by  at  least  two  elements.)  Thus,  we  might 
consider  combining  conditions  (i)  and  (ii)  from  Theorem  2.6  with  condition 
(iii)  of  Theorem  2.9  to  get 

Corollary  2.11  If  M  €  Eg  nSRn><n  and  M  U  then  there  exist 
a,  P  6  (n)  and  i,  j  £n  such  that 

(i)  a^p,  iy^j, 

(ii)  (det  Maa)(det  Mpp)  0  and  there  exists  a  nonzero  vector 
v  £  S?n  such  that  vTC(a). j  =  vTC(p).j  —  0, 

(iii)  C{a).i  and  C{P).j  are  on  opposite  sides  of  span  C(a).{  =  span  C(p).j 
[that  is,  {vT C(a).i)(vT C(P).j)  <  0  ]. 

Proof.  This  follows  immediately  from  Theorem  2.6,  Lemma  2.7  and 
Lemma  2.8.  For  (i)  and  (ii)  are  from  Theorem  2.6,  and  if  C(a).i  and  C{p).j 
were  on  the  same  side  of  span  C(a).j ,  then  condition  (iii)  of  Theorem  2.6  and 
Lemma  2.8  together  would  imply  that  int  pos  C[a )  D  int  pos  C(P)  0 ,  which 
would  contradict  Lemma  2.7. 

□ 

To  show  that  these  conditions  are  not  sufficient  for  M  not  being  in  U, 
we  have 

Example  2.12  Let 


•  0 

1 

0 

M  = 

0 

0 

1 

.  1 

0 

0  . 

This  matrix  belongs  to  Po ,  and  so  is  in  Eg .  The  only  full  complementary 
cones  are  pos  7  and  pos — M  which,  in  this  case,  meet  only  at  0.  Thus, 

int  K [M)  =  int  pos  I U  int  pos  —  M 
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and  clearly  M  6  U.  Yet  the  three  conditions  mentioned  above  are  satisfied 
as 

span  C(0),3  =  span  C(3).j 

and  C(0).8  and  {2(3). i  lie  on  opposite  sides  of  spanC(0).3  =  spanC(3).j  . 

Another  possible  variation  of  Theorem  2.6  would  be  to  make  condition 
(iii)  much  stronger.  This  would  clearly  preserve  the  sufficiency  of  the  condi¬ 
tions,  giving  us 

Corollary  2.13  Let  M  €  Eq  n  3inXn  .  If  there  there  exist  a,  /3  €  (n) 
and  i,  j  €  n  such  that 

(i)  a^/3,  i^j, 

(ii)  (det  MaQ)(det  Mpp)  0  and  there  exists  a  nonzero  vector 
v-(E  9?n  such  that  vTC(a).i  =  vTC(/3). j  =  0, 

(iii)  pos  C{a).t  C  pos  C(p).3  , 
then  M  fc  U . 

□ 

However,  this  new  condition  (iii)  is  too  strong  to  be  necessary,  as  is  shown 
by 

Example  2.14  Let 

0  0  0  —1 

0  0  0  —1 

M  = 

0  0  0  1 

l  1  0  0  0 

Obviously  ME  Po  C  Eg  ,  and  the  only  full  complementary  cones  are 
posC(0)  and  posC({  1, 4}).  They  intersect  only  on  the  respective  faces 
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posC(0).4  and  pos  C[{  1, 4  }).j  .  (Notice  that  spanC(0).4  =  span  C({  1, 4  }).j  .) 
For  x  =  (1,2, 1)T  we  have 


c({  1,4 


so  the  two  faces  do  have  (relative)  interior  points  in  common.  Hence  M  j?  U 
by  Theorem  2.6.  Now 


and 


But 


and 


C(0).4  = 


1 
0 
0 
L  0 


0  0 
1  0 
0  1 
0  0 

0 
1 
0 
0 


1 

1 

— 1 
0 


0  0 
1  0 


0 
L  0 

1 
0 
0 
L  0 


1 

0 

0 

1 

0 

0 


1 

1 

— 1 


0 

0 

1 


0  ■ 

0 

£  pos 

L  l  . 

0  1 
0 
1 

0  J 

1 

0 

0 


0 

1 

0 

0 


r 1 

0 

£  pos 

L  o  . 

0  0 
1  0 
0  1 
0  0 


0 
0 
1 

0  J 

1 
1 

— 1 
0  J 


so  neither  face  contains  the  other. 

We  now  examine  the  state  of  affairs  for  vectors  q  6  dint  K(M) 
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Theorem  2.15  If  M  6  Eon9inXn  and  q  €  dintK(M),  then 
|sol(g,  M)\  —  oo .  (In  fact,  sol (q,M)  is  unbounded.) 

Proof.  We  know  q  must  lie  in  some  (n  —  1) -dimensional  face  of  dint  K(M)  . 
Since 

g€dintK(M)  C  {  d  pos  C(a)  :  det  C(a)  0  }  , 

or£(n) 

q  must  belong  to  an  (n  —  1) -dimensional  face,  C(a).j,  of  some  full  com¬ 
plementary  cone  pos  C(a)  such  that 

dim[pos  C(a).j  n  d  K (M)  ]  =  n  —  1 . 

The  union  of  all  points  in  pos  C(a).{  fl  9K(M)  that  are  contained  in  a 
fc -dimensional  complementary  cone  with  k  <  n  —  2,  that  are  contained  in 
the  boundary  of  an  (n  —  1) -dimensional  face  of  a  complementary  cone,  or 
that  are  contained  in  an  (n —  1) -dimensional  face,  of  a  complementary  cone, 
not  contained  in  span  C{a).%  is  a  finite  union  of  sets  of  dimension  n  —  2  or 
less.  Hence,  we  can  find  a  point  q  €  int  pos  C(a). i  fl  d  K(M)  not  in  this  union 
that  is  arbitrarily  close  to  q .  If  q  €  pos  for  some  j  6  n ,  f)  7^  a ,  with 

det  C{0)  0 ,  then  posC(/3)-3  ^  span  C(a).j  and  q  £  int  pos  C(P).3  .  So  we 

have  either  C(a).*  and  C(/?).y  on  the  same  side  of  span  C(o:).{ ,  which  by 
Lemma  2.8  implies  that 


pos  C(q)  fl  pos  C(/3)  0 

contradicting  the  assumption  that  M  €  Eq  ,  or  else  we  have  C(a).t  and 
C{fi).j  on  opposite  sides  of  span  C(a).j  which  implies  that 

q  6  int[pos  C(a)  U  pos  C(/3)\  C  int  K(M)  , 
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contradicting  the  fact  that  q  G  dK(M)  .  So  q  is  contained  in  only  one 
(n  —  1) -dimensional  face  of  one  full  cone.  From  Lemma  3.2  of  Saigal  (1972a), 
we  see  that  q  must  be  contained  in  some  complementary  cone,  pos  C(/3) , 
where  C(/3 )  x  =  0  for  some  nonzero  i  >  0 .  As  there  are  finitely  many  such 
c ones,  and  as  q  was  arbitrarily  close  to  q,  we  can  find  a  sequence  qv  —*  q 
in  some  such  cone.  As  all  such  cones  are  closed,  we  may  assume  without  loss 
of  generality  that  q  G  posC(/3).  Thus  q  =  C(/3)y  for  some  y  >  0,  and  for 
each  X  >  0,  y  -f-  \x  will  give  us  a  different  solution  to  (q,  M) . 

□ 


Theorem  2.15  explains  why  we  must  define  U  with  respect  to  the  interior 
of  K(M),  rather  than  all  of  K(M).  If  ( q,M )  has  a  unique  solution  for 
q  G  K(M) ,  then  certainly  M  G  U.  But  Theorem  2.15  then  requires  that 
dintK(M)  =  0.  Thus  we  must  have  int  K(M)  =  3?”  thus  MG  Q. 
However,  U  D  Q  —  P  which  gives  us  nothing  new.  In  fact,  the  proof  of 
Theorem  2.15  shows  that  if  q  G  d  int  K(M)  then  q  is  in  a  strongly  degenerate 
cone.  Thus  we  have 


Corollary  2.16  If  M  G  Ej  ,  then  9intK(M)  is  contained  in  the 
union  of  the  strongly  degenerate  cones. 

□ 

Corollary  2.17  If  M  G  U  and  M  is  nondegenerate  or  weakly  degen¬ 
erate,  then  MG  P  • 

□ 
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2.4  Eq  fl  Qo  -matrices  and  U  D  Qo  -matrices 


In  this  section  we  confine  our  attention  to  those  matrices  within  Eg 
and  U  which  are  also  in  Qo .  (Recall  that  M  G  Qo  if  and  only  if  K(Af)  is 
convex.)  We  start  off  with  a  lemma  used  to  prove  the  next  (familiar)  theorem. 
It  expresses  the  underlying  structure  of  <9K (M)  for  Uf)  Qo  -matrices. 

Lemma  2.18  Suppose  M  G  Uf|  Qo  n9?,lXn  and  let  pos  (7(a)  be  a  full 
complementary  cone  relative  to  M .  Define  the  index  set  =  a  A  { i }  .  Then 
span  (7(a). j  is  a  supporting  (boundary)  hyperplane  of  K(M)  if  and  only  if 
C{0).i  lies  in  span  C(a). j  . 

Proof.  If  C(/3).i  G  C(a).i  ,  then  pos  C(j3)  is  a  degenerate  cone.  Therefore 
pos  C(j3)  G  5K(M)  as  M  G  U  .  Since  pos  C(/3)  C  span  C(a).; ,  we  see  that 

dim[span  C(a).j  fl  d  K  (M)  ]  =  n  —  1 . 

Thus  span  C(a).{  is  a  supporting  hyperplane  of  the  finite  convex  cone  K(M)  . 

Conversely,  suppose  C(/3).i  £  span  C(a).j  .  If  C'(a).,  and  C(/?).j  were 
on  the  same  side  of  span  C(a).{  ,  then  by  Lemma  2.8,  the  interiors  of  the 
full  complementary  cones  posC(a)  and  pos  C(/3)  would  intersect.  This 
contradicts  that  M  6  Eq  ■  Thus,  C(a). j  and  C(/3)  i  are  on  opposite  sides 
of  span  (7(a).{  .  Hence  we  have 

intpos  C(a).}  C  int(posC(a)Upos(7(/3)j  C  int  K(M)  , 

so  span  C{a).\  cannot  be  a  supporting  hyperplane  of  K (M) . 

□ 
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Without  the  assumption  that  M  6  U  we  find  that  both  directions  in 
Lemma  2.18  fail  to  hold.  The  matrix 

0  — 1  ' 

10.’ 

see  Figure  2.3  again,  is  in  Qo  (T  Eq  but  is  not  in  U.  As  always,  pos  (7(0) 
is  a  nondegenerate  complementary  cone,  and  C({2}).2  €  span  (7(0). j  .  But 
span  C(0).2  is  not  a  supporting  hyperplane  to  K(M) .  For  the  other  direction, 
consider 


This  matrix  is  also  in  Qo  but  not  in  U.  In  this  instance,  K(M)  =  . 

Lemma  2.18  would  make  each  boundary  hyperplane  contain  a  degenerate 
complementary  cone.  But  this  is  clearly  not  the  case. 

Notice  that  the  second  of  these  matrices  is  not  in  Eq  ,  and  cannot  be  as 
the  second  part  of  the  proof  only  needed  M  £  Eq  .  Notice,  also,  that  the 
first  of  these  two  matrices  belongs  to  Pq  ,  but  not  the  second.  In  fact  we 
prove 

Corollary  2.19  Suppose  M  £  Po  ft  Qo  D3?nXn  and  let  pos  (7(a)  be 
a  full  complementary  cone  relative  to  M.  Define  the  index  set 
/3  =  a  a  { i }  .  If  span  (7(a). {  is  a  supporting  (boundary)  hyperplane  of 
K(Af) ,  then  C(P).i  €  span  C(a). j  . 

Proof.  Suppose  C7(/3).»  £  span  (7(a).i  .  If  (7(a). <  and  C((3).i  are  on  the 
same  side  of  span  (7(a).f  =  span  C(/3).{  ,  then  det  (7(a)  and  det  C(ft)  are 
not  zero  and  have  the  same  sign.  Thus  (det  Mo,a)(det  Mpp)  <  0.  This 
is  impossible  when  M  €  Po  •  Thus,  (7(a).,  and  C(/?).,  are  on  opposite 
sides  of  span  C(a). j  .  As  in  the  proof  of  Lemma  2.18,  we  have 
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intpos  C(a). j  C  intfpos  C(a)  U  pos  C[p)\  C  int  K(M)  ,  so  span  C(a).$  can¬ 
not  be  a  supporting  hyperplane  of  K (M) . 

□ 

With  Lemma  2.18,  we  can  show  that  for  M  G  Qo  the  conditions  of 
Corollary  2.11  are  sufficient  as  well  as  necessary  for  M  6  Eq  not  to  be  in  U. 
We  have 

Theorem  2.20  If  M  G  Eq  n  Qo  fl  5?”Xn ,  then  M  £  U  if  and  only  if 
there  exist  a,p  6  (n)  and  i,j  G  n  such  that 

(i)  a^/3,  iy^j, 

(ii)  (det  Maa)(det  Mpp)  0  and  there  exists  a  nonzero  vector 
v  G  5R”  such  that  vTC(a).i  =  vTC(P).j  —  0, 

(iii)  C'(a).*  and  C(P).j  are  on  opposite  sides  of  span  C(a).{  =  span  C(P ).j 
[that  is,  (vT C{a).i)(v'rC(P).j)  <  0  |. 


Proof.  The  necessity  part  of  this  theorem  follows  from  Corollary  2.11.  Now 
suppose  that  the  conditions  are  satisfied.  We  know  that  pos  C(a)  is  a  full 
complementary  cone.  As  i  =£  j  and  span  C(a).s  =  spanC(/?).j  we  would 
have  C(a).j  C(P).i  and  C(P).i  G  span  C'(a).}  .  So  if  M  were  a  U-matrix, 
Lemma  2.18  would  imply  that  span  (7(a). {  is  a  supporting  hyperplane  of 
K(M) .  But  this  is  impossible  if  C{a).i  and  C[P).j  lie  on  opposite  sides  of 
span  C(a). j  . 

□ 


Condition  (iii)  in  Theorem  2.20  is  non-trivial.  Figure  2.4  shows  K(M) 
for  the  matrix 


M  = 


0 

1 
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which  is  in  Eq  f|  Qo  •  This  matrix  satisfies  (i)  and  (ii)  with  (a,  /3,  i,  j )  = 
(0,2, 1,2).  However,  (iii)  is  not  satisfied,  and,  indeed,  M  £  U. 


Notice  that  Theorem  2.20  implies  Example  2.12  must  have  used  a  matrix 
M  not  in  Qo  which,  in  fact,  it  did.  However,  Example  2.10  used  a  Qo  -matrix 
so  we  cannot  strengthen  Theorem  2.9  for  Eq  D  Qo  -matrices.  Example  2.14 
does  not  use  a  Qo -matrix,  but  we  still  cannot  strengthen  Corollary  2.13  as 
seen  by 

Example  2.21  Let 


M  = 


0  0  11 

.1  0  0  0  . 

The  full  complementary  cones  are  C(0),  C({3}),  C({1,4})  and 

C({1,2, 4}).  Suppose,  for  the  sake  of  contradiction,  that  M  j ?  E0  .  Then 
there  is  a  nonzero  x  >  0 ,  so  that  for  any  i  £  4 ,  if  x,  >  0 ,  then  (Mi),  <  0 . 
As  we  will  always  have  {Mx)^,{Mx)^  >  0,  we  must  have  X3  —  z<  =  0 .  But 
then  we  will  have  (M  1)2  >  0,  thus  requiring  that  12  =  0.  Cumulatively 
these  conditions  will  cause  (Mi)i  to  be  nonnegative,  leading  us  to  conclude 
that  =  0,  giving  a  contradiction.  Thus  M  £  Eo  •  M  has  three  non¬ 
trivial  principal  transforms  which  correspond  to  block  pivots  on  Maa  where 
a  can  be  {  3  }  ,  {  1,  4  } ,  or  {  1,  2, 4  } ,  and  are,  respectively, 


■  0 

—1 

-1 

0  • 

0 

0 

0 

1  ■ 

0 

0 

0 

1  ■ 

0 

0 

—1 

0 

1 

1 

0 

0 

—1 

1 

0 

0 

0 

0 

1 

—1 

f 

-1 

—1 

0 

0 

t 

0 

—1 

0 

0 

.  1 

0 

0 

0  . 

.  -1 

—  1 

-1 

0  . 

0 

—1 

—1 

0  . 

Similar  arguments  will  show  that  these  three  matrices  are  all  in  Eq  .  Thus 
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M  €  Eq  .  Now,  it  is  clear  that  K (M)  C  pos[  /  |  — M  j .  Given  the  following 
fact  (see,  for  example,  Proposition  4.2  of  Doverspike  and  Lemke  (1979)) 

M  £  Qo  <=>  pos(  1 1  — M  ]  =  [J  {  pos  C(a) :  det  C[a)  ^  0  } , 

and  noting  that 

pos[  1 1  —  M  ]  ==  {  x  6  9f4  :  zi,  x2,  X3  >  0  and  xx,  x2  >  —*3  } , 

pos  (7(0)  =  {  x  6  SR4  :  xx,  x2, 13,14  >  0  } , 
pos(7({3})  =  {x  6  SR4  :  xx,x2,  —13,14  >  0  and  xx,x2  >  — 13}, 
pos  C({  1, 4  })  =  {  x  6  SR4  :  xx,  x2,  — 14  >  0  and  x2  >  xx  >  — X3  } , 

pos  C({  1, 2, 4  })  =  {  x  €  !R4  :  xx,x2,  — X4  >  0  and  xx  >  x2  >  — 13  }  , 

we  find  that  M  €  Qo  .  As  in  Example  2.14,  for  x  =  (1,  2, 1)T  ,  we  have  that 

span  C({  1, 4  }). j  =  span  <7(0)4  and  C({  1, 4  }).j  x  =  C(0).4 1 ;  thus  M  ^  U . 

However,  there  are  only  four  candidiates  for  the  4' tuple  ( a,/3,i,j )  in  the 
conditions  of  Theorem  2.6,  and  checking  them  shows  that,  for  each,  we  have 
some  q  €  pos  C(a).j  \  pos  C(/3). j  and  some  g  G  pos  C(/?).j  \  pos  C(a).i , 

(C(0),  G({  1,  4  »,  4, 1)  9  =  (2, 1,  0,  Of  q  =  (1, 1,  -1,  Of 

(C(0),  C({  1,  2, 4  }),  4, 1)  q  =  (1,  2,  0,  Of  q  =  (1, 1,  -1,  Of 

(C({  3  }),  C({  1, 4  }),  4, 1)  q  =  (2, 1, 0,  Of  g  =  (1, 1, 1,  Of 

(C({  3  }),  C({  1, 2, 4  }),  4, 1)  g  =  (1,  2, 0,  Of  g  =  (1, 1, 1,  Of  . 

Hence,  M  is  an  example  showing  that  Corollary  2.13  cannot  be  strengthen 
to  say  that  its  three  conditions  are  necessary  for  a  matrix  in  Eq  fl  Qo  not 
to  be  in  U. 

We  now  come  to  a  result  which  says  that  when  K(M)  is  convex  and 
(g,  M)  has  a  unique  solution  for  all  g  €  int  K(M)  ,  the  matrix  M  cannot 
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have  any  negative  principal  minors.  The  proof  sheds  light  on  the  conical 
structure  of  Qo  f|U -matrices. 

Theorem  2.22  Qo  nU  C  P0  . 

Proof.  Let  ME  Qo  DUn9in><n.  There  are  two  cases.  If  M  E  Q,  then 
M  E  Qnu  =  P  C  Po-  Assume  therefore  that  M  E  Qo  \Q.  Thus, 
K(M)  7^9?n.  Suppose  we  have  a  collection  of  index  sets  ai,...,ak  €  (n) 
for  which 

detMajaj  >  0,  j  —  .  (2.9) 

We  know  fc  >  1  since  ai  =  0  belongs  to  the  collection.  Now  consider 

k 

Ck  =  U  pos  > 

j—i 


and  suppose  Ck  K(M)  . 

As  M  €  Qo  ,  K (M)  is  a  closed  convex  finite  cone.  The  cone  Ck  is 
closed  and  polyhedral;  by  our  assumption,  it  is  a  proper  subset  of  K(M) . 
Thus  there  must  exist  a  point  q  €  int  K(M)  \  Ck  •  Let  p  E  int  Ck  •  (Note: 
Ck  contains  and  so  has  a  nonempty  interior.)  Let 

L  =  {  r  :  r  =  (1  —  \)p  X?,  0  <  X  <  1 } , 

and  so  intL  =  {r  :  r  =  (1— X)p-f  \q,  0  <  X  <  1 }  .  As  L flint  K(M)  7^0, 
K (M)  is  convex,  and  (hence)  L  C  K(M) ,  we  have 

int  L  C  int  K( M)  and  int  L  (1  dCk  ^  0 . 


Thus, 


dCk  D int  K(M)  7^0. 
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Now  dCk  is  (n  —  1)  -dimensional  and  contained  in 


U  pos  C(<Xj).% . 

j€k 

i€n 

Hence  as  the  union  of  the  boundaries  of  the  pos  C(aj).i  is  (n — 2) -dimensional, 
and  as  we  can  slightly  perturb  the  position  of  the  point  p  we  selected  and 
still  keep  it  within  int  Ck ,  then  we  may  assume  that 

int  pos  B  fl  dCk  H  int  K(  M)  0  (2.10) 

where  B  =  C(aJ).A  €  sjjnxfn— i)  for  some  aj.  [n  the  given  set  satisfying 
(2.9).  Let 

Pj  =a,A{n}- 

Now  posfl  C  pos  C{(3j).  If  det  C(Pj)  =  0,  then  there  exists  a  point 
q  €  int  pos  C(Pj)0  int  K(M)  ,  and  (g,  M )  has  infinitely  many  solutions.  This 
contradicts  the  hypothesis  that  M  €  U.  Thus  det  C(Pj)  ^  0,  and  accord¬ 
ingly,  det M0Jp}  5^0.  If  det MpJp]  <  0,  then  as  det Maj0lj  >  0,  we  have 
(detC(aj))(detC(^))  >  0  implying  that  C[aj).n  and  C(Pj).n  lie  on  the 
same  side  of  span  2? .  Thus,  by  Lemma  2.8,  int  pos  C(aj)n  int  pos  C(Pj)  ^  0, 
which  contradicts  the  assumption  that  M  €  U.  Thus  det  Mp} >  0,  and 
we  have  I.n  and  —  M.n  lying  on  opposite  sides  of  spanB.  Hence 

int  pos  B  C  int[pos  C(aJ)  U  pos  C{Pj)\ . 

From  this  and  (2.10)  we  have  pj  g  {ai,...,a*}.  Let  ctfc+i  =  Pj  and 
adjoin  it  to  the  collection  of  known  index  sets  for  which  the  corresponding 
principal  minor  is  positive.  We  repeat  this  construction  until  l  index  sets  are 


found  and 


i 

C,  =  (J  pos  C(a})  —  K(M) . 
i-i 

If  P  6  (R)  and  p  f?  {  alt . . . ,  aj  } ,  then  det  Mpp  —  0 ;  otherwise 

int  pos  C(P)  0  int  K(M)  ^  $ 
and  this  implies  there  exists  a,  (1  <  j  <  l )  such  that 

int  pos  C(P)  H  int  pos  C(ctj)  ^  $ 
which  contradicts  our  assumption  that  M  6  U . 

□ 

The  next  theorem  sharpens  the  ideas  concerning  the  structure  of  d  K(M )  , 
for  ME  Qo  H  U ,  that  we  developed  in  the  proof  of  the  last  theorem. 

Theorem  2.23  If  M  €  (Qo  \  Q)nU0RnXn,  then  there  exists  a 
nonnegative  m  X  n  matrix  A  such  that 

K(Af)  —  {q  :  Aq  >  0} , 

and  the  number  m  is  minimal.  Moreover,  if 

a*  =  supp  Ak-  for  all  k  €  W 

then  det  Makak  =  0 .  If  det  Mpp  =  0,  for  some  p  E  (W) ,  then  there  exists 
k  €  WI  such  that  a*  C  P . 

Proof.  From  Theorem  2.22,  we  know  that  M  6  Po  \  P .  The  cone 
K(M)  being  convex  and  finitely  generated  can  be  expressed  as  a  polyhedral 
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convex  cone  (see  Weyl  (1935)).  Thus  there  exists  a  matrix  A  £  5RmXn  such 
that 

K  (M)  ={q-.Aq>  0}. 

The  matrix  A  can  be  chosen  so  that  none  of  its  rows  is  redundant.  Since 
5R!}_  C  K(M)  ,  it  follows  that  A  >  0  (and  A*.  ^  0  for  all  A:  £  m ).  Each  of 
the  hyperplanes 

H{Ak.)  =  {  x  €  0?n  :  A k,  x  =  0  }  A ;  £  m 

is  the  boundary  of  a  half-space 

H+(Ak.)  =  {xS9tn  :Alx>0}  kern 

and  has  an  (n  —  1) -dimensional  intersection  with  dK(M)  .  For  each  k  £  m, 
there  exists  an  a  €  (n)  such  that 

dim[pos  C(a)nH(Ak.)]  —  n  —  1 . 

If  det  Maa  =  0 ,  then 

pos  C{a)  Q  H(Ak.).  (2.11) 

If  det  Maa  then  by  (2.11)  there  must  exist  an  index  t  €  n  such  that 

dim[pos  C(a).j  fl  H{Ak.)\  =  n  — •  1 ,  (2-12) 

and 

C(a).itH(Ak.). 

Let  (3  —  a  &  { t } .  If  det  Mpp  ^  0 ,  as  M  £  Po  ,  we  have 

(det  C(a))(det  C(0))  <  0. 
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So  C(a). i  and  C(/3).i  lie  on  opposite  sides  of  span  C(a). »  =  span  C(/3).t . 
Thus 

int  pos  C(a).\  C  int(pos(7(a)UposC(/?))  C  int  K (M)  , 

which  contradicts  (2.12).  So,  det  Mpp  =  0. 

Hence  for  every  k  6  m ,  there  exists  a  /3 k  €  {fn) ,  with  det  Mphpk  —  0 , 
and  such  that  (2.11)  holds  with  a  =  (3 k-  Then 

j  €  ak  =*  I.k£H(Aic.)  =»  ^  pos (7(/?fc)  =»  —  Mfc  €  //(.A*.), 

and 

j  €  otk  =>  /.fc€H(Afc.). 

Thus,  the  columns  of  C(«fc)  are  all  in  H(Ak.)  which  implies  detM0kaj|1  =  0. 
In  fact,  if  P  €  (**)  and  det  Mpp  —  0,  then  pos(7(/3)  C  dK(M) ,  so 
pos(7(/3)  C  H{Ak.)  for  some  fc,  and  as  above,  j  6  ak  implies  I.k  g  H(Ak.), 
so  C(/3).j  =  —  M.y  .  This  implies  a*  C  /?. 

□ 

We  now  examine  a  situation  which  could  be  viewed  as  a  partial  converse 
to  Theorem  2.22.  It  involves  matrices  belonging  to  a  special  subclass  of  Po . 
We  shall  show  that  these  matrices  belong  to  U  and  that  they  give  rise  to 
cones  K (M)  of  a  special  form.  To  this  end,  we  introduce 

Definition  2.24  If  M  6  Po  n!Rn><n,  then  M  €  Pi  if  and  only  if 
there  exists  a  unique  index  set  a  6  (n)  such  that  det  Maa  =  0 . 

Thus,  M  €  Pi  if  and  only  if  it  has  nonnegative  principal  minors 
precisely  one  of  which  is  zero.  A  Pi  -matrix  may  or  may  not  belong  to  Q. 
For  instance 

1  -1 

M  =  €  Pi  nQ, 

1  0  J 
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tee  Figure  2.5,  whereat 


M 


€  Vi  \Q, 


tee  Example  2.2.  In  the  former  cate,  the  matrix  doet  not  belong  to  U,  but  in 
the  later  cate  it  doet. 


Theorem  2.25  If  M  e  (Pi  \Q)n3inXn,  then  M  €  U,  and  K(M) 
it  a  half-tpace.  Furthermore,  let  a  €  8?n  be  the  normal  to  the  hyperplane 
dK(M)  .  If  detMaa  —  0,  then  a  can  be  choten  to  that  aa  >  0  and 
114  =  0. 

Proof.  Let  K  t{M)  be  the  union  of  the  full  complementary  conet  attociated 
with  M .  Then  3K  r(M)  it  contained  in  the  union  of  the  boundaries  of  the 
full  cones.  Suppose  pot  C(a)  it  a  full  complementary  cone  and 


dim[0Kf(M)  n  posC(a).|J  =  n  —  1 . 


Let  j3  as  a  a  { i } .  If  pot  C{f3)  it  a  full  cone,  we  may  ask:  where  it 
C{0)  i  with  respect  to  span  C(a).j  =  span  C(0). j  ?  If  is  on  the  tame 

side  of  tpanC(a).i  as  C(a).<,  then  (det  C(a))(det  C(p))  >  0  giving  us  the 
contradiction  that  (det  Maa)(det  Mpp)  <  0.  If  C(/3).j  is  on  the  opposite 
side,  then 

int  pot  C(a).j  C  int(pos  C(a)  U  pot  C(/3) J , 


to 

dim|dKf(M)  D  pot  C(a).  1}  <  n  —  2 , 

a  contradiction.  Hence  pot  C(0)  it  degenerate  and  contains  C(a).|  .  There* 
fore  dK  f(M)  is  contained  in  the  union  of  the  degenerate  complementary 
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cones.  But,  by  hypothesis,  there  is  only  one  degenerate  complementary  cone. 
Since  M  <£  Q>  we  have  dKf(M)  ^0.  Thus,  dKf(M)  is  contained  in  this 
one  degenerate  complementary  cone. 

Let  L  =  {x  :  aTx  =  0}  be  the  affine  hull  of  this  degenerate  com¬ 
plementary  cone.  (Both  are  (n  —  1) -dimensional.)  Being  the  boundary  of 
an  n-dimensional  polyhedral  cone  contained  in  L,  dK.((M )  cannot  have  a 
boundary  relative  to  L.  Hence  dKf(JVf)  =  L,  and  K(Mj  is  a  half-space 
{ x  :  dTx  >  0}  with  0t^O'>'0  as  3?!).  C  K(A/)  . 

If  det  Maa  =  0 ,  then  pos  C(a)  is  the  only  degenerate  complementary 
cone.  Thus  Li  £  L  if  and  only  if  i  €  a.  This  implies  aa  >  0  and  a&  =  0 . 

Moreover  M  6  Po  C  Eq  ,  and  the  fact  that  the  only  degenerate  cone 
is  dKf(M)  forces  the  three  conditions  in  Theorem  2.6  to  fail  to  be  satisfied, 
so  we  have  M  6  U . 

As  final  remark  before  leaving  this  chapter,  lest  the  impression  be  given 
that  Eq  is  made  up  of  only  matrices  that  are  Po ,  U,  or  Qo ,  wc  give  an 
example  of  a  matrix  that  is  in  Eq  \  ( Po  U  U  U  Qo )  • 

Example  2.26  Let 

0  0—10 

0  0  0  1 

M  = 

10  0  0 

.01  00 

Clearly  Mg  P0  as  det  Maa  <  0  for  a  =  {  2, 4  } .  M  has  exactly  four  non¬ 
degenerate  cones:  C(0) ,  C(n) ,  C({  1,3}),  and  (7({  2,  4  }) .  Each  is  a  different 
orthant  in  !R4 ,  so  the  interiors  of  these  four  cones  are  pair-wise  disjoint,  and 


hence  M  €  Eq  .  However,  with  ( a,(3,i,j )  =  (0,  {  1, 3  },  3, 1)  we  can  satisfy 
the  conditions  of  Theorem  2.6  -  in  fact,  pos  C(0).g  =  pos  C({  1, 3  }).£  -  and  so 
M0.  Finally,  we  have  (0,2,0,0)T  €  K (M)  and  (0,0,0,-2)T  €  K (M) , 
but  (0, 1, 0,  — 1)T  g  K(M)  ,  so  K(M)  is  not  convex.  Hence  M  £  Qo  .  Thus 
M  6  Eq  \  (Po  UUU  Qo)  as  claimed. 
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Figure  2.4 
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CHAPTER  3. 

INS-MATRICES:  CHARACTERIZATION  RESULTS 


3.1  Introduction  to  INS-matrices 

We  have  now  defined  and  studied  the  class  U  which  generalizes  the  class 
P.  We  are  led  to  wonder  about  possible  larger  classes  containing  U.  As  before, 
we  must  decide  what  properties  we  wish  this  larger  class  to  inherit  from  U 
and  what  properties  we  wish  to  relax.  The  one  essential  property  of  U  is  the 
uniqueness  of  the  solution  to  ( q ,  M)  where  q  is  in  the  interior  of  K(1W) . 
However,  the  main  properties  of  the  combinatorial  and  geometric  structure 
of  K(M) ,  that  is  peculiar  to  those  M  6  U ,  is  derived  more  from  having  the 
same  number  of  solutions  everywhere  within  the  interior  of  K(M)  than  from 
that  number  being,  in  particular,  one.  With  this  in  mind,  we  focus  attention 
on  understanding  this  structure.  We  have 

Definition  3.1  For  any  fc  £  Z+,  a  matrix  A  is  said  to  be  an 
INS  it  -matrix,  A  £  INS* ,  if  and  only  if 


A  £  (J  {  M  €  9?nXn  :  |sol(9,  M)j  =  k,  for  all  q  £  int  K(M)  } . 

n 
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Definition  3.2  A  matrix  A  is  said  to  be  an  INS-matrix,  A  £  INS 
(Invariant  Number  of  Solutions),  if  and  only  if 

A£  \J  INS*. 
k€Z+ 

As  before,  and  as  will  be  shown  in  Theorem  4.7,  we  must  define  these  classes 
with  respect  to  q  in  the  interior  of  K (M) ,  not  all  of  K(M) ,  otherwise  these 
classes  will  contain  only  the  P-matrices.  Notice  that  we  have 


U  =  INSx  C  INS  . 


Thus  the  INS-matrices  seem  a  natural  extention  of  the  U-matrices,  but  are 
strictly  larger  as  seen  by 

Example  3.3  Let 


L  — 2  1  J 


As  illustrated  in  Figure  3.1,  M  £  INS2 .  Notice  that  the  full  complementary 
cones  can  be  partitioned  into  two  groups,  {  C({  1,2})}  and  {  C(0),  C{{  2})}, 
such  that  the  union  of  the  cones  in  each  group  covers  the  interior  of  K(M), 
and  the  interiors  of  the  cones  in  each  group  are  pairwise  disjoint.  We  also  see 
that  |sol(g,  M)f  for  q  £  dK (M)  is  one  or  infinity  -  never  two  -  for  points 
in,  respectively,  posC(n).2  and  int  pos  C(n).i . 

In  the  last  chapter  we  noticed  that  U  D  Q  =  P .  A  result  of  Murty’s 
shows  that  a  similar  result  holds  for  the  class  INS. 

Theorem  3.4  INS  fi  Q  =  P . 

Proof.  If  M  £  INSnQ,  then  int  K(M)  =  3Jn,  so  |sol(g,  M)\  is  constant 
for  all  q  £  3?n .  Theorem  7.10  from  Murty  (1972)  states  that  this  constant  is 
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equal  to  one.  Hence  M  E  U,  and  we  have  ME  P  as  desired. 

□ 

Before  continuing  on  to  the  next  sections,  where  we  look  at  what  goes 
into  making  an  INS-matrix,  there  are  a  few  concepts  which  should  be  brought 
up  first. 

Definition  3.5  Let  M  E  3?nXn,  we  then  define 
K(M)  =  (J  pos  C(a).{ . 

»e(n> 

ten 

K(M )  is  the  union  of  the  faces  of  the  complementary  cones.  It  contains,  in 
some  cases  equals,  3K(M)  .  In  Example  3.3,  (1,0)T  £  K{M)  \9K(Af)  , 
while  with  M  —  0  we  have  K(M)  =  dK (M)  .  K(M)  is  -the  set  of  all 
q  E  9?n  that  are  degenerate  with  respect  to  M .  Being  the  union  of  a  finite 
collection  of  sets  with  dimension  n — 1  or  less,  K[M)  has  zero  n-dimensional 
volume.  It  is  a  closed  cone  in  3in . 

We  will  be  interested  in  the  open  set  9?n\  K(M)  .  Let  £  be  the  collection 
of  the  connected  components  of  3?n\  K{M)  .  As  is  locally  path  connected 
and  as  3Rn  \  K(M)  is  open,  the  path  components  of  9?n  \  K{M)  are  the 
same  as  the  (connected)  components.  See,  for  example,  Munkres  (1975).  £ 
“almost”  partitions  SRn,  in  that  it  partitions  3Jn  \  K(M)  which  is  “almost” 
3in .  If  T  E  £ ,  then  T  is  an  open  polyhedral  cone,  i.e.,  dr  is  a  finite  collection 
of  (n —  1) -dimensional  finite  cones.  It  is  not  necessarily  true  that  T  =  intT, 
although  it  will  be  shown  later  that  T  C  int  T .  For  example 

Example  3.6  Let 

'  0  0  0 

M  =  10  0 

.000 
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then  £  contains  fcv  o  components:  Ti  =  int  and 


T2  =  int[SR3  \  SR?|_]  \  { i  €  SR3  :  £i  —  0, 13  >  0  }  . 

r2  7^intF2  =  int(S?3  \  9J^_].  This  also  shows  that,  if  T  6  £,  then  T  is  not 
necessarily  convex.  For  another  example  of  this 

Example  3.7  Let 

‘ 0  0  —  1  ' 

M=  0  0—1 
.0  0  — 1 . 

Then  £  contains  three  components: 
r\  =  int(3J3  \  JR3.] , 

r2  =  int  pos  C({  3  })  =  {1  €  5R3  :  ii  >  x3,  x2  >  x3,  x3  >  0  } , 
r 3  =  int  \  pos  C({  3  }) . 

Here  only  T2  is  convex,  although  I\  =  intTi  for  i  —  1, 2, 3 .  We  will  return 
later  to  the  subject  of  convexity  and  the  . 

We  now  discuss  necessary  conditions  for  a  matrix  to  be  INS. 

3.2  Necessary  Conditions  for  INS*matrices 

In  the  last  section  we  introduced  the  partition  of  3?n  \  K(M)  by  open 
polyhedral  cones  T  E  £ .  The  importance  of  this  structure  is  contained  in 

Theorem  3.8  If  T  €  £ ,  and  q,  q  €  T ,  then 
1 8°1  (?,  M)|  =  |sol(ijj,  M)|. 
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Proof.  Fix  q,q  £  T  £  E.  As  q,q  £  K (M)  ,  we  know  that  q  and  q  are 
not  contained  in  any  degenerate  complementary  cone,  and  are  not  contained 
in  the  boundary  of  any  nondegenerate  cone.  From  Chapter  1,  we  know  that 
any  solution  to  (q,  M)  is  associated  with  a  complementary  cone  containing 
q .  We  also  know  that,  if  the  cone  is  nondegenerate,  there  is  only  one  solution 
associated  with  it.  Now  if  q  £  pos  C(a)  then  q  £  int  pos  C(a ) .  Letting 
x  =  C(a)~1  q  >  0,  the  solution  associated  with  this  cone  is  (w,  z) ,  where 
za  —  xa  >  0  and  W&  =  >  0 .  As  in  Lemma  2.7,  any  other  solution  (w,  z) 

is  associated  with  another  complementary  cone  pos  C[/3)  containing  q.  Also, 
any  other  complementary  cone  containing  q  is  associated  with  a  different 
solution.  We  therefore  see  that  |sol(g,  M)\  is  the  number  of  complementary 
cones  that  contains  q .  The  same  holds  for  q . 

Suppose  that  for  some  a  £  (n)  we  have  q  £  pos  C(a)  and  q  £  pos  C(ot) . 
Then  any  path  from  q  to  q  must  contain  a  point  in  d  pos  C(a)  C  K(M)  ,  so 
q  and  q  are  not  in  the  same  path  component  of  3?"  \  K(M)  ,  i.e.,  not  in  the 
same  T ,  a  contradiction.  Thus  any  complementary  cone  containing  q  con¬ 
tains  q ,  and  vice  versa.  Thus  they  are  in  the  same  number  of  complementary 
cones,  and  so  |sol(g,  M)\  =  |sol(<j,  JVf)| . 

□ 

The  proof  just  given  shows  that  for  any  complementary  cone,  pos  C(a) , 
and  any  T  £  E, 


r  fl  pos  C(a)  7^  0  <=>  TC  pos  C(a) . 


The  main  result  of  this  section  is 
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Theorem  3.9  If  M  G  INS ,  then 


3K(M)  =  (J  (poaC(a).{  :  pos  C(a).i  is  not  proper}. 

«€(») 

•  €** 


Proof.  Let  pos  C(a)  be  a  degenerate  cone.  Suppose  pos  C(a)  fl  int  K(M)  0 . 
Then  there  will  exist  a  q  such  that 


q  G  int  pos  C(a)  fl  int  K(M)  . 


From  Proposition  1.6,  we  know  that  |sol(g,  M)  |  =  00 .  As  MG  INS ,  there 
must  be  infinitely  many  solutions  for  each  point  in  the  interior  of  K(M) . 
From  the  proof  of  Theorem  3.8,  we  see  that  for  any  point  in  5Jfn  \  K{M)  the 
number  of  solutions  it  has  to  the  LCP  is  equal  to  the  number  of  complemen¬ 
tary  cones  containing  it,  which  is  finite.  Hence  int  K(M)  C  K(M)  ,  but  this 
is  impossible  as  the  set  on  the  left  is  n -dimensional  and  the  set  on  the  right 
is  (n —  1) -dimensional.  Thus  all  degenerate  cones  are  contained  in  d¥L (M)  . 
(This  also  shows  that  INSqo  =  0,  so  our  definitions  cover  just  what  we  want 
without  any  technical  problems.) 

Suppose  now  that  pos  C(ar)  is  a  full  cone,  pos  C,(a).j  is  a  reflecting  face, 
and  pos  C^a).!  D  int  K(M)  7^  0 .  Then  there  is  a  q  G  int  pos  C(a). f  fl  int  K(M) 
such  that  for  any  /3  G  (n),  j  €.n,  we  have 


q  G  pos  C(/?).3 


dim[posC(/?).j]  =  n  —  1 
q  G  int  pos  C  span  C'(a).} 


and  any  small  enough  open  ball  around  q  is  bisected  by  int  pos  C(a).{  with 
the  two  open  half-balls  contained  in  rcTi  G  £,  respectively.  (We  are  not 


r 
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assuming  To  7^  I\ .)  Refer  to  Figure  3.2.  To  see  this  more  clearly,  notice 
that  the  set  of  points  that  are  in  either 

(i)  the  boundary  of  an  (n  —  1) -dimensional  face  of  a  complementary  cone, 

(ii)  a  A; -dimensional  complementary  cone  where  k  <  n  —  1, 

(iii)  the  intersection  of  pos  C(a).i  with  an  (n — 1) -dimensional  face 
(of  a  complementary  cone)  not  in  span  C(a).t  , 

is  a  set  of  dimension  less  than  n  —  1,  while  dim  [int  pos  C(a).c]  =  n  —  1. 
Furthermore,  as  all  the  A: -dimensional  facets  of  K(M)  are  closed  and  Unite 
in  number,  we  know  that  for  an  open  ball  around  q ,  that  has  a  small  enough 
radius,  we  will  have  a  A: -dimensional  facet  of  K(M)  intersecting  the  open 
ball  if  and  only  if  that  facet  contains  q . 

Since  pos  C{a). *  is  a  face  of  the  full  complementary  cone  pos  C(a) ,  then 
either  r0  H  pos  C(a)  7^  0  or  Ti  D  pos  C(a)  7^  0 ,  but  not  both  as  pos  C(a)  lies 
entirely  on  one  side  of  pos  C(a).j  .  Thus  without  loss  of  generality  we  assume 

To  C  pos  C(a)  and  Tx  D  pos  C(a)  =  0 . 

(Thus,  indeed,  Tot^Fi  .)  Let  Ho  and  Hi  be  the  two  closed  half-spaces 
with  spanC(a).j  as  boundary,  where  Tq  C  Hq  and  Ti  C  Hi  .  Suppose 
that  there  is  some  complementary  cone,  pos  C(0) ,  that  contains  T 1  but 
not  To .  Then  it  must  be  a  full  cone  and  have  some  face,  say  pos  C(P).j , 
containing  q .  By  (3.1)  this  face  lies  in  span  C(a).{ ,  hence  lies  in  Hi . 

However,  as  ppsC(a).j  is  reflecting  we  have  both  I\  and  — M.%  in  int/foi 
a  contradiction.  Thus  no  complementary  cone  contains  F 1  and  not  To .  But 
pos  C(a)  contains  To  and  not  Ti .  Hence, 

IsoK?1,  M) |  <  |sol(<?°,  M) |  +  1  q°  €  To,  ql  6  Tx .  (3.2) 
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Since  q  €  int  K(M) ,  we  have  Isol^1,  M)\  >  0,  so  ToUTi  C  int  K(Af)  . 
Hence  (3.2)  implies  that  M  £  INS ,  a  contradiction,  ^'hus 

M  6  INS  =>  dK(M)  D  { pos C(a).f :  pos C(a).{  is  not  proper). 

•>€(*) 

»€» 

Now  suppose  that  q  €  dK(M)  .  Clearly  q  is  not  interior  to  any  full 
cone.  Suppose  that  it  is  not  contained  in  a  degenerate  cone.  Then  it  is  on 
the  boundary  of  some  full  cone,  hence  q  €  dint  K(M)  .  As  int  K(M)  is 
an  n -dimensional  polyhedral  cone,  dint  K(M)  is  the  union  of  finitely  many 
(n  —  1) -dimensional  finite  cones,  each  contained  in  some  degenerate  cone  or 
a  face  of  a  full  cone.  If  pos  C(a).j  is  a  proper  face,  then  we  know  I.i  and 
— M.i  are  on  opposite  sides  of  span  C(a).i  .  Thus 

int  pos  C(a). j  C  intfpos  C(a )  U  pos  C(a  a  { t })]  C  int  K(M)  , 

giving 

dimfpos  (7(a).i  f|  d  K( M)  ]  <  n  —  1 . 

Thus  posC(a).i  is  not  a  face  containing  one  of  the  (n  —  1) -dimensional  finite 
cones  of  dint  K(M)  .  Thus  dint  K (M)  is  contained  in  the  reflecting  faces 
and  the  degenerate  cones,  and,  hence,  so  is  dK(Af) . 

□ 

Corollary  3.10  Let  M  £  SRnXn,  then 

dK(Af)  C  |J  {posC(a),i  :posC'(a).j  is  not  proper). 

«€<K) 

*€K 
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Proof.  Simply  notice  that  in  the  last  part  of  the  proof  of  Theorem  3.9  we 
never  used  the  fact  that  M  E  INS  when  showing  this  result. 

□ 

Saigal  (1972b)  uses  the  concept  of  a  “regular  pseudomanifold.”  We 
borrow  the  terminology  for  the  similar,  but  stronger,  concept  embodied  in 

Definition  3.11  Let  M  E  3fcn><n,  then  K(M)  is  said  to  be  regular  if 
and  only  if 

dK(M)  —  U  {posC(a).i  :  posC(a).{  is  no t  proper}. 

*€« 

Theorem  3.9  then  says  that 

M  E  INS  =4  K(M)  is  regular. 

This  is  the  general  necessary  condition  for  a  matrix  to  be  in  INS.  In  the  next 
section  we  take  up  the  question  of  this  condition’s  sufficiency. 


3.3  Sufficient  Conditions  for  INS-matriees 

We  now  know  that  if  a  matrix  M  is  in  INS  then  K(M)  is  regular.  The 
natural  question  is  to  ask  whether  this  is  a  sufficient  condition.  To  this  end, 
we  prove  the 

Lemma  3.12  Assume  M  €  3Jn><n  and  K(M)  is  regular.  Assume  also 
that  r0,r!  E  £  are  subsets  of  K(M)  —  and  hence  its  interior.  Suppose  for 
some  i  6  To  and  y  ET  i  there  is  a  path  L  €  int  K  (M)  from  i  to  y.  Then 
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L  can  be  chosen  to  have  the  following  “nondegeneracy”  properties: 

(i)  Lf)  K(M)  is  a  finite  set; 

(ii)  if  q  is  a  point  in  Lfl  K(M)  , 

then  q  is  in  the  interior  of  any  face  containing  it; 

(iii)  all  faces  containing  q  lie  in  the  same  hyperplane. 

Proof.  If  To  =  Ti ,  by  definition,  we  can  construct  a  path  L*  within 
To  from  x  to  y .  The  above  are  then  vacuously  true.  If  To  7^  Ti ,  we  can 
construct  the  path  L* -from  L  as  follows.  We  know  that  L  is  the  image  of 
some  continuous  function 

/:[o,i]-»»n,  /(o)  =  *er0,  /(ljf-yer!. 

Since  r0  is  closed,  we  have 

0  <  X  =  max{  /“ x(ro) }  <  1  • 

Let  q  —  /(X) .  Then  q  E  c)Fo .  Let  B  be  an  open  ball  in  int  K (M)  around 
q  E  L  C  int  K (M)  .  Since  all  the  facets  are  closed  sets,  we  may  assume  that 
B  is  so  small  that  any  facet  of  K(M),  of  any  dimension,  intersecting  B  must 
contain  q .  See  Figure  3.3  for  a  picture  of  the  local  situation  around  q . 

To  is  a  component  so  we  may  construct  a  path  L*  from  x  to  q  where 
L*\q  C  To  •  Let  q  E  B(l(L* \q).  We  claim  that  for  each  point  in  B\  K(M) 
there  is  a  path  in  B  from  that  point  to  q  that  satisfies  the  conditions  of  the 
lemma.  Clearly,  if  such  a  path  exists  between  q  and  some  point  in  £?nl\, 
then  one  exists  between  ?  and  all  points  in  BflT,  .  (This  does  not  follow 
from  what  has  been  set  up  as  it  could  be  that  B  0  Ti  is  not  connected.  In 
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this  case  we  may  temporarily  take  the  T,  as  the  connected  components  of 
B\  K(M)  and  all  will  go  through.  It  will  turn  out,  in  the  next  chapter,  that 
this  precaution  is  not  necessary.  However,  we  do  need  to  know  that  the  path 
can  be  built  within  B  for  later  reference.)  The  set  of  points  in  B  that  can  be 
connected  to  q  by  a  path  satisfying  the  given  conditions  is,  then,  the  closure 
of  the  union  of  some  of  the  B  f|  I\  .  Call  this  set  S .  S  is  the  intersection  of 
B  with  a  polyhedral  cone  with  vertex  translated  to  q.  It  is  n -dimensional 
as  Bn  To  C  S .  If  S  ?£B ,  then  S  has  a  boundary  in  B .  We  may  then 
find  a  point  q  £  B ,  in  the  interior  of  one  of  the  (n  —  1) -dimensional  faces 
making  up  dS ,  such  that  the  faces  of  K(M)  containing  q  all  lie  on  the  same 
hyperplane  and  all  contain  q  in  their  interiors.  (These  restrictions  will  remove 
a  set  of  points  that  is  (n  —  2) -dimensional  at  most,  and  we  have  a  set  that  is 
(n  —  1) -dimensional  from  which  to  choose.)  A  sufficiently  small  line  segment, 
L,  with  q  as  midpoint  and  orthogonal  to  the  (unique)  boundary  face  of  S 
through  q ,  will  make  a  path  from  some  tq  6  int S  to  some  ri  6  B\5  where 

Ln  K{M)  =q. 

The  conditions  of  the  lemma  are  satisfied  for  this  path.  Since  r0  6  S,  we 
have  a  path  to  ro  from  q  satisfying  the  conditions.  Combining  the  paths 
gives  a  path  from  q  to  ri  G  B  \  S  satisfying  the  conditions,  a  contradiction. 
Thus  B  =  S . 

Now,  let 

X'  =  max{  /-l(I\)  :flnr^0}. 

Clearly  X  <  X' ,  as  L  did  not  end  at  q .  Let 

?'  =  /(x')er2. 
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Pick  a  point  r  in  BnTa.  There  will  exist  a  path,  from  x  to  q,  from  q 
to  r,  and  from  r  to  q' ,  satisfying  the  conditions  of  the  lemma.  If  X'  =  1 , 
then  q'  =  y  and  we  are  done.  If  X'  <  1,  an  open  ball  around  q'  can  be 
made  small  enough  to  repeat  the  above  arguments,  extending  the  path  into 
some  new  component  T3 .  As  there  are  finitely  many  components,  we  will 
eventually  have  a  path  from  x  to  y  in  int  K(M)  that  satisfies  the  lemma’s 
conditions. 

□ 

Consider  a  point  q  6  L(~)  X’(JW)  .  The  previous  lemma  shows  that  for 
a  small  enough  open  ball,  B,  around  q,  there  is  a  hyperplane  H  such 
that  q  €  Bf)  K[M)  =  BC\H ,  and  B(~)H  splits  B  into  two  open  hemi- 
hyperspheres,  contained  in,  say,  and  1^  respectively.  (See  Figure  3.4.) 
Since  q  €  int  K(M)  ,  all  faces  containing  q  are  proper.  Suppose  that  a  full 
complementary  cone,  pos  C(a ) ,  contains  T 2  but  not  P 3  .  Hence  for  some 
t  €  n,  we  must  have  q  6  posC(a).j  C  H .  The  previous  lemma  allows 
us  to  assume  that  int  pos  C(a).j  bisects  B  into  the  aforementioned  hemi- 
hyperspheres.  As  posC(a).j  is  proper,  /.<  and  — M.i  lie  on  opposite  sides 
of  H.  Thus  pos C(a  A{i})  contains  Ts  but  not  I^.  Since  we  could  have 
assumed  at  the  start  that  pos  C(a)  contained  T 3  and  not  1^ ,  we  have  a 
bijective  correspondence  between  complementary  cones  containing  T2,  not 
r 3 ,  and  complementary  cones  containing  T 3 ,  not  T2 .  So  the  number  of 
complementary  cones  containing  T2  is  the  same  as  the  number  containing 
r3 .  Thus 

q  €  r2,  q  e  r3  =>  |sol(?,  M) I  =  jsol(g,  M)\ . 

Therefore,  if  we  start  at  x  and  follow  the  path  L,  we  will  pass  through  a 
finite  sequence  of  I\  €  E  where  \so\(q,  M)\  is  invariant  for  all  q  in  the  I\. 
Hence 
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|sol(i,  M) I  =  |sol(y,  M) I . 


We  have  been  assuming  that  x,y  fc  K(M)  .  Now  suppose  we  have 
y  E  int  K(M)  f)  K{M) .  As  in  the  proof  of  Lemma  3.12,  we  can  find  an  open 
ball  B  C  int  K(M) ,  with  y  E  B ,  so  small  that  B  n  K(M)  is  the  intersection 
of  B  with  the  union  of  finitely  many  (n  —  1) -dimensional  finite  cones  with 
vertex  translated  to  y.  (See  Figure  3.5.)  Since  y  is  contained  in  only  full 
complementary  cones,  each  cone  containing  y  is  associated  with  exactly  one 
solution  in  sol(y,  M ) .  Suppose  that  Ti  C  pos  C{ci) .  Then  y  E  pos  C(a) 
and  let  the  associated  solution  be  ( w ,  z) .  We  will  show  that  no  other  cone 
containing  Fi  has  ( w,z )  as  the  associated  solution  to  (y,  M) . 

We  may  assume  that  a  =  0  as  wc  can  always  block  pivot  on  Maa  to 
get  the  principal  transform  M .  As  shown  in  Chapter  1,  the  cone  structure 
is  preserved  and  working  with  posC7^-(0)  is  equivalent  to  working  with 
posCjvr(ar).  Thus  (w,  z)  =  (y,  0).  We  may  assume  that  supptu  =  suppy  = 
n\k,  where  0  <  k  <  n.  Thus  a  full  cone,  pos C(/3),  has  (y, 0)  as  its 
associated  solution  to  (y,  M)  if  and  only  if  0  6  (fc) .  However,  for  all 
P  (E  ( k )  and  for  all  i  E  k,  we  have  y  6  pos  (?(/?).; .  Hence  pos C(fl). 5  is 
a  proper  face  of  K(M)  as  y  6  int  K(M)  .  Therefore,  if  /3, 7  6  ( k ) ,  then 
(det  det  M-,7)  >  0 .  As  0  6  [k) ,  we  then  see  that  M ^  €  P .  If  P  6  (k) , 
and  there  is  some  0  <  y  6  pos<7(/3),  then  posC(/3)^  is  associated  with  a 
solution  to  (yp  Mpg) .  As  M E  P  and  y^  >  0 ,  there  is  only  one  such 
solution  and  it  is  associated  with  only  the  positive  orthant.  Thus  P  —  0.  Yet 
Pi  C  int  pos  C(0)  =  int  3i^  .  We  may  conclude,  as  claimed,  that  pos  C[a ) 
is  the  only  complementary  cone  containing  Ti  with  (tu,*)  as  the  associated 
solution  to  (y,  0).  Hence  jsol(y,  M)j  is  at  least  as  large  as  the  number  of 
complementary  cones  containing  Ti . 
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Suppose  now  that  y  £  pos  (7(a).  Thus  looking  at  Figure  3.5  again 
we  have  that  some  T *  containing  y  is  contained  in  pos  (7(a) .  In  fact, 
i?n  pos  (7(a)  is  the  closure  of  the  union  of  sets  in  the  form  i?nr\.  Select 
two  points  9  €  Tj  and  q  £  rfc .  Let  L  be  a  path  in  B  between  q  and 
q  satisfying  the  conditions  of  Lemma  3.12.  (In  the  proof  we  noted  that 
such  a  path  can  be  made  within  B .)  Suppose  we  cross  a  boundary  of 
pos  (7(a)  moving  along  L  from  q  to  q.  We  will  leave  the  cone  at  some 
point  interior  to  a  face,  say  pos  C(a). j .  This  face  must  be  proper,  as  it 
contains  a  point  in  int  K(M)  .  We  then  have  that  and  — M,  lie  on 
opposite  sides  of  span  (7(a). j .  Let  /3  =  ai{i}.  Thus  we  enter  pos  (7(/3) 
when  we  leave  pos  (7(a) .  Moreover,  L  C  B  so  pos  C(a).{  C\B  ^  0  implying 
y  £  pos  (7(a).i  =  pos  (7(/3).j .  Hence  the  solutions  to  ( y ,  M )  associated  with 
both  pos  (7(a)  and  pos  C(/3)  are  the  same,  both  using  only  the  columns  in 
pos  (7(a).{.  Thus  we  will  eventually  reach  a  full  cone  containing  Ti  such 
that  the  solution  to  ( y ,  M)  it  is  associated  with  and  the  solution  to  (y,  M) 
that  pos  (7(a)  is  associated  with  are  the  same.  Hence  |sol(y,  M) |  equals  the 
number  of  complementary  cones  containing  Ti .  We  have  thus  shown 

Theorem  3.13  Let  M  €  9?nXn.  If  K(M)  is  regular,  and  S  is  a(! 
connected  component  of  int  K(M)  ,  then 

q,qeS  =>  |soI(9,  M)\  =  |sol($,  M)\ . 


□ 


We  get  the  partial  converse  to  Theorem  3.9 

Corollary  3.14  Let  M  £  3RnXn.  If  K(M)  is  regular,  and  int  K(AT) 
is  connected,  then  M  £  INS  . 

□ 
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Example  2.3  shows  an  INS-matrix  for  which  int  K (M)  is  disconnected. 
In  fact,  there  exist  points  in  K(M),  for  example  (1,  l)r  and  (— 1, — l)r, 
which  can  be  connected  in  K(M)  only  with  paths  containing  the  origin.  Since 
K (M)  is  a  cone,  any  two  of  its  points  can  be  connected  by  a  path  through 
the  origin,  so  this  particular  K(M)  is  just  “barely”  connected.  However,  we 
note 

Theorem  3.15  Let  M  6  SRnXn,  n  >  1.  If  no  complementary  cone  is 
strongly  degenerate,  then  any  two  nonzero  points  in  K(M)  can  be  connected 
by  a  path  in  K(M)  not  containing  the  origin. 

Proof.  Define  the  map  F  :  3?n  — *•  as 


F(x)  =  y^(max(xt,  0)  •  I.i  -f-  min(xt,  0)  •  M.i)  . 

»=i 

Thus  K(M)  =  F(5Rn) .  Clearly  F  is  continuous.  Define  the  continuous 
radial  projection  P  :  \  {0}  — ►  Sn—1  as  P{x)  =  x/||x||,  where 

5n—1  =  {  x  €  3?"  :  ||x||  =  1 } 

is  the  unit  sphere  in  n -space.  Since  no  complementary  cone  is  strongly 
degenerate,  i^x)  =  0  implies  that  x  =  0.  So  0  £  F(Sn~l) ,  hence 
P  o  F  :  Sn~~l  -*  Sn~ 1  is  a  continuous  mapping.  Furthermore 

P  o  F{S n~1)  =  Sn— 1  H  K(M)  . 

This  and  the  path  connectedness  of  5n— 1  imply  that  5n—1  fl  K (M)  is  path 
connected.  However,  any  nonzero  point  in  K [M)  can  be  connected  by  a 
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path  to  Sn  1  D  K (M)  ,  i.e.,  the  ray  through  that  point  from  the  origin.  The 
theorem  follows. 

□ 

Example  3.16  The  matrix 

0  1  0  ' 

10  0 

—1  —1  —1  . 

belongs  to  INS2 .  However,  int  K(M)  is  not  connected,  and  no  complemen¬ 
tary  cone  is  strongly  degenerate.  This  example  shows  that  we  cannot  have  a 
result  similar  to  the  previous  one  concerning  the  connectedness  of  int  K(M) 
in  the  weakly  degeneracy  case.  However,  in  the  case  of  nondegeneracy  we 
have 

Theorem  3.17  Let  M  €  9JnXn  .  If  no  complementary  cone  is  degener¬ 
ate,  then  int  K(M)  is  connected. 

Proof.  Take  q,  q  E  int  K (M)  .  We  can  find  full  complementary  cones  so 
that  q  E  pos  (7(a)  and  q  E  pos  C(f3) .  If  a  =  /3 ,  then  q  and  q  are  path 
connected  within  int  pos  C(a)  C  int  K(M)  ,  even  though  q  and  q  may  be 
the  only  points  of  the  path  not  in  int  pos  (7(a) . 

Suppose  a  A  /?  =  { »}.  If  pos  (7(a).  1  ;s  reflecting,  then  /.*  and 
—  M.i  lie  on  the  same  side  of  span  (7(a). f .  By  Lemma  2.8, 
int  pos  (7(a)  ft  int  pos  C(/3)  7^  0 .  We  can  thus  build  a  path  in  int  pos  (7(a) 
from  q  to  a  point  in  this  intersection,  and  then  to  q  through  int  pos  C(/?) .  If 
pos(7(a).{  is  proper,  then  I.i  and  — M.i  lie  on  opposite  sides  of  spanC(a).|. 
So 

int  pos  (7(a). j  C  intfpos  C(a)  U  pos  (7(/3)]  C  int  K(M)  . 
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The  path  can  then  be  constructed  from  q  though  intposC(<a)  to  a  point 
within  int  pos  C(a).\ ,  and  from  there,  through  intpos  C(/3),  to  q. 


In  general,  if  a  A  (3  =  { i’i, . . .  } ,  let  71  =  a  and  for  1  <  j  <  k 

let  7y+i  =  7j  a  { ij  } .  Then,  posC'(7J)  and  pos  C(~fj+i)  are  adjacent  for 
1  <  j  <  k ;  moreover,  ^k+i  —  P  •  This  and  the  previous  arguments  show 
that  int  K(A/)  will  contain  a  path  from  9  to  9.  That  is,  int  K(M)  is  path 
connected. 

□ 

We  conclude  this  chapter  with  a  partial  characterization  of  the  class  INS. 

Corollary  3.18  Let  M  6  3?nXn,  and  suppose  that  M  has  no  zero 
principal  minors.  We  then  have 

M  €  INS  K(M)  is  regular. 

□ 
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CHAPTER  4. 

INS-MATRICES:  FURTHER  RESULTS 


4.1  Convexity  of  the  T 

The  partition  £ ,  defined  in  the  last  chapter,  was  seen  to  be  an  important 
object.  We  noted  in  Example  3.7  that  a  component  T  £  E  need  not  be 
convex,  even  if  T  €  K(M)  .  The  matrix  used  in  the  example  was  a  degenerate 
matrix,  but  degeneracy  was  unnecessary  as  the  matrix 


M  = 


0 

— 1 
0 


is  nondegenerate  and  has,  geometrically,  the  same  £  as  the  matrix  in 
Example  3.7.  However,  we  do  have  the  result 

Theorem  4.1  If  M  6  INS  n  9inXn ,  then  all  T  £  E  contained  in  K(Af) 
are  convex. 

Before  starting  the  proof,  we  will  need  the  following  lemma. 
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Lemma  4.2  If  M  6  9Jn><n  and  T^intT  for  some  T  €  E,  then 
T  C  int  T  and  x  £  int  T  \  T  implies  that  x  is  in  a  degenerate  cone. 

Proof.  T  is  open  and  contained  in  T.  As  intT  is  the  largest  open  set  in  T, 
it  follows  that  T  C  int  T . 

Now  I  €  r  \  r  =  ar  C  K(M )  .  Thus  x  is  contained  in  the  boundary 
of  some  complementary  cone,  say  pos  (7(a).  Suppose  pos  (7(a)  is  a  full  cone. 
Then  either  F  C  pos  (7(a) ,  or  T  D  pos  C(a )  —  0 . 

In  the  first  case,  x  £  <9[9?n  \  pos  C(a)]  C  !Rn  \  pos  C(a) .  Notice 
r  C  posCfa)  as  pos '(7(a)  is  closed.  Hence  9?n  \  pos (7(a) n int T  =  0,  a 
contradiction. 

In  the. second  case,  x  6  dpos<7(a).  As  pos  C(a)  is  full,  for  all  e  >  0, 
the  set  B(x,  e)  n  pos  C(a)  is  n -dimensional.  Now  3r  C  K{M)  is  (n  —  1)- 
dimensional  at  most,  so  B(x,  e)  \  T  0  for  all  e  >  0 .  Thus  x  G  \  T ,  and 
so  x  £  int  r ,  a  contradiction. 

We  have  shown  that  pos  (7(a)  is  a  degenerate  cone,  as  required. 

□ 

Proof  of  Theorem  4.1.  Suppose  there  exists  a  nonconvex  T  6  K (M)  .  Then 
there  exist  two  points  x,  y  £  T  such  that  the  line  segment  between  them, 

L  =  {  Xx  +  (1  —  \)y  :  0  <  X  <  1 }, 

is  not  contained  in  T.  Thus  there  must  exist  a  point  q  €  Lfldrn  K{M)  . 
We  may  assume 

q  =  Lf]  K{M)nB{q,e),  (4.1) 

for  some  small  e  >  0.  To  see  this,  notice  that,  for  small  e  >  0, 
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K(M)  f]B(q,t)  is  the  intersection  of  B(q,e )  with  a  finite  collection  of  finite 
cones  with  vertex  q  and  dimension  less  than  or  equal  to  n  —  1 .  Since  T  is 
open,  we  may  take  ex,  ey  >  0  small  enough  so  that 

Bx  —  B( x,  tx)  C  T  and  By  =  B(y,  ey)  C  T . 

We  may  thus  take  x  to  be  any  point  in  Bx  and  y  to  be  any  point  in  By . 
This  means  we  may  “perturb”  x  and  y ,  and  hence  the  line  segment  L ,  with 
n -dimensional  “freedom.”  We  can  thus  perturb  L  so  that  it  contains  q  and 
satisfies  (4.1).  See  Figure  4.1. 

For  the  moment  assume  that  q  £  int  K(iW)  .  Thus  q  is  not  in  any  de¬ 
generate  cone,  so  we  know  from  the  previous  lemma  that  q  £  T  n  int[5R”  \  F) . 
Thus,  for  all  e  >  0,  K(M)  C\B(q,  e)  must  be  ( n —  1) -dimensional.  Since 
we  can  perturb  L  with  n -dimensional  freedom,  we  may  assume  that  for 
q,  and  some  e  >  0  small  enough,  X"(M)  f|B(?,«)  —  H  r\B[q,e)  for  some 
hyperplane  H ,  see  Figure  4.2,  and  that  any  face  of  any  complementary  cone 
containing  q  is  (n  —  1) -dimensional  and  contains  q  in  its  interior.  (The 
argument  here  is  similar  to  several  given  before.  We  are  selecting  from  a  set 
that  is  (n  —  1) -dimensional  and  eliminating  a  set  that  is  at  most  (n  —  2)- 
dimensional.)  Now  let  pos  (7(a). j  be  a  face  containing  q .  q  £  int  K (M) 
implies  that  this  is  a  proper  face,  so  as .  q  6  dT  we  may  assume  T  C  pos  (7(a) , 
for  otherwise  T  C  pos  C(a  {  i  }) .  But  pos  C(a).i  C  H ,  so  pos  C(a ) ,  and 
hence  T ,  lies  entirely  on  one  side  of  H .  But  L  crosses  H  with  i  6  T  on 
one  side  and  y  £  T  on  the  other.  Contradiction. 

Now  assume  q  £  dK (M)  .  This  implies  q  £  3 int  K(M)  .  Again,  by  the 
perturbation  argument  given  above  and  the  fact  that  d  int  K(M)  is  a  finite 
set  of  (n  —  1) -dimensional  finite  cones,  we  can  assume  q  is  contained  in  the 
interior  of  some  face  pos  C(a).j  of  which  L  is  a  transversal.  As  q  is  in  3r 
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and  in  d  int  K(M)  ,  there  is  some  full  cone,  say  pos  C(a) ,  that  has  a  face  in 
the  affine  hull  of  pos  C{a).\  and  contains  T .  Thus  T  is,  again,  only  on  one 
side  of  the  affine  hull  of  pos  C(o).j .  Contradiction. 

□ 


As  a  side  result,  Lemma  4.2  implies 

Corollary  4.3  If  M  £  3tnXn  is  nondegenerate,  then  for  all  T  €  E 
we  have  F  =  int  T . 

□ 

We  remark  that,  even  for  nondegenerate  M  6  INS  ,  if  T  K(M)  then 
r  may  not  be  convex.  For  example,  in  !R2  if  we  let  M  =  — I  6  INS4  then 
we  get  |E|  =  2  where  one  component  is  int  K(M)  =  int  9?^  and  convex, 
with  the  other  component  being  5R2  \  3J?|_  and  nonconvex. 

Failing  to  show  for  nondegenerate  M  that  all  the  r  E  E  are  convex, 
one  might  consider  showing  that  some  particular  Y  is  convex.  With  this  in 
mind,  we  prove  the  next  theorem  before  leaving  this  section.  Recall  that,  by 
Theorem  3.8,  the  number  Isol^,  M)|  is  invariant  over  q  6  T  for  each  r  €  E . 

Theorem  4.4  Let  M  £  3inXn  be  nondegenerate.  There  exists  least 
one  T*  6  E  such  that  for  all  T  6  E 

|sol(g*,  M)1  >  jsol(g,  M)\ ,  for  q*  £  T* ,  q  €  T ,  (4.2) 

and  any  such  T*  must  be  convex. 

Proof.  It  is  clear  that  at  least  one  T*  exists.  As  in  the  proof  for  Theorem  4.1, 
we  assume  otherwise,  and  take  x,y  £  T*  such  that  the  line  segment  between 
them,  L,  contains  a  point  q  not  in  T* .  As  before,  using  the  nondegeneracy 
of  M ,  we  may  assume  q  £  dY*  and  that  there  is  a  hyperplane  H ,  of  which 


1 

«  • 
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L  is  a  transversal,  such  that  if  q  is  in  any  face  of  any  complementary  cone, 
then  q  is  in  the  interior  of  the  face,  the  face  is  (n  —  1)  -dimensional,  and  the 
face  is  contained  in  H .  We  also  know  at  least  one  face  of  a  complementary 
cone  contains  q .  If  a  complementary  cone,  with  a  face  containing  q ,  contains 
r* ,  then,  as  in  the  proof  of  Theorem  4.1,  we  will  have  T*  lying  entirely  on 
one  side  of  H .  As  before,  this  contradicts  the  fact  that  L  is  a  transversal 
of  H .  Thus  no  complementary  cone  with  a  face  contained  in  H  can  contain 
T* .  By  nondegeneracy  and  the  fact  that  some  face  does  contain  q  and  hence 
is  in  H ,  we  know  some  full  complementary  cone  does  have  a  face  lying  in  H . 
That  cone  must  contain  T ,  where  T  is  the  other  component  in  E  that  has  q 
on  its  boundary.  (Since  for  e  >  0  small  enough  we  know  that  B(q,e)\  K(M) 
is  two  hemi-hypersphercs,  one  on  each  side  of  H ,  we  see  that  at  most  two 
components  in  E  contain  q  on  their  boundaries.  We  know  q  €  dT*  and  we 
have  just  seen  that  another  component  must  also  have  q  on  its  boundary.) 
Hence,  every  complementary  cone  containing  T*  also  contains  T,  but  some 
cone  containing  T  does  not  contain  T* .  Thus,  with  q  6  T  and  q*  6  T* ,  we 
have  |sol(g,  M)\  >  |sol(g*,M)|.  This  contradicts  (4.2). 

□ 


4.2  The  Number  of  Solutions 

In  discussing  the  class  INS  an  important  question  to  ask  is  for  what 
values  of  k  is  INS*  empty?  We  know  INSi  =  U  is  certainly  nonempty.  It 
can  be  easily  seen  that  for  all  positive  integers  n,  — I  €  !RnXn  is  in  INS2«  . 
What  about  values  of  k  other  than  the  powers  of  two?  We  will  attempt  to 
give  evidence  suggesting  that  INS  *  =  (?  if  k  is  not  a  power  of  two.  We  begin 
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by  proving 


Theorem  4.5  Suppose  M  £  INS*  Pi  3iT‘Xn  .  If  there  exists  some  point 
in  dK(M)  that  is  not  contained  in  a  strongly  degenerate  cone,  then  k  is 
even. 

Proof.  Let  q  £  dK(M)  be  contained  only  in  full  or  weakly  degenerate 
cones.  By  dimensional  arguments  similar  to  ones  given  previously,  we  may 
assume  there  is  a  hyperplane  H  such  that  if  q  is  contained  in  a  face  of  a 
complementary  cone,  then  that  face  is  (n  —  1) -dimensional  with  q  in  its 
interior,  and  the  face  is  contained  in  H .  We  can  then  take  an  e  >  0  so 
small  that  B(q,  e)  fl  int  K (M)  C  T  for  some  particular  T  E  £  .  See  Figure 
4.3.  Any  full  complementary  cone  containing  q  must  contain  T,  and  likewise 
any  complementary  cone  containing  T  must  contain  q.-  Since  there  are  no 
strongly  degenerate  cones  containing  q,  by  Lemma  3.2  of  Saigal  (1972a)  it 
follows  that  q  is  contained  in  an  even  number  of  full  cones.  Thus  for  any 
q  E  r ,  we  have  |sol(^,  M)\  is  even,  whence  k  is  even. 

□ 

Corollary  4.6  Suppose  M  £  INS*n3inXn  .  If  there  are  no  strongly 
degenerate  cones  in  K(M),  then  k  is  even,  or  ME  P. 

□ 

We  now  reconsider  the  proof  of  Theorem  4.5.  This  time  we  will  allow 
strongly  degenerate  cones.  If  q  is  contained  in  a  degenerate  face,  then 
|sol(^,  M) |  =  oo .  Otherwise  q  is  contained  only  in  reflecting  faces,  as 
q  £  dK (M)  .  Thus  q  is  contained  only  in  full  cones.  Let  (w,  z)  £  sol(g,  M) 
be  the  solution  associated  with  a  full  cone  posC(a)  that  contains  q,  and 
so  there  is  an  i  E  n  such  that  q  £  intpos(7(a).i .  Thus  2a\{i}  >  0, 
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w&\{i}  >  0  and  Zi  —  wx  =  Hence,  (w,  z)  is  also  the  solution  associated 
with  the  full  cone  pos  C(a  A  {  i  ,  and  is  associated  with  no  other  full  cone. 
Thus,  as  we  had  q  contained  in  k  full  cones,  it  follows  that  |sol(9,  M) |  =  $ . 
In  any  case,  |sol(g,  M) j  7^  k .  This  reasoning,  along  with  Theorem  3.4,  proves 
the  following  assertion  which  was  mentioned  at  the  beginning  of  Chapter  3. 

Theorem  4.7 

p  =  U  U  {M  e  senxn  :  |sol(<7,M)|  =  A;,  for  all  <7  €  K(M)  } 

□  . 

At  the  start  of  this  section  it  was  suggested  that  INS*  =  0  if  k  is  not 
a  power  of  two.  As  will  be  shown  later,  this  would  follow  from 

Conjecture  4.8  Let  M  e  INSfc  nSRnXn .  If  K(M)  has  no  reflecting 
faces,  then  k  <  2. 

The  author  has  examined  many  INS-matrices,  and  studied  their  general 
structure  in  the  case  where  all  boundary  faces  are  degenerate.  No  coun¬ 
terexample  to  Conjecture  4.8  has  been  found.  To  obtain  some  feeling  for 
why  the  conjecture  should  be  true,  let  us  consider  trying  to  construct  K(M) 
for  an  INSfc  matrix,  k  >  3,  with  all  boundary  faces  degenerate.  Clearly 
dint  K(M)  0,  otherwise  ME  P.  Let  H  be  a  hyperplane,  let 
C  =  H  n  d  int  K(M)  ,  and  suppose  that  dim  C  =  n—  1 .  Since  only  degener¬ 
ate  faces  are  in  dK[M)  ,  each  such  face  acts  as  the  “base”  of  at  most  one  full 
complementary  cone.  We  would  then  find  that  every  point  in  C  that  is  not 
in  a  m -dimensional  facet  of  a  complementary  cone,  where  m  <  n  —  2,  i.e., 
“almost  all”  the  points  in  C ,  must  be  contained  in  exactly  k  degenerate  faces 
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which  act  as  bases  for  k  full  complementary  cones.  In  building  K (M)  we  find 
there  is  a  “tradeoff’  in  our  placement  of  the  column  vectors  of  [  I  |  — M  ] . 
The  more  we  place  in  C ,  the  more  degenerate  faces  we  will  have  to  form 
bases  of  full  cones,  which  can  be  used  for  this  multiple  covering  of  C .  The 
more  we  place  outside  of  C ,  the  more  full  cones  we  can  actually  form  on 
these  degenerate  faces.  There  are  other  tradeoffs  in  the  construction.  For 
instance,  the  more  boundary  hyperplanes  H  that  <9intK (M)  has,  i.e.,  the 
more  possible  C ’s  that  exist,  the  more  we  must  worry  about  putting  the 
column  vectors  of  [  I  |  — M  ]  on  the  boundary  of  each  C  to  “spread  them 
around”  to  the  different  C’s.  The  fewer  the  number  of  boundary  hyper¬ 
planes,  the  more  likely  the  C’s  will  contain  lower  dimensional  linear  spaces 
(linealities),  requiring  many  degenerate  faces  for  our  multiple  coverings  of 
the  C ’s,  and  the  previously  mentioned  tradeoff  becomes  more  critical.  With 
these  and  other  requirements  on  the  structure  of  K(M),  including  the  way 
in  which  the  columns  vectors  of  [  /  ]  —  M  j  form  the  complementary  cones, 
it  seems  certain  that  the  2 n  column  vectors  of  [  I  |  — M  ]  would  not  permit 
fc  to  exceed  two.  If  this  is  so,  we  have 

Theorem  4.9  If  M  €  INS^nS?71*71  and  Conjecture  4.8  is  true,  then  k 
is  a  power  of  two. 

Proof.  The  proof  uses  induction  on  n.  If  n  =  1,  then  there  are  at  most 
two  complementary  cones.  Thus  k  <  2  and  the  theorem  is  true.  Suppose 
the  theorem  is  true  for  n  —  1 .  If  no  faces  are  reflecting,  then  k  <  2  by 
Conjecture  4.8  and  the  theorem  holds. 

Thus  suppose  posC(a).{  is  a  reflecting  face.  Then,  posC(a).j  C 
3intK(M)  .  Let  H  be  the  hyperplane  spanC(a).j.  Let 
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s=  (J  {pos  :  pos  COS).,  C  H  > . 

*e(!T) 

pos  C(0). i  is  in  S  if  and  only  if  the  columns  of  C(/3). j  are  all  in  S .  (Notice 
that  we  use  •*  here  as  both  /.,•  and  — M.,-  are  on  the  same  side  of  H ,  but 
not  in  H .)  We  can  now  think  of  the  vectors  of  and  — M,  as  forming 
an  (n  —  1) -dimensional  LCP.  (For  notation,  say  that  the  matrix  associated 
with  this  new  LCP  is  M .)  The  correspondence  is  as  follows: 

H  takes  the  place  of  SRn— 1 ; 

S  takes  the  place  of  K(M); 

pos  C(a).j  takes  the  place  of  the  identity  matrix  as  pos  C(a).j  is  a  known 
full  cone  in  H ; 

if  pos  C[a).j  €  H ,  then  — M.j  is  represented  by  pos  C(a).j  ,  otherwise 
— M.j  is  represented  by  the  zero  vector.  (Here  we  index  on  t,  so  we 
have  j  €  » •) 

We  will  refer  to  this  LCP  in  if  as  the  reduced  LCP.  We  claim  that  K(M) 
is  regular. 

Suppose  q  €  int  S  is  contained  in  a  reflecting  face  of  the  reduced  LCP. 
By  dimensional  arguments  similar  to  earlier  ones,  we  may  assume  there 
is  an  (n  —  2) -dimensional  hyperplane  H  C  H  such  that  if  a  face  of  a 
complementary  cone  in  the  reduced  LCP  contains  q ,  that  face  is  (n  —  2)- 
dimensional,  contains  q  in  its  interior,  and  is  contained  in  H .  Thus  for 
e  >  0  small  enough,  B{q,  z)C\H  is  a  hypersphere  divided  into  two  hemi- 
hyperspheres  by  H ,  with  one  hemi-hypersphcre  contained  in  dr  and  the 
other  contained  in  dr' .  Here  T  and  T'  are  in  the  E  of  the  original  LCP. 
See  Figure  4.4.  Let  pos  C(a).t^  be  a  reflecting  face  in  S  containing  q .  Thus 
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both  full  cones  of  the  reduced  LCP  which  contain  that  face  lie  on  the  same 
side  of  H  in  H .  We  may  assume  they  both  contain  dFf~)H ,  and  both 
intersect  dV  C\H  only  on  H .  But  then,  as  both  I.j  and  — M.j  lie  in  H  on 
the  same  side  as  dFC\H ,  no  full  cone  of  S  can  contain  dV  OH  and  without 
containing  dF  f)H .  Thus  more  cones  of  the  reduced  LCP  contain  dT  (~)H 
than  contain  dF'  C\H .  But  each  full  cone  of  the  reduced  LCP  is  a  face  for 
exactly  two  cones  of  the  original  LCP.  (The  cones  you  get  by  adding  in  /.,• 
and  — M.i ,  respectively,  as  another  generator  of  the  cone.)  Also,  each  full 
cone  of  the  original  LCP  with  a  face  in  H  has  that  face  as  a  full  cone  of  the 
reduced  LCP.  Hence  as  q  £  int  S ,  we  have  some  full  cone  of  S  containing 
dF'C\H  and  so  if  iff  and  y  €  F'  then  |sol(x,  M)\  >  |sol(y,  M)|  >  0. 
This  contradicts  the  assumption  that  M  £  INS . 

Now  suppose  q  £  int  S  is  contained  in  a  degenerate  cone,  say  pos  Cj^[a) , 
of  the  reduced  LCP,  where  i  £  a  £  (n) .  If  all  the  columns  in  C^(a)  are  from 
the  original  LCP,  i.e.,  none  of  them  are  zero  columns  made,  as  mentioned 
before,  because  the  associated  — M.j  was  not  s7.  H ,  then  pos[  |  C^[a)  ] 
is  a  degenerate  cone  of  the  original  LCP.  What’s  more,  as  q  6  int  S ,  this 
degenerate  cone  contains  points  in  the  interior  of  the  convex  hull  of  S  and 
I.i ,  which,  in  turn,  is  contained  in  K(M) .  This  is  impossible  since  M  £  INS . 
We  thus  assume  Cj^(a)  contains  columns  which  were  made  zero,  as  described 
before.  Now  substitute  for  all  but  one  of  these  columns  that  were  made  zero, 
say  all  but  — M.j ,  the  associated  complementary  column  from  Let 

this  new  matrix  be  C^(/?),  and  notice  that  q  £  pos  CM-  if  ca(0*U}) 
is  a  degenerate  cone  in  the  reduced  LCP,  then  as  none  of  its  columns  were 
made  zero  in  the  way  initially  described,  we  would  be  back  to  the  previous 
case.  Thus  assume  that  pos  C^(/3  &■{]})  is  a  full  cone  in  the  reduced  LCP, 
and  thus  dim(pos  =  n  —  2 .  We  can  now  use  the  same  argument  as 
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in  the  case  when  we  assumed  q  was  in  a  reflecting  face  of  S .  However,  here 
we  use  pos  C^(/3).j  instead  of  pos  C(a)t-  .  Also,  there  is  one  full  cone,  not 
two,  of  the  reduced  LCP  with  posC^(/?)j  as  a  face,  but  we  still  have  this 
one  cone  containing  dr  C\H  -  and  not  containing  dr'  C\H  -  so  the  argument 
remains  valid.  We  conclude  that  S  =  K(M)  is  regular,  as  claimed. 

Let  q  be  in  a  connected  component  of  the  interior  of  S .  By  familiar 
dimensional  arguments,  we  may  assume  that  if  a  complementary  cone  of 
the  reduced  LCP  contains  q,  then  it  is  a  full  cone  containing  q  in  its 
interior.  Thus,  for  an  e  >  0  small  enough,  B(q,  e)  fl  int  K(M)  C  T  for 
some  particular  T  €  S.  (See  Figure  4.3  again.)  Since  q  £  dK(M)  ,  the 
complementary  cones  of  the  reduced  LCP  that  contain  q  are  reflecting  faces 
of  the  original  LCP.  (They  can’t  be  degenerate  faces  as  both  /.»  and  — Mj 
are  not  in  H .)  Thus  each  cone  of  the  reduced  problem  that  contains  q  is  the 
face  of  two  distinct  full  cones  of  the  original  problem,  and  these  two  cones 
will  contain  T .  Also,  any  cone  containing  T  must  contain  q .  As  we’ve  seen 
before,  the  number  of  cones  containing  T  must  be  k,  hence  the  number  of 
full  cones  in  the  reduced  LCP  containing  q  must  be  $ .  Since  q  could  be  in 
any  connected  component  of  int  S ,  using  Theorem  3.13  we  find  M  £  INSjt/j  . 
By  induction  on  the  dimension  of  the  LCP  we  see  that  $  is  a  power  of  two. 
Thus  A:  is  a  power  of  two. 


The  previous  theorem  makes  it  seem  almost  certain  that 

00 

INS  =  |J  INS  a* . 

P“0 

However,  there  is  a  large  class  of  matrices  for  which  we  can  show  the  result 
of  the  theorem  holds  without  recourse  to  Conjecture  4.8.  We  see  this  in  the 
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following 


Theorem  4.10  Let  M  e  INSfcn3inXn.  Suppose  that  for  all  a  6  ( ri ) 
we  have  det  C(a)  =  0  if  and  only  if  C(a).i  —  0  for  some  t  €  n.  It  is  then 
the  case  that  A;  is  a  power  of  two. 

Proof.  The  point  here  is  to  show  that  the  proof  of  Theorem  4.9  goes  through  - 
without  using  Conjecture  4.8  and  with  only  minor  changes  -  when  we  restrict 
ourselves  to  the  matrices  described  in  the  statement  of  Theorem  4.10.  (We 
use  here  the  notation  of  the  proof  of  Theorem  4.9.) 

In  the  case  where  we  have  a  reflecting  face,  the  proof  is  the  same.  The 
only  thing  needing  commentary  is  the  induction  step  where  we  must  now  show 
the  reduced  LCP  satisfies  the  hypothesis  of  this  theorem.  Suppose  pos  C^(a) 
is  a  degenerate  cone  in  S ,  where  t^a£(n).  Assume  no  column  of  CM 
is  zero.  Thus  all  the  columns  in  Cfa(a)  come  from  the  original  LCP,  i.e., 
are  not  “artifical”  zero  columns  as  described  before,  and  so  pos  C(a)  is  a 
degenerate  cone  of  the  original  LCP  with  no  zero  columns.  This  contradicts 
the  fact  that  the  original  problem  satisfied  the  hypothesis  of  the  theorem. 
Hence  the  reduced  problem  satisfies  the  hypothesis  of  the  theorem. 

Now  suppose  there  are  no  reflecting  faces.  If  M  €  Q,  then  M  €  P 
and  we’re  done.  Otherwise  Sint  K(M)  ^0  and  so  must  be  made  up  of 
degenerate  cones.  Thus  M  must  have  at  least  one  column  that  is  all  zeros, 
say  M.j  =  0 .  Thus 


{?€»$.:  ?,  =  0}  C  Sint  K(M)  , 


and  we  can  let 


f/  =  <9  €»":*  =  0}. 


We  can  now  go  through  with  the  proof  of  Theorem  4.9  for  the  case  of  a 
reflecting  face.  The  reduced  LCP  is  made  in  the  same  way.  J.|  represents 
C*(i)  taking  the  place  of  the.  identity  matrix  for  the  reduced  LCP.  For 
j  G  i,  if  — M.j  €  H  then  — M.j  represents  — M.j ,  otherwise  — M.j  =  0. 
The  difference  is  that  each  full  cone  in  S  is  the  face  of  one  full  cone  in  the 
original  LCP  -  which  will  contain  /.,•  -  and  for  each  full  cone  of  the  original 
LCP  with  a  face  in  H ,  that  face  will  be  a  full  cone  in  S .  We  will  finally  get 
the  reduced  LCP  in  INS* ,  which,  by  induction,  will  mean  k  is  a  power  of 
two.  (The  reduced  LCP  satisfies  the  hypothesis  of  this  theorem  by  the  same 
reasoning  as  given  in  the  second  paragraph  of  this  proof.)  We  thus  arrive  at 
the  same  conclusion  as  in  Theorem  4.9. 

□  • 

We  leave  this  section  with  the  following  immediate  corollary  to  the  last 
theorem. 

Corollary  4.11  If  M  €  INS*  fi  SRn><n  is  nondegenerate,  then  k  is  a 
power  of  two. 

□ 


4.3  The  Structure  of  K(M)  and  9K(M) 

The  purpose  of  this  section  is  to  build  a  link  between  the  combinatorial 
and  geometric  representations  of  K(M )  for  nondegenerate  INS -matrices.  The 
main  result  is  to  show  K (M)  and  dK(M)  can  be  divided  into  several 
disjoint  pseudomanifolds.  For  this  purpose  we  review  some  of  the  basic 
definitions  related  to  pseudomanifolds.  For  a  more  detailed  discussion  of 
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these  topological-combinatorial  constructs  discussion  see,  for  example,  Eaves 
(1972,  1976),  Freund  (1980),  and  Spanier  (1966). 

Definition  4.12  Let  V  be  a  finite,  non-empty  set  of  elements  (vertices). 
We  say  that  a  collection  P  of  subsets  of  V  is  an  n- dimensional  pseudo¬ 
manifold  if  and  only  if 

(i)  S  €  P  implies  that  |S|  =  n  -j-  1 .  The  subsets  S  are  referred  to  as 
n-simplexes. 

(ii)  F  C  V  and  |F|  =  n  implies  that  F  is  a  subset  of  at  most  two  elements 
in  P.  (F  is  an  (n  —  1) -simplex.) 

(iii)  For  every  pair  S,S  €  P ,  there  is  a  finite  sequence  S  =  So,  Si,... ,  Sm  = 
S  of  elements  of  P  such  that  (S*  ft  S.+il  =  n,  for  0  <  i  <  m. 

The  boundary,  OP ,  of  the  pseudomanifold  P  is  the  collection  of  subsets 
F  C  V  which  have  n  elements  and  are  contained  in  exactly  one  element  of 
P. 

Definition  4.13  Let  S  be  a  simplex  of  the  n -dimensional  pseudo¬ 
manifold  P.  Let  (so,  3i, . . . ,  an)  be  some  fixed  ordering  of  the  elements  of 
S .  Any  ordering  of  these  elements,  say  (sjB,  Sj, , . . . ,  sJn) ,  is  then  defined  to  be 
a  positive  (negative)  orientation  if  and  only  if  the  permutation  (jo,ji, ... ,  jn) 
is  even  (odd).  In  this  way  we  say  we  have  oriented  the  simplex  S .  We  say 
two  distinct  simplexes  in  P  are  adjacent  if  they  have  n  elements  in  com¬ 
mon.  Thus,  if  S  and  S  are  adjacent,  we  can  write  S  =  (s,si, . . . ,  sn)  and 
S  =  (5, si,..., sn).  If  these  particular  orderings  for  S  and  S  are  given 
different  signs  by  the  orientations  on  S  and  5 ,  then  we  say  S  and  S  are 
coherently  oriented.  Finally,  we  say  P  is  orientable  if  we  can  specify  an 
orientation  for  all  S  6  P  such  that  any  two  adjacent  simplexes  are  coherently 
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oriented. 

Example  4.14  For  any  matrix  M  €  5RnXn,  K (M)  can  be  viewed  as 
the  geometric  representation  of  an  orientable  (n  —  1) -dimensional  pseudo¬ 
manifold  without  boundary,  i.e.,  the  boundary  is  an  empty  set.  (Notice  the 
combinatorial  dimension  is  one  less  than  the  geometric  dimension.)  Let  V 
be  the  set  of  column  vectors  in  the  matrix  [  I  |  —  M  ] .  The  elements  of 
the  pseudomanifold  are  the  sets  of  column  vectors  of  the  complementary 
matrices.  The  geometric  representation  of  C(a).p  is  then  pos  C(a).p  for 
a,  f}  6  (n).  For  any  a  6  (n),  let  the  orientation  of  (C(a).i, . . . ,  C(a).n)  be 
determined  by  the  sign  of  ( — l)l“l .  It  is  not  hard  to  see  this  will  orient  the 
pseudomanifold. 

Doverspike  and  Lemke  (1981)  showed  that  for  a  large  class  of  non- 
degenerate  matrices  M  6  Qo  >  it  is  possible  to  find  a  collection  of  complemen¬ 
tary  cones  whose  union  is  K(lWj ,  and  forms  a  pseudomanifold  P  in  such  a 
way  that  the  geometric  union  of  the  faces  in  dP  is  d  K(M)  .  Furthermore, 
there  will  be  exactly  one  other  collection,  disjoint  from  the  first,  of  com¬ 
plementary  cones  whose  union  is  also  K(M) ,  which  also  is  a  pseudomanifold 
whose  boundary  is  dP .  As  we  will  be  building  somewhat  similar  pseudo¬ 
manifolds  from  INS-matrices,  eventually  to  prove  Theorem  4.18  -  which  the 
reader  may  wish  to  glance  at  now  -  it  will  be  useful  at  this  point  to  go  over 
the  proof  of  the  Doverspike-Lemke  result  before  proceeding.  The  basic  idea 
of  the  proof  is  explained  in  the  following  paragraph.  (The  full  details  of  the 
proof  would  require  many  pages  and  are  omitted.) 

Consider  the  geometric  structure  of  K(M) .  For  each  1-dimensional  facet 
of  K(M)  we  find  a  column  from  [  /  |  — M  ]  whose  “pos”  spans  it.  (We  have 
our  choice  of  any  column  vector  which  is  in  the  facet  when  there  is  more  than 
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one.)  In  this  way  we  build  up  (trivial)  pseudomanifolds  for  the  1-dimensional 
facets  of  K(M) .  From  here  on,  we  assume  we  have  built  up  pseudomanifolds 
for  the  r  -dimensional  facets  of  K(M),  1  <  r  <  n.  The  boundary  of 
any  (r  -|-  1) -dimensional  facet  is  the  union  of  r -dimensional  facets.  We  can 
take  the  union  of  the  pseudomanifolds  of  these  r -dimensional  facets  as  a 
boundaryless  pseudomanifold  over  the  geometric  boundary  of  our  selected 
(r  -f- 1) -dimensional  facet.  (They  will  “fit”  together  as  their  boundaries  were 
made  from  the  same  pseudomanifolds  on  the  (r  —  1) -dimensional  facets.)  We 
then  give  a  construction  to  show  there  will  be  exactly  two  pseudomanifolds, 
as  previously  described,  on  the  (r  -f-  1) -dimensional  facet  whose  boundary 
pseudomanifold  is  the  pseudomanifold  we  pieced  together  on  the  geometric 
boundary  of  the  (r  +  l) -dimensional  facet.  We  continue  this  until  r-f-1  =  n , 
at  which  point  we  have  the  result. 

The  concept  we  wish  to  use  from  this  is  the  family  of  pseudomanifolds  on 
the  r -dimensional  facets  of  K(M) ,  with  the  r -dimensional  pseudomanifolds 
forming  the  boundaries  of  the  (r-j-1) -dimensional  pseudomanifolds.  However, 
we  will  be  working  from  the  higher  dimensions  to  the  lower  dimensions, 
whereas  Doverspike  and  Lemke  do  the  opposite.  Notice  we  are  able  to  start 
our  constructions  since  K(M)  is  regular,  which  implies  that  for  any  face  in 
K(M),  say  pos  C(a).j,  we  have 

dim[posC(a).}n^K(M))  =  n  —  1  =>  pos  C(a).t  C  dK(M ) . 

The  following  lemma  will  prove  useful. 

Lemma  4.15  Let  M  €  INS  f)  9Rn><n  be  nondegenerate.  It  is  then  the 
case  that  the  r -dimensional  facets  of  K (M)  are  regular.  (That  is,  if  the 
(r— 1)  -dimensional  cone  pos  C(a).p ,  where  a,  /?  €  (n)  and  |/?|  =  r— 1,  is  the 
common  face  of  two  r -dimensional  complementary  cones  in  an  r  -dimensional 
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facet  of  K(M),  then  it  is  either  a  boundary  face  of  the  facet,  or  it  is  a  proper 
face.  If  the  two  r -dimensional  cones  are  not  in  the  same  r -dimensional 
facet  of  K(M)  then,  clearly,  pos  C(a).p  is  on  the  common  boundary  of  both 
r  -dimensional  facets  which  contain  the  r  -dimensional  cones,  so  we  need  not 
worry  about  this  case.) 

Proof.  This  is  easily  seen  by  reverse  induction.  It  is  true  for  dimension 
n,  as  K(M)  is  regular  by  assumption.  Suppose  it  is  true  for  dimension 
r  1,  1  <  r  <  n .  Suppose  that  it  fails  in  dimension  r.  We  may  assume 
some  q  in  the  interior  of  an  r -dimensional  facet  is  contained  in  a  reflecting 
(r  —  1) -dimensional  face,  pos  ,  which  is  the  common  face  of  the 

two  cones  pos[  |  — M.r  ]  and  pos{  |  I.T  ]  contained 

in  the  r -dimensional  facet.  (As  M  is  nondegenerate,  there  cannot  be  any 
degenerate  faces  here.)  Some  (r  +  1) -dimensional  facet  will  contain  this 
t  -dimensional  facet  in  its  boundary,  and  thus  must  contain  some  column  vec¬ 
tor  from  [  1 1  — M  ]  which  is  not  in  [/.?,  — M.y] .  Say  it  contains  J.n .  As  the 
r -dimensional  complementary  cones  covering  the  r -dimensional  facet  must 
be  generated  from  column  vectors  of  \I.7)  — M,?]  -  due  to  non  degeneracy  of 
M  -  then  the  interior  of  the  cone  pos[  q  j  I.n  ]  is  contained  in  the  interior  of 
the  (r-(-l)  -dimensional  facet.  Hence  the  reflecting  face,  pos[  C(a)  ~ j-  |  J.n  ] , 
which  is  the  common  face  between  the  cones  pos[  I  I  I-n  1  and 

pos[  C{a)  ~ j  |  — M.r  |  I.n  ],  contains  points  in  the  interior  of  the  (r  +  1)- 
dimensional  facet,  contradicting  the  regularity  of  that  facet.  This  completes 
the  induction.  Thus  all  the  r  -dimensional  facets  are  regular,  for  1  <  r  <  n. 

□ 

We  can  now  start  building  up  our  pseudomanifolds. 


Definition  4.16  For  any  complementary  cone,  C ,  in  K(M)  define  the 
pseudomanifold  P  =  P(C )  to  be  C  and  all  complementary  cones  C*  for 
which  there  exists  a  finite  sequence  of  complementary  cones  C  —  C\}Ci, . . ., 
Cm  —  C* ,  where,  for  1  <  *  <  m,  Ci  and  C,^ i  are  adjacent  cones  whose 
common  face  is  proper. 

Let  M  €  INSfc  fl  3fn><n  be  a  given  nondegenerate  matrix.  Thus  by 
Theorem  3.17  we  know  that  int  K (M)  is  connected.  Fix  some  C  C  K(M)  . 
Let  P  =  P(C) .  Let  q  and  q  be  two  points  in  int  K(M)  such  that  if  a 
complementary  cone  contains  one  of  these  points,  then  it  contains  that  point 
in  its  interior,  i.e.,  q,  q  £  K(M)  \  K(M)  .  We  can  now  use  Lemma  3.12 
to  get  a  path  L  from  q  to  q  satisfying  the  conditions  of  Lemma  3.12. 
Suppose  s  members  of  P  contain  q.  Now  move  along  L  from  q  to  q. 
When  L  crosses  a  face  of  a  complementary  cone,  that  face  must  be  proper 
as  LC  int  K(M)  and  K(M)  is  regular.  Thus  L  leaves  one  complementary 
cone  and  enter  another  one.  If  the  first  cone  was  a  member  of  P,  then 
the  second  cone  will  also  be.  Hence,  for  points  in  L  \  K(M)  ,  the  number 
of  members  of  P  that  contain  any  given  point  is  independent  of  the  point 
selected.  Thus  q  is  contained  in  s  members  of  P,  as  was  q.  Thus  every 
point  in  int  K(M)  \  K(M)  is  contained  in  s  members  of  P. 

Before  continuing  on,  let  us  digress  momentarily  to  point  out  a  simple 
fact  about  P.  Suppose  that  C*  is  a  cone  in  P.  By  definition  we  have 
the  sequence  of  cones,  C  —  Ci,  C2,  •  •  • ,  Cm  —  C * ,  adjacent  on  proper 
faces.  Suppose  Ci  —  pos  C(a)  and  C2  =  pos  C(/3).  By  the  definition  of 
a  proper  face,  we  have  (det  Maa)(det  Mpp)  >  0.  If  we  have  C*  =  pos  C^) , 
then  it  is  easily  seen  that  continuing  the  reasoning  in  the  last  sentence 
along  the  sequence  of  cones  C\,  C3, . . .  ,Cm  will  lead  us  to  conclude 
(det  Mac)(det  M^)  >  0.  Hence,  the  sign  of  the  determinant  of  the  prin- 


cipal  submatrix  of  M  associated  with  every  member  of  P  is  invariant  over 
P.  This  fact  will  be  useful  momentarily.  For  notational  purposes,  we  say 
a  full  complementary  cone  pos.C'(a)  is  positive  ( negative )  if  det  Maa  >  0 
(detMaa  <  0). 

Returning  to  the  main  discussion,  we  have  shown  every  point  in 
K(M)  \  K(M)  =  int  K (M)  \  K(M)  is  contained  in  s  members  of  P,  where, 
clearly,  s  >  1 .  Since  these  complementary  cones  are  closed,  and  by  non¬ 
degeneracy  K(M)  =  int  K(M)  \  K(M)  ,  it  follows  that  the  geometric  union 
of  the  members  of  P  is  K(M) . 

Now,  if  pos  C(a).i  is  a  face  of  exactly  one  member  of  P ,  i.e,  is  in  dP , 
then  it  must  be  in  dK (M)  .  For  if  it  contains  a  point  in  int  K(M)  ,  then  it 
must  be  a  proper  face,  which  would  imply  either  both  or  neither  of  the  cones 
containing  it  are  in  P.  Hence  the  geometric  union  of  the  members  of  dP 
must  be  contained  in  d  K(M)  . 

Let  q  €  dK{M)  be  such  that  if  any  complementary  cone  contains  it, 
the  cone  contains  it  in  the  interior  of  one  of  its  faces,  and  that  face  must  be 
contained  in  d  K(M)  .  (This,  as  usual,  allows  q  to  be  “almost  all”  the  points 
in  dS .)  Hence,  for  e  >  0  small  enough,  any  point  in 

B{q,  f)  D  int  K(M)  C  K(Af)  \  K{M)  (4.3) 

will  be  in  the  same  complementary  cones  as  q.  As  the  points  in  (4.3)  are 
contained  in  s  members  of  P ,  it  follows  that  q  is  contained  in  s  members 
of  P.  Suppose  that  pos  (7(a)  €  P  contains  q  in  its  face  pos  C(a).j.  Let 
/3  =  a  A  { i } ,  thus  pos  C(/3)  is  the  one  other  complementary  cone  with 
pos  C(a).i  as  a  face.  The  face  cannot  be  proper  as  it  is  in  d  K (M)  .  Thus  the 
face  is  reflecting,  and  so  (det  Maa)(det  Mpp)  <  0.  Since  pos<7(a)  €  P,  the 
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digression  above  shows  that  pos  C((3)  £  P .  Hence,  any  face  of  a  complemen¬ 
tary  cone  in  dK(M)  is  a  face  of  at  most  one  cone  in  P .  Thus  q  is  in  s 
members  of  dP . 

In  the  previous  paragraph  we  showed  that  any  cone  pos  C[a).p ,  where 
\/3\  =  n  —  1 ,  is  contained  in  at  most  one  member  of  P .  We  now  assume  that 
this  proper  holds  for  lower  dimensional  cones. 

Assumption  4.17  Let  M  £  INSn3inXn  be  nondegenerate,  and  let  C 
be  any  complementary  cone  in  K(M).  Let  F  be  any  r -dimensional  facet  of 
K(M),  1  <  r  <  n,  and  pos  C(a).p  C  3K(M)  ,  be  any  cone  in  dF  where 
a,0£  {n) ,  pos  C(a)  £  P(C) ,  and  \fl\  =  r  —  1 .  Then  there  exists  at  most 
one  7  £  (n)  with  I7)  =  r  such  that,  posC^a).^  C  F  and  [3  C  7. 

(This  assumption  is  essentially  the  “consistency”  assumption  used  in  the 
previously  cited  work  of  Doverspike  and  Lemke.) 

We  can  now  state  the  main  theorem  of  this  section. 

Theorem  4.18  Let  M  £  INSjtn9inX”  be  nondegenerate.  If  Assump¬ 
tion  4.17  holds,  then  the  complementary  cones  of  K(M)  can  be  partitioned 
into  k  disjoint  collections  where  each  collection  is  an  orientable  (n  —  1)- 
dimensional  pseudomanifold  by  the  representation  described  in  Example  4.14. 
Furthermore  if  P  is  one  of  these  pseudomanifolds,  then  the  geometric  union 
of  the  cones  in  P  equals  K(Af).  Also,  if  pos  C[a),  pos  C(/3)  £  P,  then 
int  pos C(a) H  int  pos  C(/3)  =  0.  (In  this  way  each  pseudomanifold  partitions 
K(A/).)  In  addition,  the  (n  —  l)-faces  making  up  the  boundary  of  P ,  call  it 
dP ,  also  have  disjoint  interiors  and  their  union  is  geometrically  dK{M)  .  (It 
is  known  that  dP  will  be  an  orientable  (n —  1) -dimensional  pseudomanifold 
without  boundary.) 
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Proof.  Most  of  the  work  has  been  already  done.  We  will  use  reverse 
induction.  Suppose  by  induction  we  have  a  sequence  of  facets  of  K(M),  say 
Fr,Fr+i, . . .  ,Fn_i ,  and  dimF*  =  t  for  r  <  *  <  n.  In  addition,  suppose 
for  each  F,-  there  is  a  collection,  F, ,  of  i -dimensional  facets  of  members 
of  P  =  P[C) ,  each  facet  being  in  Fj ,  such  that  for  each  q  6  F*  which  is 
not  contained  in  any  (i  —  1) -dimensional  facet  of  any  complementary  cone, 
(which  is  “almost  all”  of  F, ,)  there  are  exactly  s  members  of  Pi  containing 
?• 

We  have  already  shown  we  may  start  the  induction  by  taking  Fn_x  as 
any  (n  —  1) -dimensional  facet  of  K(M),  and  selecting  as  Fn_i  all  those 
faces  of  members  of  P  that  lie  in  Fn_i .  Now  suppose  we  are  at  the  general 
induction  step.  Take  as  Fr_i  any  boundary  facet  of  Fr .  For  any  q  E  Fr_i 
that  is  not  contained  in  any  (r  —  2)  -dimensional  facet  of  any  complementary 
cone  there  is  an  e  >  0  small  enough  so  that  each  member  of  Pr  either 
contains  or  is  disjoint  from  B(q,  e)C\Fr.  By  induction,  for  small  e  >  0,  we 
have  each  point  in  B{q,  e)  n  int  Fr  must  be  in  exactly  s  members  of  Pr  ,  thus 
q  is  in  s  members  of  Pr .  By  Assumption  4.17,  each  (r  —  1) -dimensional 
facet,  of  a  complementary  cone,  contained  in  Pr—i  must  be  the  face  of  no 
more  than  one  member  of  PT .  Hence,  if  we  define  Fr— l  as  all  the  faces  of 
members  of  Pr  contained  in  Fr_i ,  the  point  q  is  contained  in  exactly  s 
members  of  Pr~i .  Noticing  that  the  members  of  Pr—i  must  be  facets  of 
members  of  P  completes  the  induction. 

The  only  “catch”  in  the  induction  is  where  we  assume  that  Fr  has 
a  boundary  facet.  Suppose  it  doesn’t,  and  hence  Fr  is  an  r -dimensional 
subspace.  By  the  nondegeneracy  of  M,  we  may  assume  I.i  and  — 
are  not  in  Fr,  for  r  <  i  <  n.  Also,  notice  Fr  must  be  covered  by  the 
r-dimensional  facets  of  complementary  cones  that  are  in  Fr .  If  there  is 
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some  column  vector  of  [  I.7  |  — M.f  ]  that  is  not  in  Fr ,  say  I.x ,  then  every 
r  -dimensional  facet,  of  a  complementary  cone,  that  is  in  Fr  must  contain 
— Mi  and  so  by  nondegeneracy  must  not  contain  M.\ .  This  contradicts 
the  fact  that  F ,  is  covered.  Thus  Fr  contains  all  the  column  vectors  of 
[  I.r  |  — Mf  ] .  We  can  view  Fr  as  5Rr ,  and  [  I.?  |  — My  ]  as  defining  a 
LCP,  where  the  matrix  of  the  LCP  is  Mr  =  — Myy.  We  know  Mr  €  Q. 
By  Lemma  4.15  K(Mr)  is  regular  so,  by  Theorem  3.18,  M  6  INS.  Hence 
Mr  6  P  •  Thus  s  =  l. 

Now  suppose  that  we  can  always  continue  the  induction  down  to  Fi , 
no  matter  what  choices  we  make  along  the  way.  By  nondegenracy,  F\ 
can  contain  at  most  two  column  vectors  from  [  /  |  — M  ] .  If  it  contains 
only  one  such  vector  then  we  have  a  =  1  as  before.  If  it  contains  two 
such  vectors  then  they  are  /.,  and  — M.,  for  some  *  6  n.  (In  this  case 
a  =  2.)  Hence  — M,  =  X/.,  for  some  X  >  0.  As  the- boundary  of  K(M) 
contained  no  lineality  (no  linear  subspace),  there  must  be  a  minimum  of  n 
1-dimensional  facets.  Hence,  for  some  i  G  n,  each  such  facet  must  contain 
— M.%  and  /.* .  Each  facet  must  be  associated  with  a  different  i.  Thus  M 
is  a  diagonal  matrix  with  negative  diagonal  entries.  It  is  easily  seen  that 
for  each  complementary  cone  C ,  we  have  P(C )  =  {C}.  Hence  s  =  1,  a 
contradiction. 

In  all  cases  we  have  have  3  =  1.  Thus  any  two  cones  in  P  must  have 
disjoint  interiors,  otherwise  the  intersection  would  be  n -dimensional  which 
would  mean  some  of  the  points  in  the  intersection  are  in  K(M)  \  K(M) 
and  would  have  to  be  in  only  one  member  of  P.  The  same  can  be  said  for 
the  members  of  dP .  Since  any  point  in  K(M)  \  K(M)  is  contained  in 
k  complementary  cones,  as  M  €  INS* ,  and  each  complementary  cone  C 
is  contained  in  some  pseudomanifold  P,  for  example  P{C),  then  the  corn- 
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plementary  cones  can  be  partitioned  into  k ,  clearly  disjoint,  collections  with 
each  collection  forming  a  pseudomanifold.  Each  pseudomanifold,  P[C ) ,  can 
be  oriented,  as  in  Example  4.14,  by  giving  the  ordering  (C7(a).i, . . . ,  C(a).n) 
of  pos  C(a)  the  sign  ( — l)lal .  This  induces  orientations  for  the  boundary 
pseudomanifolds.  See,  for  example,  Freund  (1980). 

□ 


As  a  final  remark,  it  should  be  mentioned  that  the  boundary  pseudo¬ 
manifolds  need  not  be  distinct.  For  example,  Figure  4.5  shows  K(M )  for 


In  this  case  M  £  INS4 ,  and  the  four  pseudomanifolds  are  the  four  com¬ 
plementary  cones.  Each  has  a  different  boundary  pseudomanifold.  Figure  4.6 
shows  K(M)  for 


M  = 


1 

1 


Here  M  €  INS2 ,  and  the  two  pseudomanifolds  are  {  pos  C(0),  pos  C({  2  }) } 
and  {  pos  C({  1  }),  pos  C(%) } .  These  both  have  the  pseudomanifold  {  pos  1. 2, 
pos  — M. 2  }  as  boundary. 


4.4  A  Simple  Class  of  INS-matrices 


The  relation  between  INS-matrices  and  other  matrix  classes  will  be  dis¬ 
cussed  in  Chapter  5,  however,  it  seems  appropriate  at  this  point  to  introduce 
a  simple  subclass  of  INS.  So  saying,  we  have 
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Definition  4.19  We  say  that  a  matrix  A  is  in  the  class  GNI  (Gener¬ 
alized  Negative  Identity)  if  and  only  if 

A  6  {ME  3?lXn  :  |supp  Af.*|  <  1,  for  all  t  6  n  } , 
ngZ+ 

i.e.,  each  column  of  the  matrix  contains  at  most  one  non-zero  entry,  and,  if 
it  exists,  this  non-zero  entry  is  negative. 

Suppose  M  E  GNI  and  a  €  (n).  If  pos  C(a)  is  a  full  cone  then 
pos  C(a)  =  .  Otherwise,  pos  G(a)  C  <99?^  .  Thus  no  face  of  any 

complementary  cone  intersects  the  interior  of  K(M)  =  9?I}_ .  Accordingly, 
int  K(M)  =  int3?^_  is  itself  one  of  the  connected  components  in  £.  So,  by 
Theorem  3.8, 

GNI  C  INS . 

GNI-matrices  satisfy  the  conclusion  of  the  theorems  of  Section  4.2,  and 
the  proof  sheds  light  on  the  combinatorial  aspect  of  the  subject.  In  fact, 
the  theorem  essentially  follows  from  the  next  lemma  which  is  an  interesting 
combinatorial  result  by  itself. 

Lemma  4.20  Suppose  we  are  given  n  boxes,  labelled  1, 2, . . . ,  n ,  and 
2 n  balls,  labelled  1, 2, . . . ,  n,  i,  2, . . .  ,  n  .  Suppose  also  that,  for  all  i  E  n, 
ball  i  is  in  box  t ,  whereas  ball  t  may  be  in  any  one,  or  none,  of  the  boxes. 
Say  that  (/j,  I2, . . . ,  ln)  is  a  list  if  for  all  *  £  n,  U  equals  j  or  j  for  some 
j  E  n,  and  ball  /,  is  contained  in  box  >.  Say  that  a  list  is  proper  if  for  all 
*  E  n  there  exists  a  j  En  such  that  lj  E  { *,  i  }  .  Then  the  number  of  proper 
lists  is  a  power  of  two. 

Proof.  This  will  be  by  induction  on  n .  If  n  =  1 ,  the  number  of  proper 
lists  is  2  or  1  depending  on  whether  ball  I  is,  respectively,  in  or  not  in  box  1. 
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The  lemma  is  true  in  this  case. 


Assume  the  lemma  is  true  for  1, . . . ,  n  —  1 .  We  will  show  it  true  for  n . 

Suppose  some  bail  is  in  no  box.  We  may  assume  it  is  ball  i .  Then  any  proper 

list  (lh  1%, ...  i  In)  has  ln  =  n.  Notice  that  (/*,  I2, . . . ,  ln— 1)  is  a  proper 

list,  and  any  such  proper  list  can  be  extended  to  a  complete  proper  list  by 

adjoining  ln  =  n.  Also  the  distribution  of  the  balls  1, . . . ,  n —  1, 1, . . . ,  n  —  1 
in  the  boxes  1  ,...,n  —  1  satisfies  the  conditions  of  the  lemma.  Thus  by 
induction  the  number  of  proper  lists  (l\, . . .  ,ln— 1)  is  a  power  of  two,  and 
this  equals  the  number  of  complete  proper  lists  (li, . . . ,  ln) • 

Assume  all  the  balls  I, . . . ,  n  are  each  in  some  box.  Suppose  ball  t  is 
in  box  t.  Again,  we  may  assume  t  =  n.  Then,  as  above,  any  proper  list 
(lit  •  •  • , In— 1)  can  be  made  into  a  complete  proper  list  by  adjoining  either 
ln  —  n  or  by  adding  in  ln  —  h.  (Notice  that  either  is  possible.)  Also,  any 
proper  list  will  have  ln  =  n  or  ln  —  n,  hence  (£1, ...» /n — 1)  is 

a  proper  list.  As  above,  we  may  use  induction  to  show  that  the  number  of 
proper  lists  (£j_, . . . ,  in_i)  is  a  power  of  two.  Thus  we  have  twice  that  number 
of  complete  proper  lists.  This  is  still  a  power  of  two. 

Now  suppose,  for  all  ign,  that  we  have  j  in  some  box,  but  not  box  i . 
Let  i'i,  t'2,  »3, ...  be  a  sequence  defined  by  letting  ii  =  1  and  saying  that  ball 
ij  is  in  box  ij+i  for  all  j  €  Z+ .  Then  the  sequence  must  clearly  repeat  a 
number  at  some  point,  say  ij  =  ik ,  such  that  j  <  k,  and  ij,  ij+t, ...» 1 
are  all  distinct.  We  may  assume  that  j  =  1 ,  that  3  <  k  <  n  -(-  1 ,  and  the 
sequence  at  the  end  of  the  last  sentence  is  1,2,..., A; —  1.  Let  (li,...,ln) 
be  a  proper  list.  If  l\  =  1 ,  then  as  ball  k  —  1  is  in  box  1,  we  need  to  have 
lk—i  =  k  —  1 .  As  ball  k  —  2  is  in  box  k  —  1 ,  we  need  to  have  Ik— 2  =  k  —  2 . 
Continuing  in  this  fashion  we  find  that  (llt...,lk— 1)  is  (1, . ...  ,fc —  1).  If 
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l\  7^  1 ,  then  we  need  I2  —  1 .  As  Z2  7^  2 ,  we  need  that  Z3  =  2 .  Continuing 
in  this  fashion  we  have  that  (h, lk~  1)  is  (fc  —  1, 1, 2, . . . , k  —  2) .  Thus 
we  must  have  that  (Zj,  equals  (1, . . . ,  k  —  1)  or  (A:  —  1,1 ,k  —  2) . 

Notice  that  either  one  of  these  two  will  do,  since  having  either  one  of  these 
begin  the  complete  proper  list  will  force  the  rest  of  the  proper  list,  (Ik, . . n) , 
to  be  selected  from  the  set  {  k,  .  ,n}.  Hence  given  an  “ending” 

to  the  proper  list  that  works  with  either  of  the  previous  two  “beginnings,” 
the  “ending”  will  work  with  both  of  the  “beginnings.”  Furthermore,  we 
see  the  “ending”  is  just  a  proper  list  for  the  boxes  fc,  ...,n  using  balls 
fc, . . . , n,  fc, . . . , n ,  and  any  such  proper  list  will  do.  By  induction,  the  number 
of  such  proper  lists  for  the  “ending”  is  a  power  of  two.  Since  there  are  two 
possible  “beginnings,”  the  number  of  complete  proper  lists  is  also  a  power  of 
two.  This  completes  the  induction,  and  the  lemma  follows. 

□ 


The  above  lemma  translates  almost  immediately  into  the 
Theorem  4.21 

OO 

GNI  C  (J  INS2, . 
v— 0 

Proof.  Let  M  £  GNln3?nXn.  With  reference  to  Lemma  4.20,  ball  * 
corresponding  to  and  ball  t  corresponding  to  — M.,  for  *  £  n.  We  say  a 
ball  is  in  box  i  if  and  only  if  the  ith  component  of  the  corresponding  vector  is 
nonzero.  Thus  there  is  a  bijective  correspondence  between  full  complementary 
cones  and  proper  lists,  where  the  elements  of  a  proper  list  correspond  to  the 
columns  of  a  nondegenerate  complementary  matrix.  Each  of  these  full  cones 
is  equal  to  =  K (M) ,  and  as  Lemma  4.20  now  tells  us  the  number  of 
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such  cones  is  a  power  of  two,  we  have  M  €  INS2P  for  some  nonnegative 
integer  p. 

□ 

It  seems  bothersome  to  require  that  the  nonzero  entries  in  a  GNI-matrix 
be  negative.  It  would  be  preferable  to  work  with  the  following  class  of 
matrices. 

Definition  4.22  We  say  that  a  matrix  A  is  in  the  class  GI  if  and  only 
if 

A€  {  M  €  3?nXn  :  |supp  \  <  1,  for  all  *  6  n  } , 

*»€  Z+ 

i.e.,  each  column  of  the  matrix  contains  at  most  one  non-zero  entry. 

Unfortunately,  as  seen  at  the  end  of  Section  2.1,  the  matrix 

0  -1 
.10. 

is  not  INS,  but  is  in  GI.  See  Figure  2.3.  Thus  GI  INS .  However,  it 
is  “close”  enough  to  warrant  investigation,  and  so  we  look  at  the  following 
combinatorial  lemma  which  is  an  extension  of  Lemma  4.20.  (The  proofs  are 
almost  identical  so  the  proof  of  Lemma  4.23  will  be  given  in  less  detail  than 
necessary,  but  familiarity  with  the  reasoning  in  the  proof  of  Lemma  4.20  will 
be  assumed.) 

Lemma  4.23  Suppose  we  are  given  2n  boxes  which  are  labelled  1,2,..., 
n,  1/,  2 ',...,  n' ,  and  2n  balls,  labelled  1, 2, . . . ,  n,  1, 2, . . . ,  n .  Suppose  also 
that,  for  all  i  €  n,  ball  i  is  in  box  t,  but  ball  i  can  be  in  any  one  box,  or 
no  box  at  all.  Say  that  (61, 62  >  •  •  •  ,  6n)  is  a  box  list  if,  for  each  i  €n,  bi  is 
either  i  or  i' .  Furthermore,  given  a  box  list,  we  say  that  (li,h,  •  -  ■ ,  U)  is  a 
list  for  the  box  list  if,  for  each  t  €  n,  li  equals  j  or  j  for  some  j  6  n,  and 
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li  is  contained  in  box  bi .  Say  that  a  list  is  proper  if,  for  all  i  6  n,  there  is  a 
j  €  n  such  that  lj  €  •{»,*}.  Then,  there  exists  a  nonnegative  integer  p  such 
that  the  number  of  proper  lists,  associated  with  any  box  list  is  either  zero  or 
2P . 

Proof.  We  will  use  induction  on  n .  For  n  =  1  either  1  is  in  no  box  (box  1), . 
in  which  case  box  1  has  one  (two)  proper  list(s)  and  box  V  has  none,  else  1 
is  in  box  1/ ,  in  which  case  both  boxes  have  one  proper  list.  The  lemma  holds 
here. 

Now  assume  the  lemma  holds  for  1, . . . ,  n  —  1 .  We  will  show  it  true  for 
n .  Suppose  some  ball,  say  n ,  is  in  no  box.  Then  any  box  list  with  at  least 
one  proper  list  must  have  bn  —  n .  Also,  any  proper  list  must  have  ln  =  n. 
Similar  to  before,  we  find  that  the  number  of  proper  lists  that  are  associated 
with  box  lists  of  the  form  (6j,, . . . ,  6n_i,  n)  equals  the  number  of  proper  lists 
(lit  •  •  • ,  In— l)  for  the  associated  box  lists  (&i, . . . ,  6n— l)  when  we  consider 
the  embedded  smaller  problem  for  n  —  1 .  The  lemma  then  holds  here  by 
induction. 

Assume  all  the  balls  1,2 ,  ...,n  are  in  some  box.  If  for  some  »,  say 
t  =  n ,  we  have  ball  n  in  box  n ,  then  we  will  have  a  situation  similar  to  the 
above.  Any  box  list  will  have  no  proper  lists  if  bn  =  n' ,  and  the  number  of 
proper  lists  of  box  lists  in  the  form  (&i, . . . ,  6n__i,  n)  equals  twice  the  number 
of  proper  lists  [1%, . . . ,  ln— i)  for  the  box  lists  (6 1, . . . ,  bn— l)  when  we  consider 
the  smaller  embedded  problem.  (We  would  just  add  on  ln  =  n  and  ln  —  n 
to  get  the  two  complete  proper  lists  from  the  box  lists  of  the  smaller  problem.) 
The  lemma  then  holds  here  by  induction. 

If  we  had  n  in  box  n'  in  the  previous  paragraph,  then  each  box  list 
(&i, . . . ,  6n)  would  have  the  same  number  of  proper  lists  as  (&i, . . . ,  6n_i)  in 
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the  smaller  embedded  problem.  The  reasoning  is  the  same  as  before,  only 
now  we  complete  the  smaller  proper  lists  by  adding  on  ln  =  n  if  bn  —  n, 
and  ln  =  h  if  bn  =  n' .  The  lemma  will  still  hold  by  induction. 

So  now,  finally,  assume  that  all  the  balls  1, ...,n  are  in  some  box  but, 
for  each  i  En,  ball  t  is  neither  in  box  i  nor  in  box  .  Thus,  as  in  the 
proof  of  Lemma  4.20,  we  may  assume  for  some  k ,  where  3  <  k  <  n  -{-  1 , 
that,  for  all  t  6  k  —  2 ,  the  ball  ?  is  in  either  box  t  -j-  1  or  box  (t  +  1)' , 
and  the  ball  k  —  1  is  in  either  box  1  or  box  1' .  Now  suppose  the  box  list 
(bj, . . . ,  6n)  has  a  proper  list  (l:  . . . ,  L) .  Suppose  1^1.  Thus  bi  =  1 ,  so 
no  other  ball  in  the  proper  list  be  from  box  1  or  box  1' .  Hence  we  need 
lk—i  =  k —  1 .  Continuing  on  in  this  fashion,  as  in  the  proof  of  Lemma  4.20, 

we  get  (£i, ...» _ i)  is  (1,  ..  .,&  —  1).  If  li  ^  1 ,  then  we  need  i2  .==  1  • 

Thus  &2  equals  whichever  of  the  two  boxes,  2  or  2' ,  contains  1 .  In  either 
case,  we  cannot  select  any  other  ball  from  either  of  the  two  boxes  for  the 
proper  list,  hence  we  need  f3  =  2 ,  and  continuing  on  we  have  (£i , . . . ,  lk—i) 
is  (fc  —  1, 1, 2, . . . ,  k  —  2) .  Each  of  these  cases  determines  the  “beginning” 
of  the  box  list,  i.e.,  (6i,  b2, . . . ,  bfc— i).  If  the  list  in  each  case  is  the  same, 
then  any  box  list  which  has  at  least  one  associated  proper  list  must  have 
this  “beginning.”  The  ending,  as  in  Lemma  4.20,  can  be  any  proper  list 
(Ik,  •  ■ . ,  ln)  from  the  smaller  embedded  problem.  By  induction  this  is  a  fixed 
power  of  two,  say  2P ,  and  so  the  number  of  complete  proper  lists,  for  any  box 
list  having  proper  lists,  is  2P+1 .  If  the  two  “beginnings”  are  different,  then 
any  box  list  having  proper  lists  must  have  one  of  these  two  “beginnings,”  and 
the  number  of  proper  lists  it  will  have  will  be  2P . 

In  all  cases,  all  box  lists  with  proper  lists  have  the  same  number  of  proper 
lists,  and  that  number  is  a  power  of  two.  The  induction  is  now  completed. 

□ 
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Now,  with  this  lemma,  we  can  finish  this  section  with  the  following 

Theorem  4.24  If  M  6  Gln3ln><n,  there  exists  a  nonegative  integer 
p,  such  that  for  all  q  £  K(M) 

jsupp  q\  =  n  =>  |sol(g,  M)\  =  2P  . 


Proof.  With  reference  to  Lemma  4.23,  let  ball  t  correspond  to  I.i ,  and  ball  t 
correspond  to  — M.i .  We  say  a  ball  is  in  box  t  if  and  only  if  the  corresponding 
vector  has  it’s  ith  component  positive.  We  say  a  ball  is  in  box  i'  if  and 
only  if  the  corresponding  vector  has  it’s  ith  component  negative.  Each 
full  complementary  cone  must  be,  geometrically,  an  orthant  in  3?n .  Each 
degenerate  complementary  cone  must  be,  geometrically,  contained  in  the 
union  of  the  boundaries  of  the  orthants.  There  is  a  bijective  correspondence 
between  orthants  and  box  lists.  Thus  the  interior  of  each  orthant  that  is 
contained  in  K(M)  must  be  an  element  in  £ .  Since  there  is  a  bijective 
correspondence  between  full  cones  covering  an  orthant  and  proper  lists  of  the 
orthant’s  associated  box  list,  Lemma  4.23  implies  that  each  orthant  contained 
in  K(M),  i.e.,  with  some  associated  proper  list,  must  be  covered  by  the  same 
number  of  full  cones  as  the  other  orthants  in  K(M),  and  that  number  must 
be  a  power  of  two.  Thus,  by  Theorem  3.8,  |sol(g,  M)\  is  this  power  of  two 
for  any  q  belonging  to  K(M)  and  the  interior  of  an  orthant. 


□ 


CHAPTER  5. 

MATRIX  CLASSES  AND  LCP  THEORY 


5.1  Matrix  Classes 

Much  of  the  literature  concerning  the  LCP  deals  with  the  study  of  matrix 
classes.  Some  classes  arc  defined  using  the  LCP  itself  and  so  we  seek  more 
constructive  characterizations.  Other  classes  are  defined  using  more  simple 
and  testable  criteria  and  results  are  found  concerning  the  nature  of  the  LCP 
(q,  M )  when  M  is  in  one  of  these  classes.  The  relationships  among  the  classes 
has  also  been  a  rich  subject  of  study,  and  much  work  has  been  devoted  to 
trying  to  understand  which  basic  properties  of  importance  to  the  LCP  are 
common,  or  different,  among  the  matrix  classes.  In  Figure  5.1  we  have  listed 
the  seven  matrix  classes  defined  in  this  work  along  with  some  of  the  more 
well-studied  matrix  classes  in  the  field.  This  figure  should  be  referred  to 
throughout  this  section.  (The  arrows  indicate  inclusion  relationships  among 
the  classes,  with  the  larger  classes  tending  to  be  at  the  top  of  the  page.) 
The  purpose  of  this  section  is  to  define  the  classes  in  this  “family  tree,”  and 
to  discuss  just  where  U  and  INS  fit  into  it.  There  is  no  attempt  to  give  a 
detailed  review  of  these  classes,  however  references  are  given  showing  where 
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more  information  can  be  obtained.  Some  basic  references  of  general  value 
are  Lemke  (1970),  Karamardian  (1972),  Kostreva  (1976),  Mohan  (1978),  and 
Cottle  (1983).  The  classes  are  presented  in  alphabetical  order  by  the  symbols 
used  in  Figure  5.1.  At  times  it  will  be  necessary  to  refer  to  the  definition  of 
a  matrix  class  not  yet  given. 

(A)  A  matrix  M  6  D?nXn  is  said  to  be  adequate,  M  6  A,  if  and  only 
if  M  6  Po  and  for  a  €  (n)  we  have  detMaa  =  0  implies  the  column 
vectors  M.a  are  linearly  dependent  and  the  row  vectors  Mas  are  linearly 
dependent.  See  Ingleton  (1966),  Cottle  (1968)  and  Eaves  (1971). 


(BG)  A  matrix  M  €  9inXn  is  said  to  be  a  bimatrix  game  matrix, 
M  E  BG,  if  and  only  if  for  some  m  E  (n  —  1 )  there  are  matrices 
A  e  3r*X(*-m)  and  Q  g  gj(n  m}Xm  where  A>  B  >  Q  and 


See  Lemke  and  Howson  (1964),  Cottle  and  Dantzig  (1968)  and  Eaves  (i971). 


(CP)  A  matrix  M  E  SR71*71  is  said  to  be  copositive,  M  E  CP,  if  and 
only  if  for  all  x  E  9?n ,  z  >  0  implies  xT Mx  >  0.  This  matrix  class  has  also 
been  denoted  as  Co .  Copositive  matrices  are  important  in  combinatorics  and 
other  fields  aside  from  complementarity.  There  is  a  large  literature  about  this 
class,  for  example,  see  Gaddum  (1958),  Cottle  and  Dantzig  (1968),  Cottle, 
Habetler,  and  Lemke  (1970b),  Pereira  (1972),  Hoffman  and  Pereira  (1973), 
and  Evers  (1978). 


(C+)  A  matrix  M  E  3?nXn  is  said  to  be  copositive-plus,  M  E  C+,  if 
and  only  if  M  is  copositive  and  for  all  x  E  ,  x  >  0  and  xT M x  =  0  imply 
( M  +  MT)x  —  0 .  Like  the  copositive  matrices,  there  is  a  large  literature 
concerned  with  these  matrices.  See  the  papers  given  as  references  for  the 
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copositive  matrices. 


(Eo)  A  matrix  M  £  3inXn  is  said  to  be  semi-monotone,  M  £  Eo  , 
if  and  only  if  for  all  x  £  3in ,  C  5^  x  >  0  implies  there  is  some  k  £  n 
such  that  Ik  >  0  and  ( Mx)k  >  0.  (This  class  has  also  been  denoted 
as  Li .)  If  M  is  symmetric,  then  M  is  semi-monotone  if  and  only  if 
M  is  copositive.  We  have  used  these  matrices  previously,  with  their  other 
characterization  of  being  the  class  of  matrices  M  for  which  |sol(g,  M)\  —  1 
for  all  0  <  q  £  5Rn  .  Like  the  copositive  matrices,  the  semi-monotone  matrices 
have  been  extensively  studied.  See,  for  example,  Lemke  (1970),  Eaves  (1971), 
Pereira  (1972),  Karamardian  (1972),  and  Garcia  (1973). 

(Eg)  A  matrix  M  £  8inXn  is  said  to  be  fully  semi-monotone  if  and 
only  if  all  principal  transforms  of  M  are  semi-monotone.  This  matrix  class 
was  introduced  in  this  work,  and  was  she  -  n  to  contain  the  matrix  classes 
U  and  Po .  (It  is  clearly  contained  in  Eo .)  As  seen,  it  can  be  characterized 
as  the  class  of  matrices  such  that  for  all  q  £  9?" ,  if  (w,  z)  £  sol (<7,  M)  and 

t 

w  -f-  z  >  0  then  { (to,  z) }  =  sol(g,  M) .  To  see  why  Eg  has  been  placed 
where  it  is  in  Figure  5.1,  consider  the  matrices 


(5.1)  (5.2)  (5.3)  (5.4)  (5.5) 

None  of  these  matrices  are  in  Eg.  However,  (5.1)  is  in  SCP,  E,  and 
(N  D  Q)—1 ,  (5.2)  is  in  Z ,  (5.3)  is  in  N  D  Q ,  (5.4)  is  in  N  \  Q ,  and  (5.5)  is  in 
GNI.  Consider  now 

0  0  12 
0  0  2  1 

12  0  0 

2  10  0 
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M  is  in  BG,  but  is  not  in  Eg  .  This  can  be  seen  as  we  have 


and  the  inverse  of  a  matrix,  if  it  exists,  is  always  a  principal  transform. 
However,  with  x  =  (1,0,  l,0)r,  we  note  that  there  is  no  index  A:  €  4  for 
which  it  >  0  and  (M~1x)/e  >  0.  Hence  M  g  Eg  . 

(E)  A  matrix  M  6  3RnXn  is  said  to  be  strictly  semi-monotone  if  and 
only  if  for  all  x  6  R” ,  0  x  >  0  implies  there  is  some  k  €  n  such  that 
I*  >  0  and  (Mi)i  >  0.  (This  class  has  also  been  denoted  as  L* .)  If 
M  is  symmetric,  then  M  is  strictly  semi-monotone  if  and  only  if  M  is 
strictly  copositive.  Similar  to  the  semi-montone  matrices,  these  matrices  can 
be  characterized  as  being  the  class  of  matrices  M  for  which  |sol(g,  M)\  =  1 
for  all  0  <  q  €  9tn  .  See  the  papers  given  as  references  for  the  semi-monotone 
matrices.  This  matrix  class  is  also  the  class  of  completely  Q-matrices,  which 
are  defined  to  be  those  Q-matrices  all  of  whose  principal  submatrices  are  also 
Q-matrices.  This  equivalence  was  shown  by  Cottle  (1979). 

(GI)  A  matrix  M  €  5Rn><n  is  said  to  be  in  GI  if  and  only  if  for  all  tgn 
we  have  |suppM.j|  <  1.  This  class  was  brought  up  in  Chapter  4  due  to  its 
combinatorial  nature,  and  because  it  is  "almost”  in  the  class  INS.  For  such 
a  simple  class  of  matrices,  it  seems  surprising  that  it  is  contained  in  none 
of  the  other  matrix  classes  in  Figure  5.1.  Still,  Example  2.3  is  a  Gl-matrix 
that  is  not  in  Qo ,  and  the  1  X  1  matrix  (these  are  usually  referred  to  as 
“numbers”)  [—1]  is  not  in  So .  As  mentioned  in  Chapter  4,  GI  2  INS.  The 
class  GI  is  contained  in  no  other  matrix  class  in  Figure  5.1,  since  every  other 
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matrix  class  shown  there  is  a  subclass  of  Qo ,  So ,  or  INS. 

(GNI)  A  matrix  M  6  SfnXn  is  said  to  be  in  GNI  if  and  only  if  M  €  GI 
and  M  <  0.  It  was  shown  in  Chapter  4  that  this  class  is  in  INS.  In  fact, 
M  €  INS/cflGNI  implies  k  =  2P  for  some  nonnegative  integer  p.  Also 
INS2P  D  GNI  7*^  0  for  all  nonnegative  integers  p,  as  the  zero  matrix  is  in  INSi 
and  —  J€SRpXp  is  in  INS2>> . 

(INS)  A  matrix  M  6  SUnXn  is  said  to  be  in  INS,  for  Invariant  Number 
of  Solutions,  if  and  only  if  there  is  some  positive  integer  k  such  that  for  all 
q  €  int  K(M)  we  have  |sol(g,M)|  =  k.  We  have  studied  these  matrices 
a  great  deal.  Notice  now  where  they  fit  into  Figure  5.1.  We  know  from 
Theorem  3.4  that  INS  fl  Q  =  P .  Also,  it  is  shown  in  Garcia  (1973)  that 
M  €  L(d)  with  d  >  0  implies  |sol(d,  M)|  =  1 ,  and  hence  we  have 

U  L(d)nDMS=U. 

d>  0 

We  see  that  Eo  nINS  =  U,  since  for  M  €  Eo  we  have  |sol(<jr,  M)\  —  1  for 
all  q  >  0 .  More  will  be  said  about  INS  matrices  in  relation  to  some  of  the 
other  classes,  but,  before  moving  on,  notice  that  the  matrix  given  in  (5.5)  is 
in  INS  but  is  not  in  So ,  so  INS  So  . 

(K)  A  matrix  M  €  3ln  * n  is  said  to  be  in  K  if  and  only  if  M  6  P  n  Z . 
(These  matrices  have  also  been  referred  to  as  the  Minkowski  matrices  and 
denoted  as  the  class  M.)  These  matrices  have  a  great  deal  of  structure, 
both  geometric  and  algebraic.  It  is  interesting  to  note  K  =  Z  fl  Q ,  i.e.,  the 
complementary  cones  of  a  Z-matrix  cover  DJn  if  and  only  if  they  partition  3?n . 
(The  meaning  of  “partition”  allows  the  cones  to  intersect  on  their  boundaries.) 
The  classic  reference  for  these  matrices  is  Fiedler  and  Ptak  (1962).  See  also 
Cottle  and  Veinott  (1972). 
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(Kq)  A  matrix  M  G  9inXn  is  said  to  be  in  Ko  if  and  only  if 
M  G  Po  fl  Z .  Again,  the  classic  reference  here  is  Fiedler  and  Ptak  (1962). 
In  Mohan  (1980),  it  is  shown  that  the  boundary  of  a  Ko  matrix  is  the  union 
of  the  degenerate  faces.  Since,  for  MG  Ko  ,  there  are  no  reflecting  faces  in 
K(M) ,  as  Ko  C  P0  ,  it  follows  that  K (M)  is  regular.  In  Chandrasekaran 
(1970)  it  was  shown  that  Z  G  Qo  >  hence  K(M)  is  convex  for  a  K0  -matrix, 
and  so  int  K(M)  will  be  connected.  Thus  Ko  C  INS  by  Corollary  3.14. 
As  Kq  C  P0  C  Eo  ,  we  see  we  must  have  Ko  C  U.  In  Mohan  (1980), 
other  results  are  derived  about  Ko  which  can  be  viewed  as  consequences  of 
some  of  the  theorems  presented  here  concerning  U-matrices.  See  also  Mohan 
(1978)  for  more  on  Ko -matrices. 

(L)  A  matrix  M  G  3?nXn  is  said  to  be  in  L  if  and  only  if  M  G  E0  , 
and  for  all  (tu,  z )  G  sol(0,  M) ,  where  z  ^  0 ,  there  is  a  x  G  3ftn ;  0  x  >  0 , 
with  z  >  x  and  w  >  — MTx  >  0.  This  is  one  of  the  largest  classes  of 
matrices  that  Lemke’s  algorithm  using  e  =  (1, 1, . . . ,  1)  is  known  to  process. 
The  standard  reference  for  this  class,  which  is  also  the  reference  defining  the 
class,  is  Eaves  (1971). 

(L (d))  A  matrix  Af  G  SnXn  is  said  to  be  in  L[d )  if  and  only  if  for  all 
(w,  z)  G  sol(Xd,  M) ,  where  z  ^  0  and  X  >  0 ,  there  is  a  xG9Rn,  07^x>0, 
with  z  >  x  and  w  >  —MTx  >  0.  The  standard  reference  for  these  classes 
is  Garcia  (1973).  It  should  be  pointed  out  that  L  =  D£j>oL(^)* 

(L*(d))  A  matrix  M  G  ?Rn><n  is  said  to  be  in  L*(d)  if  and  only  if 
for  all  X  >  0  we  have  that  ( w,z )  G  sol(Xd,  M)  implies  (t u,z)  =  (Xd,0). 
For  d  >  0,  L*(d)  is  the  class  of  all  matrices  M  where  in  K (M)  the  only 
complementary  cones  containing  d  are  pos  C{a)  where  a  fl  supp  d  =  0 ,  and 
there  are  no  strongly  degenerate  cones  in  K(M).  For  d  ^  0,  L *(d)  is 
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the  class  of  all  matrices  M  where  K(Af)  has  no  strongly  degenerate  cones 
and  does  not  contain  d .  These  classes,  as  well  as  the  L(d) ,  were  introduced 
in  Garcia  (1973).  There  it  is  shown  if  M  £  L(d),  with  d  >  0,  then  for 
all  X  >  0,  we  have  {(Xd,  0)}  =  sol(Xd,  M) .  Hence,  for  d  >  0,  we  have 
L*(d)  =  L*(0)nL(d).  While  before  we  had  L  =  fldxjMd),  we  can  only  say 
here  that  E  C  rid>oL*(d).  For  example,  the  matrix 

0  1  ' 

.  — 1  1  . 

is  in  nd>0L*(d),  but  is  not  in  E.  We  will  have  more  to  say  about  these 
classes  later  on. 

(N)  A  matrix  M  £  5RnXn  is  said  to  be  in  N  if  and  only  if  all  principal 
minors  are  negative.  Two  standard  references  for  this  class  are  Saigal  (1972a), 
and  Kojima  and  Saigal  (1979).  More  will  be  said  about  this  class  in  what 
follows. 

(NflQ)  A  matrix  M  £  9inXn  is  in  this  class  if  and  only  if  it  is  both 
in  N  and  in  Q.  It  is  shown  in  Kojima  and  Saigal  (1979)  that  if  M  £  N, 
then  M  £  Q  if  and  only  if  M  0 .  It  is  also  shown  if  M  £  N  fl  Q ,  then 
|sol(g,  M) |  equals  1  for  q  ^  0 ,  and  equals  2  for  0  <£  q  >  0 .  According 
to  Theorem  3.3  of  Kojima  and  Saigal  (1979),  if  M  £  NnQ  and  q  >  0 
then  |sol(g,  Af)|  equals  3  if  all  solutions  to  ( q,M )  are  nondegenerate,  and 
equals  2  otherwise.  (Actually,  what  it  means  for  a  solution  of  ( q ,  M)  to  be 
“nondegenerate”  is  never  defined  in  that  paper,  however,  it  can  be  inferred 
from  context  and  the  cited  references  that  the  intended  definition  is  the 
one  given  here  in  Chapter  1.)  While  it  is  true  that  there  will  be  exactly 
three  solutions  for  all  q  >  0  having  only  nondegenerate  solutions,  it  is  true 
that  there  are  exactly  three  solutions  for  all  q  >  0.  The  last  line  of  the 
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proof,  given  in  Kojima  and  Saigal  (1979),  says  a  solution  is  “lost”  because  of 
degeneracy.  This  will  be  the  case  when  0  Q  >  0  and  it  is  on  a  reflecting 
face;  however,  for  q  >  0  we  are  only  contending  with  proper  faces,  and  no 
solutions  are  “lost.”  Consider 

Example  5.1  Let 


M  = 


4 

—  1 
—  1 


It  can  be  easily  checked  that  M  €  N  and,  as  M  <};  0,  M  €  Q.  Let 
q  =  (2, 4,  \)T  .  Then  ( q ,  M)  has  three  solutions 

(tu1,*1)  =  (2,4, 1,0, 0,0) 

(vj2,z2)  =(0,6, 5, 2, 0,0) 

(in3,*3)  =(3,0, 0,0, 0,1) 


and  (iw3,  z3)  is  degenerate. 

(N  \  Q)  A  matrix  M  €  3JnXn  is  in  this  class  if  and  only  if  M  is  in 
N  but  not  in  Q.  In  Kojima  and  Saigal  (1979)  it  is  shown  that  this  is  the  set 
of  matrices  M  £  N  for  which  M  <  0,  hence,  as  pointed  out  in  the  paper, 
we  will  have  K(M)  =  .  Therefore  these  matrices  are  in  Qo ,  hence  all 

of  N  is  in  Qo .  Kojima  and  Saigal  (1979)  also  shows,  for  all  q  >  0 ,  i.e,  all 
q  6  int  K(M)  ,  we  have  |sol(9,  M)\  =  2 .  This  means 


N  \  Q  C  INSa . 


((NnQ)- 1)  A  matrix  M  6  !RnXn  is  in  this  class  if  its  inverse  is 
in  NnQ.  This  is  equivalent  to  saying  the  matrix  is  in  Q,  its  determinant 
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is  negative,  and  all  of  its  proper  principal  submatrices  are  in  P.  Thus  all 
proper  principal  submatrices  are  in  Q,  along  with  the  matrix  itself.  Hence 
these  matrices  are  completely-Q,  which  is  to  say  they  are  in  E.  For  more  on 
these  matrices  see  Saigal  (1972b).  The  following  example  helps  to  justify  the 
placement  of  these  matrices  in  Figure  5.1.  Let 


0  ' 
4 

2  . 


Notice  M  1  G  N  n  Q .  However,  M  $£  CP  as  letting  x  =  (3, 2, 0)T  we  have 
xTMx  <  0. 


(P)  A  matrix  M  6  SnXn  is  in  the  class  P  if  and  only  if  all  principal 
minors  of  M  are  positive.  This  is  one  of  the  most  studied  classes  of  matrices 
related  to  the  LCP.  There  are  many  equivalent  characterizations  of  these 
matrices,  for  example:  M  G  P  if  and  only  if  for  all  <7  G  9?n  we  have 
|sol(<7,  M)|  =  1,  see  Samelson,  Thrall  and  Wesler  (1958),  and  Murty  (1972); 
M  G  P  if  and  only  if,  for  x  G  3?”,  we  have  Zj(Mz)<  <  0  for  all  *  6  n 
implies  x  =  0 ,  see  Fiedler  and  Ptak  (1962),  also  Gale  and  Nikaido  (1965); 
and  M  G  P  if  and  only  if,  for  A  G  3?n  x  n ,  we  have  det(/  —  A  -j-  AM)  ^  0 
for  all  0  <  A  <  /,  see  Aganagic  (1978).  The  middle  characterization  gives 
some  intuition  behind  the  definition  of  E,  as  it  states  a  matrix  belongs  to  P  if 
and  only  if  for  every  non-zero  x  G  3?n  >  (not  just  x  G  ),  we  have  an  index 
k  G  n  for  which  x *(Mi)fc  >  0.  An  interesting  characterization  by  Habetler 
and  Kostreva  (1980)  is  as  follows.  Say  a  point  x  G  is  a  complementary 
point  of  ( q ,  M )  if  and  only  if  there  is  a  :  G  K",  where  for  all  *  €  ft  we 
have  ( Mz  q)i  z+  =  Q ,  such  that  x  =  z  +  ( Mz  -f-  q).  It  is  then  the  case 
that  MG  P  if  and  only  if  there  is  some  q  E  9 Jn  such  that  the  interior  of 
each  orthant  in  5Rn  contains  exactly  one  complementary  point  of  (q,M). 


For  more  on  P-matrices,  see  the  references  mentioned  and  also  Fiedler  and 
Ptak  (1966). 

(Po)  A  matrix  M  €  5Rn><n  is  in  Po  if  and  only  if  all  principal  minors 
of  M  are  nonnegative.  Like  P,  this  class  has  been  extensively  studied.  In 
fact,  the  question  of  exactly  what  structure  and  properties  are  lost  when 
dealing  with  Po  as  opposed  to  dealing  with  P  was  one  of  the  questions 
leading  to  the  present  work,  and  to  other  works  in  the  field.  Again,  major 
references  to  this  class  are  the  papers  by  Fiedler  and  Ptak  (1962,  1966).  An 
interesting  characterization  of  Po ,  giving  insight  into  the  definition  of  Eo , 
comes  from  Fiedler  and  Ptak  (1966)  and  states  M  6  Po  if  and  only  if  for  all 
O^z  6  Sn,  (not  just  x  €  9?!). ),  we  have  an  index  k  6  n  for  which  ijt  7^ 0 
and  Xk(Mx)k  >  0.  We  move  on  to  a  special  class  of  Po -matrices  which  were 
defined  earlier  in  this  work. 

(Pi )  A  matrix  M  G  3RnXn  is  in  Pi  if  and  only  if  M  £  P0  and  exactly 
one  principal  minor  of  M  is  zero.  This  class  fits  into  Figure  5.1  in  about  the 
same  position  as  Po .  However,  we  do  know 

Theorem  5.2  Pi  CL. 

Before  starting  the  proof,  we  introduce  a  lemma. 

Lemma  5.3  If  M  €  Eo  f|9inXn  and  for  some  t'£n  we  have  M<  >  0 
with  Mu  =  0 ,  then  Mg  Q . 

Proof.  Suppose  we  have  a  matrix  M  satisfying  the  hypothesis  of  the  lemma. 
Take  some  q  €  pos  C(0). j .  By  reasoning  similar  to  previous  dimensional 
arguments,  we  may  assume  q  lies  in  no  A: -dimensional  faces  of  K(M),  for 
k  <  n  —  1 ,  and  any  (n  —  1) -dimensional  face  of  K(M)  that  contains  q  is 
contained  in  the  hyperplane 
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H  =  span  C(0).j  =  {z€9RTl:Xj  =  0}. 


Let  H +  be  the  closed  half-space  with  H  as  boundary  that  contains  I\ , 
and  let  H~  be  the  other  closed  half-space.  Let  posC(a).j  be  a  (n  —  1)- 
dimensional  face  of  K(M)  that  contains  M.i  and  is  contained  in  H .  We  see 
pos  (7(a)  j  cannot  contain  I.i  £  H .  Also,  the  vectors  of  C(a ).}  are  linearly 
independent,  as  pos(7(a)j  is  (n  —  1) -dimensional,  hence  the  face  cannot 
contain  — M.i .  Thus  j  —  i.  As  M\  E  pos  (7(0). j,  and  as  the  (n —  1)- 
dimensional  faces  of  K(M)  are  finite  in  number  and  closed,  we  can  select 
q  close  enough  to  M.i  such  that  we  have  the  additional  property  that  any 
(n —  1) -dimensional  face,  pos  C(a).j,  of  K (M)  that  contains  q  must  have 
j  =  i.  Now  for  all  c  >  0  small  enough,  B[q,  e)  D  K (M)  =  B{q,t)C\H ,  and 
no  face  of  K(M)  whose  dimension  is  smaller  than  n  —  1  intersects  B(q,  e) . 
Thus  B(q,  c)  fl  pos  C(0).j  =  B(q,  e)C\H .  Hence 

0  7^  B(q,  e)  fl  int  pos  (7(0)  C  H+  . 

Since  M  £  Eo  ,  no  other  full  cone  can  intersect  the  interior  of  pos  (7(0) .  Thus 
any  full  cone,  pos  C(a) ,  containing  B(q,  e)  fl  int  H~  must  have  a  boundary 
face  in  H.  This  face  will  then  contain  q,  and  so  this  face  is  posC(a).f. 
As  both  I.i  and  —M.i  are  in  H+ ,  then  we  have  pos  (7(a)  C  H+ ,  giv¬ 
ing  us  a  contradiction.  Hence  no  full  cone,  and  hence  no  cone,  contains 
B{q,e)nH-  Thus  Mi  Q. 

□ 

Proof  of  Theorem  5.2.  Let  M  6  Pi  ft  5?nXn .  We  know  M  €  Eo  as 
Pi  Q  ?o  Q  Eq  C  Eo  .  If  sol(0,  M )  =  { (0, 0) }  then  ME  L .  Thus  assume 
there  is  a  non-trivial  solution,  say  (w,  z)  with  z  0,  to  (0,  M).  Thus,  letting 
y  =  tv  -f*  z ,  we  have  for  some  a  E  (n)  that  ya  =  za ,  y&  =  w& ,  and 
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C(a)  y  =  0.  Thus  we  know  Maa  ya  —  0 .  In  addition,  we  must  have  ya  >  0 
for  otherwise  some  principal  submatrix  of  Maa ,  and  hence  of  M ,  is  singular, 
but  Maa  is  the  only  singular  principal  submatrix  of  M .  In  the  same  way, 
we  know  y&  >  0 .  Else,  for  some  t'6a,we  have  iyaya  =  0 .  Thus,  with 
P  —  a  U{  t } ,  we  have  Mpp  zp  =  0  which,  again,  contradicts  the  fact  that 
Maa  is  the  only  singular  principal  submatrix.  Thus  y  >  0. 

We  now  show  Af^Q.  If  jo)  =  1 ,  then  it  must  be  that,  for  some  i  E  n 
where  a  —  { i } ,  we  have  M.t  >  0  and  M»,  =  0 .  Thus,  by  Lemma  5.3, 
M  £  Q.  Suppose  |o|  >  1.  Pick  some  ft  C  a  with  \f3\  =  |o|  —  1,  and 
let  M  be  the  principal  transform  of  M  gotten  by  block  pivoting  on  Mpp . 
(Again,  we  know  Mpp  is  nonsingular  as  Maa  is  the  only  singular  principal 
submatrix.)  Since  we  have  Maa  Va  —  0  and  M&a  ya  >  0 ,  then,  letting 
{i}  =  a  \  (3 ,  we  have  Mu  >  0  and  Mu  —  0.  Since  M  6  Eq  ,  we  have 
JVf  €  E0  ,  thus  Lemma  5.3  gives  us  ~M  &  Q .  Hence,  as  claimed,  M  £  Q . 

From  Theorem  2.25,  M  €  U  and  K  (M)  is  a  half-space.  Let  9?n 

be  a  normal  to  the  hyperplane  d  K(M)  .  As  M  £  U ,  so  K(M)  is  regular, 
we  must  have  posC(o)  C  dK(M)  thus  C(a)T x  =  0.  Since  all  other 

complementary  cones  are  full,  they  cannot  be  contained  in  3K(M)  .  Thus 

C{a)T x  >  0.  Therefore  xa  >  0  and  z&  =  0.  Also,  xT M.a  =  0 

and  iTM.i  <  0.  Hence,  we  can  choose  x  so  that  ||xj|  is  so  small  that 

2  ^  >  0,  and  w  >  —  MTx  >  0.  This  means  M  satisfies  the  conditions 

to  be  in  L,  and  the  theorem  follows. 

□ 

It  should  be  noted  that  P*  CP ,  for  consider  the  matrix 
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Clearly,  M  €  Pi .  Yet,  with  x  =  (1,  l)r ,  we’ve  xT Mx  <  0. 

(Pi  \  Q)  A  matrix  M  €  S?nXn  is  in  this  “class”  if  and  only  if  it  is  in 
Pi ,  but  not  in  Q .  This  class  has  the  same  position  in  Figure  5.1  as  does  Pi , 
except  it  is  also  contained  in  U.  More  was  said  about  Pi  \  Q  at  the  end  of 
Chapter  2. 

(PD)  A  matrix  M  6  SRnXn  is  said  to  be  positive  definite,  M  6  PD,  if 
and  only  if  for  all  0  x  6  SRn  we  have  xT  Mx  >  0 .  For  symmetric  matrices, 
being  in  P  is  equivalent  to  being  in  PD,  which  is  equivalent  to  there  being 
some  L  €  3fnXn  such  that  L  is  nonsingular  and  M  =  LTL.  For  more 
concerning  positive  definite  matrices,  see  Gantmacher  (1960),  Dantzig  and 
Cottle  (1967),  Cottle,  Habetler  and  Lemke  (1970a),  and  Cottle  (1983). 

(PSD)  A  matrix  M  6  5RnXn  is  said  to  be  positive  semi- definite, 
M  6  PSD ,  if  and  only  if  for  all  x  €  5Rn  we  have  xT  Mx  >  0 .  For  symmetric 
matrices,  being  in  P0  is  equivalent  to  being  in  PSD,  which  is  equivalent 
to  there  being  some  L  €  5RnXn  such  that  M  —  LTL.  The  class  PSD  is 
usually  thought  of  in  connection  with  convexity  as  the  quadratic  function 
F(x)  :  5Rn  — >  5R  defined  by  F( x)  =  xT Mx  -J-  cTx  +  d,  with  M  €  lRnXn, 
c  €  3Rn  and  d  6  9?,  is  convex  if  and  only  if  M  6  PSD.  See  the  references 
given  for  positive  definite  matrices. 

(Q)  A  matrix  M  €  9?nXn  is  said  to  be  in  Q  if  and  only  if  for  all 
q  €  5Rn  the  LCP  ( q ,  M)  has  at  least  one  solution.  This  is  equivalent  to  saying 
K(M)  =  SRn .  One  of  the  major,  and  perhaps  most  difficult,  problems  in 
linear  complementarity  theory  is  to  find  a  “good”  characterization  of  Q,  i.e., 
a  characterization  with  which  one  could  quickly  test  a  matrix  to  determine 
whether  or  not  it  is  in  Q.  Many  of  these  other  matrix  classes  were  studied  in 
attempts  to  find  more  classes  of  matrices  that  were  contained  in  Q,  or  Qo . 
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Two  interesting  works  concerning  Q  are  Watson  (1974),  Kelly  and  Watson 
(1979).  The  latter  contains  a  counterexample  to  a  result  of  the  former.  In 
essence,  it  shows  the  annoying  result  that  the  set  of  Q-matrices  is  neither 
open  nor  closed  in  SRnXn  for  n  >  4.  Hence,  the  class  Q  will  be  nard  to 
characterize.  See  also  Cottle,  von  Randow  and  Stone  (1981). 

(Qo)  A  matrix  M  €  3?nXn  is  said  to  be  in  Qo  if  and  only  if  for  all 
q  6  3tn  where  the  system  of  inequalities 

Mz  4-  q  >  0  z  >  0 

is  feasible,  there  exists  at  least  one  solution  to  the  LCP  [q,  M ) .  This  is 
equivalent  to  saying  K(M)  is  convex.  Like  Q,  characterizing  Q0  in  a  “good” 
way  is  a  long  standing  problem.  In  fact,  with  a  characterization  of  this  class 
we  can  just  say  Q  =  Qo  D  S .  Again,  many  of  the  works  mentioned  are 
concerned  with  Q0 .  For  a  recent  and  interesting  paper  on  this  class  see 
Doverspike  and  Lemke  (1981).  (In  other  works,  this  class  is  denoted  K;  it 
should  not  be  confused  with  the  K  used  here.) 

(R)  A  matrix  M  6  9?nXn  is  said  to  be  regular,  M  £  R,  if  and  only 
if,  with  e  =  (1,1,...,  l)T  €  9?n,  we  have,  for  all  X  >  0,  that  { (Xe,  0)}  = 
sol(Xe,  M).  Clearly,  R  =  L*(e).  The  standard  reference  for  this  class  is 
Karamardian  (1972).  It  is  of  interest  to  note,  as  shown  by  Aganagic  and 
Cottle  (1978),  that  Po  flQ  =  Po  HR.  We  cannot  do  better  than  this  in 
classifying  Po  D  Q  as  far  as  Figure  5.1  is  concerned.  For  example,  the  matrix 

0  1  ’ 

.  — 1  1  . 

is  in  Q  and  in  Po  but  is  not  in  E  which  is  the  next  matrix  class  “lower”  in 
Figure  5.1  than  R. 
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(S)  A  matrix  M  G  9JnXn  is  said  to  be  in  S  if  and  only  if  there  exists 
an  x  G  9Jn  such  that  x  >  0  and  Mx  >  9.  This  is  the  class  of  matrices 
for  which  ( q ,  M )  is  “feasible”  for  all  q  G  SRn  ,  i.e.,  for  all  q  G  SRn  there  is  an 
xg  G  !Rn  with  Xg  >  0  and  Mxq  -f-  q  >  0,  see  Lemke  (1970).  The  classic 
reference  for  these  matrices  is  Fiedler  and  Ptak  (1966).  Other  relevant  works 
to  look  at,  that  use  S-matrices  are  Saigal  (1971a)  and  Cottle  (1979). 

(So)  A  matrix  M  G  9Jn><n  is  said  to  be  in  So  if  and  only  if  there 
exists  an  i  £  S"  such  that  0  ^  x  >  0  and  Mx  >  0 .  This  is  clearly 
one  of  the  largest  matrix  classes  listed  here,  containing  many  of  the  others. 
Again,  Fiedler  and  Ptak  (1966),  Lemke  (1970),  and  Saigal  (1971a)  are  good 
references  for  this  class.  For  a  nice  reference  which  extends  the  properties 
embodied  in  the  matrix  classes  P,  Po ,  S,  and  So  to  non-linear  functions,  see 
More  and  Rhcinboldt  (1973). 

Two  of  the  inclusion  arrows  leading  to  the  class  So  in  Figure  5.1  are 
not  trivial,  and  have  not  been  found  by  the  author  in  the  literature.  The 
justification  for  these  inclusions  is  in  the  following  two  theorems. 

Theorem  5.4  U<i>0L(d)  C  S0  . 

Proof.  Suppose  for  some  0  <  d  £  we  have  M  G  L(d)n3ftn><n,  but 
M  So  .  If  (w,  z)  G  sol(0,  M),  then  .Mz  >  0  and  z  >  0.  Thus  M  £  So 
implies  z  =  0.  Hence  {(0,0)}  =  sol(0,  M).  Garcia  (1973)  shows  that 
M  G  L(d)  implies  for  all  X  >  0  that  we  have  { (Xd,  0) }  =  sol(Xd,  M) . 
We  conclude  M  G  L*(d).  But  L*(d)  C  Q  C  S  C  So  ,  which  gives  us  a 
contradiction.  Thus  MG  So  and  the  theorem  holds. 

□ 

Theorem  5.5  E0  C  S0  . 
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Proof.  Suppose  M  E  Eo  D9in><n-  Let  /  E  3RnXn  be  the  identity  matrix. 
Thus,  for  all  e  >  0,  we  have  almost  directly  from  the  definitions  that 
M  -f-  *1  €  E .  Now  E  C  Q  C  S  C  So  ,  so  for  each  e  >  0  there  is  a 
0f^xe  >0  such  that  (M  -j-  e/)x£  >  0 .  We  may  assume,  by  scaling,  that 
||x£||  =  1 .  As  the  set 

{i£  dln  :  ||x||  =  1 },  the  unit  sphere  in  3ffn , 

is  compact,  we  have  some  point  xo  €  5Rn ,  with  )|xo||  =  1 ,  that  is  a  cluster 
point  of  the  set  of  x£ .  Thus,  letting  e  — ►  0 ,  we  see  that  xo  >  0 ,  and  that 
Mx o  >  0 .  Thus  ME  So  ,  and  the  theorem  holds. 

□ 

(SCP)  A  matrix  M  E%lnXn  is  said  to  be  strictly  copositive,  M  E  SCP , 
if  and  only  if  xT  Mx  >  0  for  all  x  6  such  that  Oy^x  >  0 .  This  class  has 
also  been  denoted  as  C.  For  symmetric  matrices,  being  in  SCP  is  equivalent 
to  being  in  E.  See  the  references  given  for  copositive  matrices. 

(U)  A  matrix  M  E  3?nXn  is  said  to  be  in  U,  for  Unique  solution,  if 
and  only  if  for  all  q  E  int  K (M)  we  have  |sol(g,  M)\  —  1 .  This  matrix 
class  was  the  topic  of  Chapter  2.  If  M  €  Eq  ,  then  for  all  q  >  0  we  have 
|sol(<7,  M)\  =  1.  Garcia  (1973)  shows  that  if  ME  L  (d),  for  some  d  >  0, 
then  |sol(d,  M)\  =  1 .  Hence  we  see, 

E0  DINS  =  U  (  U  L(d) 
ld>0 

Thus,  as  U  D  Q  =  P ,  we  have 

(  U  L*(d)lnINS  =  P.  (5.8) 

U>o  ) 


jnlNSCU.  (5.7) 
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This  helps  us  understand  how  U  fits  into  Figure  5.1.  Now  consider  the 
following  matrices 

'  0  1 
.1  0  . 

(5.9) 

Notice  (5.9)  is  contained  in  U,  but  not  in  Qq  .  (This  is  Example  2.3.)  Hence, 
the  right  side  of  (5.7)  is  a  proper  inclusion.  As  for  (5.10),  it  is  not  in  U,  yet 
it  is  in  AnPSDnSCPnEflPi  •  Also,  (5.6)  showed  an  example  of  a  matrix 
that  was  in  BG  but  not  in  Eq ,  hence  certainly  not  in  U. 

(Z)  A  matrix  M  g  5JnXn  is  said  to  be  in  Z  if  and  only  if  for  all  i,j  6  n, 
where  j ,  we  have  Mij  <  0 .  These  matrices  have  been  well  studied.  See, 
for  example,  Saigal  (1971b)  and  the  references  mentioned  in  the  paragraphs 
concerning  the  classes  K  and  Ko .  In  particular,  see  Mohan  (1978). 

(Z)  A  matrix  M  €  dtnXn  is  said  to  be  in  Z  ir  and  only  if  M  €  Z 
and,  for  all  i  £  n,  we  have  Mu  >  0.  (This  class  has  also  been  denoted  by 
L.)  Sec  the  references  for  Z-matrices  and,  in  particular,  see  Saigal  (1972b) 
and  Mohan  (1978).  One  thing  that  should  be  pointed  out  is  an  error  in 
Theorem  5.4  of  Saigal  (1972b).  The  theorem  states  that  M  €  Z  implies 
K(M)  is  regular.  (Saigal’s  definition  of  regularity  and  the  definition  used 
in  this  work  are  different,  however,  all  that  need  be  known  here  is  that  the 
two  definitions  coincide  for  nondegenerate  matrices.)  This  is  not  the  case. 
Consider 
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M  = 


This  matrix  is  nondegenerate  and  contained  in  Z.  Let  q  =  (51, 11, 19, 39)T 
and  q  —  (49, 9,  21, 41)r  .  Then  (q,  M)  and  (q,  M )  both  only  have  non¬ 
degenerate  solutions,  but  Jsol(<7,  M)|  =  4  and  |sol(9,M)|  =  2.  Specifically, 
the  solutions,  ( w ,  z ) ,  of  ( q ,  M)  are 

(51,11,19, 39,0, 0, 0, 0) 

(0,0, 6^,0, 3^,  18^,0,  §|) 

while  the  solutions  of  (q,  M )  are 

(49,9,21,41,0,0,0,0)  (0,  0, 11  2$,3£,  17$,  0,  0)  . 

This  implies  M  INS  and,  as  M  is  nondegenerate,  that  K(M)  is  not 
regular  by  either  definition.  This  incorrect  result  is  referred  to  by  Mohan 
(1978)  in  several  places.  Saigal  (1972b)  uses  it  to  “show”  that  if  M  €  Z  then 
M  €  INS2 ,  which  is  clearly  not  true  as  seen  in  the  example  just  given.  More 
will  be  said  about  this  in  the  next  section. 


(4, 0, 0, 0, 0, 8§,  10  J,  9§) 
(0,0,0,0,1J,13§,5§,6$), 


5.2  Related  LCP  Theory 

In  this  section  we  will  consider  some  results  in  the  LCP  literature  that 
seem  related  to  the  material  we  have  covered.  Most  of  the  results  concerning 
the  exact  number  of  solutions  to  the  LCP  have  already  been  mentioned. 

There  is  the  classic  result  of  Samelson,  Thrall  and  Wesler  (1958)  that  P 
is  the  set  of  all  matrices  M  such  that  for  all  q  we  have  jsol(g,  M)\  =  1. 

In  Eaves  (1971),  it  is  shown  if  M  €  Po  and  9  is  contained  in  the 
interior  of  some  full  complementary  cone  then  |sol(g,  M) |  =  1 .  This  can 
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be  seen  to  follow  from  the  fact  that  Po  C  Eg  .  For  it  is  easily  shown 
that  E q  can  be  characterized  as  the  set  of  those  matrices  M  such  that, 
for  all  q,  if  q  is  contained  in  the  interior  of  a  full  complementary  cone  then 
|sol(g,  M) |  =  1 .  Related  to  this  is  Theorem  2.2  in  Saigal  (1970a)  which  states 
that  ME  Po  implies  that  if  sol(g,  M)  contains  a  nondegenerate  solution 
then  |sol(g,  M)\  =  1.  This  can  be  seen  to  be  in  error  by  considering  the 
matrix  [0]  =  M  E  Po  fl  SR1*1 .  We  have  for  X  >  0  that  (w,  z)  =  (0,  X) 
is  a  nondegenerate  solution  to  (0,  M) .  It  should  also  be  mentioned  that 
this  result  generalized  a  previous  result  of  Lemke  (1965)  which  used  positive 
semi-definite  matrices,  a  smaller  subclass  of  Po . 

There  are  several  theorems  by  Murty  (1972)  on  this  subject,  including 
another  proof  of  the  Samelson,  Thrall  and  Wesler  result.  The  main  theorems 
from  Murty  (1972)  of  interest  here  are:  |sol(<?,  M)|  <  oo  for  all  q  if  and  only 
if  M  is  nondegenerate;  if  |sol(?,  M) j  is  constant  over  all  non-zero  q ,  then 
that  constant  is  one  and  M  E  P ;  if  |sol(g,  M)\  is  constant  for  all  q  which 
are  nondegenerate  with  respect  to  M ,  then  that  constant  is  one. 

The  class  N  is  studied  in  Kojima  and  Saigal  (1979).  It  is  shown  that  for 
ME  N,  if  M  0  then  the  value  of  |sol(g,  M)|  is  one  for  q  ^  0,  is  two 
for  0  q  >  0,  and  is  three  for  q  >  0  nondegenerate  with  respect  to  M . 
This  last  part,  as  noted  earlier,  should  state  |sol(q,  M)\  =  3  for  all  q  >  0. 
(Kojima  and  Saigal  (1979)  incorrectly  state  the  value  is  two  for  q  >  0  but 
degenerate  with  respect  to  M .)  It  is  also  shown  for  M  E  N ,  if  M  <  0  then 
the  value  of  jsol(g,  M)\  is  zero  for  q  ^  0,  is  one  for  0  q  >  0,  and  is  two 
for  q  >  0. 

In  Mohan  (1980),  it  is  shown  for  M  E  Ko  ,  q  €  int  K(M)  implies 
|sol($, M)|  =  1  and  q  E  3K (M)  implies  |sol(g, M)\  =  oo. 
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In  Saigal  (1972b),  the  concept  of  a  “regular  pseudomanifold”  was  dis¬ 
cussed  with  reference  to  K (M).  K(M)  was  defined  there  as  being  “regular” 
if  and  only  if  every  face  was  either  proper  or  contained  in  d  K(M)  .  However, 
the  definition  of  “proper”  given  there  is  different  from  what  is  used  here- 
A  face  was  defined  there  as  being  proper  if  and  only  if  it  is  the  intersec¬ 
tion  of  the  two  adjacent  complementary  cones  containing  it.  (It  is  clearly 
in  the  intersection.  The  requirement  here  is  that  the  intersection  contain 
nothing  else.)  “Proper”  in  our  sense  implies  “proper”  in  Saigal’s,  but  not 
conversely.  For  instance,  if  a  full  cone  is  adjacent  to  a  degenerate  cone  the 
common  face  would  be.  considered  “proper”  by  Saigal’s  definition  but  not  by 
ours.  Hence,  our  definition  of  K(M)  being  regular  is  strictly  stronger.  We 
will  use  the  italic  proper  and  regular  to  refer  to  Saigal’s  (1972b)  definition, 
and  standard  lettering  for  our  own  definitions.  Notice  the  definitions  are 
equivalent  for  nondegenerate  matrices.  As  pointed  out  before,  Saigal  (1972b) 
incorrectly  “proves”  that  ME  Z  implies  K(M)  is  regular.  An  example  of 
a  Z  -matrix  where  K(M)  is  neither  regular  or  regular  was  given  in  the  last 
section.  However,  the  paper  also  contained  the  “theorem”  that  if  K(M) 
is  regular,  M  P ,  M  is  nondegenerate  and  soI(g,  M )  contains  only  non¬ 
degenerate  solutions,  then  |sol(g,  M)j  =  2.  This  is  also  incorrect.  For  ex¬ 
ample,  letting  M  be  the  negative  of  the  identity  matrix  in  5R2X2  we  have 
M  is  non  degenerate,  K  (M)  is  regular  and  hence  is  regular,  M  £  P  and  yet 
q  E  int  K(M)  implies  |sol(g,  M)\  —  4.  A  possible  substitute  here  could  be 
gotten  from  Corollary  4.6  which  would  state  that  if  M  is  non  degenerate, 
M  £  P,  K(M)  is  regular  (so  M  6  INS  by  Corollary  3.18),  then  if  sol(g,  M) 
contains  only  nondegenerate  solutions,  then  |sol(g,  M )|  is  even.  These  two 
errors  in  Saigal  (1972b)  cause  some  results  of  Mohan  (1978),  which  depend  on 
them,  to  be  incorrect.  These  results  are  Theorems  1.3.8,  1.5.8,  1.5.12,  3.3.3, 
3.3.4,  and  Corollary  3.3.1  of  Mohan  (1978). 
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Aside  from  questions  concerning  the  exact  number  of  solutions,  another 
concept  that  has  been  studied  is  the  constant  parity  property.  We  say  a 
matrix  M  has  the  constant  parity  property  if  and  only  if  the  parity  of 
|sol(q,  M) | ,  i.e.,  whether  it  is  even  or  odd,  is  the  same  over  all  q  where 
sol(q,  M)  contains  no  degenerate  solutions.  (Thus  if  M  ^  Q  and  has  the 
constant  parity  property  then  the  parity  is  even.  Given  any  q  K(M)  we 
have  sol(?,  M )  =  0  contains  no  degenerate  solutions  and  has  even  parity.) 
The  concept  of  constant  parity  is  a  weaker  form  of  the  concept  of  a  con¬ 
stant  number  of  solutions.  Clearly  all  INS -matrices  have  the  constant  parity 
property. 

The  classic  theorem  on  constant  parity  was  shown  by  Murty  (1972).  It 
states  that  a  nondegenerate  matrix  has  the  constant  parity  property.  Also 
in  this  paper  is  the  theorem  that  a  nonnegative  Q-matrix  has  the  constant 
parity  property  with  the  parity  being  odd. 

In  Saigal  (1970b)  we  find  the  following  theorem  on  the  constant  parity 
property:  If  — MT  £  S,  then  M  has  the  constant  parity  property  with  the 
parity  being  even.  The  final  word  on  the  subject  was  in  essence  given  by  Saigal 
(1972a).  It  states  that  a  matrix,  M  £  3}nXn ,  has  the  constant  parity  property 
if  and  only  if  it  is  true  that  for  any  collection  pos  C^aj),  pos  C(a2)»  ••• » 
pos  C(afc)  of  strongly  degenerate  complementary  cones,  where  k  is  odd  and 


dimfpos  C(ai)  fl  •  •  •  fl  pos  C(ajt)J  =  n  —  1 , 


there  exists  for  each  q  in  this  intersection  another  strongly  degenerate  com¬ 
plementary  cone,  pos  C(a*^_i) ,  such  that  q  E  pos  C(afc-^i) .  This  result 
expresses  the  basic  geometric  structure  behind  the  constant  parity  property. 
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Mohan  (1978)  proves  the  following  related  theorem  concerning  Z-matri- 
ces:  If  M  £  ZnSRnXn  and  there  is  a  x  6  3ftn  such  that  MTx  >  0,  then  M 
has  the  constant  parity  property  and  the  parity  is  odd  if  and  only  if  ME  K . 

The  last  area  of  complementarity  theory  we  will  bring  up  is  Lemke’s 
algorithm.  An  algorithm  for  solving  the  LCP  was  suggested  by  Lemkc  and 
Howson  (1964),  and  Lemke  (1965).  It  has  since  become  a  major  tool  in 
the  field,  inspiring  much  research  into  other  algorithms  based  on  the  same 
principles  and  causing  many  studies  to  determine  conditions  for  which  the 
algorithm  is  guaranteed  to  “process”  a  given  LCP.  For  a  detailed  description 
of  Lemke’s  algorithm  see  the  two  papers  mentioned  or  see  Eaves  (1971)  or 
Cottle  (1983).  The  essential  concept  is  as  follows.  Given  the  LCP  ( q ,  M ) ,  we 
take  some  vector  0  <  d  E  9?”  and  consider  the  family  of  LCPs  ( q  -f-  Od,  M), 
where  the  parameter  6  is  taken  as  a  nonnegative  number.  (In  the  canonical 
statement  of  the  algorithm,  d  is  taken  to  be  (1, 1, ... ,  1)T  .)  For  all  9  large 
enough  we  will  have  q  -f-  Od  >  0  and  hence  ( q  +  Od,  0)  £  sol (q  Od,  M) .  In 
other  words,  the  “tail”  of  the  ray 

{q  +  9d\0>0}  (5.11) 

is  contained  in  the  positive  quadrant.  We  then  move  back  along  the  ray  (5.11) 
attempting  to  get  to  q .  When  we  reach  the  face  of  a  complementary  cone  we 
continue  in  the  adjacent  cone.  Thus  a  proper  face  allows  us  to  travel  in  the 
same  direction  along  the  ray  (5.11)  as  we  had  been  traveling,  while  a  reflecting 
face  causes  us  to  change  direction.  The  problems  associated  with  reaching  a 
degenerate  face,  or  with  reaching  a  nondegencrate  face  on  its  boundary,  can 
be  taken  care  of  by  lexicographical  methods.  Again,  see  Eaves  (1971).  The 
actual  algorithm  is  carried  out  by  a  pivoting  scheme  which  gives  us  a  solution 
to  the  LCP  ( q  -f-  9d ,  M )  when  we  are  at  the  point  q  -j-  0d  of  the  ray  (5.11). 
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(This  solution  is  associated  with  the  complementary  cone  through  which  we 
are  currently  traveling.)  The  hope  is  we  eventually  reach  the  end-point  of  the 
ray  (5.11)  and  thus  find  a  solution  to  the  original  LCP  ( q,M ).  The  other 
two  possiblities  are  that  we  go  off  on  the  infinite  end  of  the  ray  (5.11)  never 
to  return,  or  we  reach  a  degenerate  face  with  no  other  full  cone  to  travel 
through  than  the  one  by  which  we  arrived.  It  is  now  clear  for  Mg  INS 
where  K(M)  is  star-shaped  at  d  >  0  that  we  will  find: 

1)  Lemke’s  method  using  d  will  process  ( q ,  M)  for  all  q  g  . 

2)  If  a  solution  is  found,  then  9  will  haye  been  monotonically  decreasing. 
That  is,  after  the  first  pivot  to  initialize  the  algorithm,  each  pivot  will 
cause  6  to  be  strictly  smaller.  (Actually,  to  prevent  degeneracy,  we  use 
lexicographical  techniques.  In  this  case  the  vector  used  in  place  of  0  is 
lexicographically  decreasing.) 

3)  If  when  running  the  algorithm  we  find  that  0  increases,  or  that  we  reach 
a  degenerate  face,  we  may  conclude  ( q ,  M )  has  no  solution. 

While  it  is  necessary  K(M)  be  star-shaped  at  d  >  0  for  these  conditions 
to  hold,  it  is  not  necessary  that  M  belong  to  INS.  For  the  matrix  (5.10), 
K(M)  is  star-shaped  at  d  —  (1,1)T,  and  the  above  three  conditions  hold. 
However,  (5.10)  is  not  in  INS. 

These  observations,  stated  with  a  different  emphasis,  are  basically  seen 
in  Theorem  4.1  of  Saigal  (1972b).  This  theorem  states  that  if  K(M)  is 
regular  and  contains  no  strongly  degenerate  cones,  then  a  necessary  and 
sufficient  condition  for  Lemke’s  algorithm  to  solve  (q,  M)  using  d  >  0  for 
all  q  g  K(M)  is  that  K (M)  be  star-shaped  at  d.  In  addition,  the  theorem 
states  that  9  will  be  monotonically  nonincreasing.  As  pointed  out  above, 
we  may  replace  ” regular ”  in  this  theorem  by  “regular.”  In  this  case,  the 
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condition  that  K(M)  contain  no  strongly  degenerate  cones  can  be  dropped. 
It  is  interesting  to  note  that  the  theorem  is  false  in  one  direction.  While  the 
star-shapedness  is  certainly  necessary,  it  is  not  sufficient.  Let 

‘  0—1  0 

M  —  — 1  0  1 

.  0  0—1 

We  find  K(M)  =  pos — M.  Also,  K(M)  contains  no  strongly  degenerate 
cones.  Notice  that  all  faces  of  all  complementary  cones  are  contained  in 
9K (M)  ,  except  for  pos  (7(0). 2  and  pos  C({  2  }).j .  However,  the  complemen¬ 
tary  cones  adjacent  on  pos  (7(0). 3  are  the  full  cone  pos  C(0)  and  the  degen¬ 
erate  cone  posC({2}).  Thus  pos  (7(0). 2  is  proper,  but  not  proper.  Similarly, 
the  cones  adjacent  on  pos  C({  2  }).  j  are  the  degenerate  cone  pos  (7({  2  })  and 
the  full  cone  pos  C{{  1, 2  }) .  Thus  K(M)  is  regular,  but  not  regular.  It  is  cer¬ 
tainly  star-shaped  at  d—  (1, 1,  l)r  .  Yet,  while  q  —  (1,-1, 2)  is  contained 
in  K(M),  in  fact  we  have  (0,  0,  0, 1, 1,  2)  6  sol(9,  M),  Lcmke’s  algorithm 
finds  no  solution  to  ( q ,  M )  using  any  d  >  0 .  Thus  the  sufficiency  part  of 
Theorem  4.1  in  Saigal  (1972b)  is  in  error. 

One  more  point  before  finishing  this  chapter  is  that  Saigal  (1972b)  defines 
K(M)  to  be  the  union  of  all  complementary  cones  of  dimension  n  —  1  or 
greater,  where  M  €  SRnXn.  While  this  is  often  the  case,  it  is  not  always 
true.  For  example,  let 

‘  0  0  1  0 

0  0  0  1 

M  = 

0  0  0  0 

.0  0  0  0 

140 


Now  q  =  ( — 1, — 1,0, 0)T  6  K(M)  .  Yet,  if  q  £  pos  C(a)  then  {3,4}  C  a 
and  hence  dimfpos  C(a) J  <  3 .  The  inclusion  relationships  of  the  matrix 
classes  discussed  in  this  chapter,  are  diagrammed  in  the  following  figure. 


CHAPTER  6. 
CONCLUSION 


The  central  emphasis  of  this  work  has  been  on  the  underlying  geometric 
stucture  of  LCP’s  with  the  global  property  of  an  invariant  number  of  solu¬ 
tions.  There  are  other  interesting  open  questions  related  to  this,  and  to  LCP 
theory  in  general.  It  seems  appropriate  to  mention  some  of  these  questions 
as  a  conclusion  to  this  study. 

Theorem  2.22  shows  that  Qo  flU  C  P0  .  In  essence,  if  we  think  of 
starting  in  the  positive  orthant,  which  is  a  positive  complementary  cone,  and 
“moving”  in  K(M)  through  a  sequence  of  adjacent  complementary  cones 
then,  if  M  (=  Qo  HU,  every  common  face  we  encounter  between  two  com¬ 
plementary  cones  is  proper,  until  we  reach  a  degenerate  face  which  must  be 
on  the  boundary.  Since  "reflecting”  isn’t  allowed,  as  those  type  of  faces  are 
forbidden  by  the  fact  that  M£U,  and  since  int  K(M)  is  path  connected  so 
we  can  reach  all  the  complementary  cones,  then  we  can  never  reach  a  negative 
cone.  (There  isn’t  enough  “room,”  and  there  are  too  many  restrictions,  to 
allow  us  to  “turn  around.”)  It  seems  that  there  isn’t  enough  “room”  even 
allowing  degenerate  faces  within  int  K(JW) .  Thus  a  problem  left  open  by 
this  study  is  to  determine  whether  or  not  Qo  ft  Eq  C  Pq  . 


Look  once  again  at  the  map,  F ,  used  in  the  proof  of  Theorem  3.15. 


If,  as  before,  we  assume  that  no  complementary  cone  of  K(M)  is  strongly 
degenerate,  then  we  can  associate  with  F,  and  hence  with  M,  a  special 
integer  referred  to  as  the  degree  of  F  (of  M).  Let  q  €  5Rn  be  any  vector 
that  is  nondegenerate  with  respect  to  M .  Then  the  degree  of  F  (of  M )  is 
the  number  of  positive  complementary  cones  containing  q  minus  the  number 
of  negative  complementary  cones  containing  q.  (It  can  be  shown  that  this 
number  will  be  invariant  over  all  q  nondegencrate  with  respect  to  M .  For 
more  on  the  concept  of  degree  see  Ortega  and  Rheinboldt  (1970).)  The  degree 
of  a  map  is  a  measure  of  the  number  of  points  in  the  domain  which  are  mapped 
to  each  point  in  the  range.  For  a  general  map  of  degree  k ,  however,  it  is 
not  necessary  that  any  point  in  the  range  have  exactly  |fc|  points  mapping 
into  it  from  the  domain.  However,  the  map  F  associated  with  an  LCP  is 
not  a  general  map.  Perhaps  it  is  the  case  for  these  special  maps,  that  when 
the  degree  of  F  is  k,  one  can  always  find  a  point  in  the  range  which  has 
exactly  |fc|  points  of  the  domain  mapping  into  it.  In  the  case  fc  =  0,  this 
would  mean  that  every  matrix  M  with  zero  degree  is  not  in  Q,  i.e.,  some 
point,  q ,  in  the  range,  9?n ,  of  F  has  no  point  in  the  domain  mapping  into 
it.  (Note  that  q  would  then  trivially  be  nondegenerate  with  respect  to  M .) 
This  is  not  the  case.  Kelly  and  Watson  (1979)  show  that  the  nondegenerate 
matrix 


'  21 

25 

—27 

-36  ' 

M  = 

7 

3 

—9 

36 

12 

12 

-20 

0 

-  4 

4 

—4 

—8  . 

is  in  Q,  yet  it  is  a  straightforward  calculation  to  verify  that  the  degree  of 
M  is  zero.  For  the  case  k  ^  0  the  question  is  still  open,  and  there  are 
reasons  to  believe  that  the  geometric  structure  of  non- zero  degree  matrices  is 
significantly  different  from  the  geometric  structure  of  zero  degree  matrices. 


144 


Hence  we  have  the  deep  question  in  LCP  theory  of  determining  whether  or 
not  there  exists  a  matrix  M  with  no  strongly  degenerate  complementary 
cones,  and  with  degree  k  such  that  |sol(g,  M)\  >  |fc|  whenever  q  is 
nondegenerate  with  respect  to  M.  Another  way  of  phrasing  this  is  to  ask 
if,  for  matrices  M  with  no  strongly  degenerate  complementary  cones,  it  is 
true  that  when  the  union  of  the  positive  complementary  cones  is  S?n  and 
the  union  of  the  negative  complementary  cones  is  ,  then  every  vector 
q,  nondegenerate  with  respect  to  M ,  is  contained  in  the  same  number  of 
positive  complementary  cones  as  negative  complementary  cones.  Indeed,  this 
is  conjectured  to  be  true  in  Garcia  and  Gould  (1980).  See  also  Howe  (1980). 
(It  should  be  pointed  out  that  the  class  “Qo”  in  Garcia  and  Gould  (1980)  is 
not  the  same  as  the  one  discussed  in  the  present  work.) 

At  the  end  of  Chapter  3  we  showed  that,  for  nondegenerate  matrices  M , 
K(M)  is  regular  if  and  only  if  M  6  INS .  It  seems  that  the  nondegeneracy 
assumption  should  be  unnecessary;  this  raises  the  question  of  whether  it  is  in 
general  true  that  K(M)  is  regular  if  and  only  if  M  £  INS . 

In  Chapters  3  and  4  we  developed  the  idea  of  the  partition  E  of 
3in  \  K(M)  .  We  noted  that  the  elements  of  E  are  not  in  genjral  convex,  not 
even  when  considering  only  those  elements  contained  in  K(M) .  The  question 
then  arises  as  to  whether  the  elements  of  E  are,  in  general,  star-shaped.  This 
question  is  open,  as  is  the  related  question  of  whether  there  will  always  be  an 
element  of  E  which  is  convex.  (Is  Theorem  4.4  valid  for  degenerate  matrices 
as  well  as  nondegenerate  matrices?) 

In  Chapter  4  we  already  have  discussed  Conjecture  4.8,  but  have  not 
spoken  about  Assumption  4.17.  This  is  a  technical  assumption  that  has 
been  used  in  another  form  by  other  authors.  The  last  open  question  we’ll 
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consider  is  the  one  surrounding  this  assumption  on  the  geometry  of  LCP’s. 
It  can  be  stated  as  follows.  Given  an  LCP  consider  the  related  map  F  as 
defined  in  the  proof  of  Theorem  3.15.  Suppose  we  let  D  be  the  union  of 
some  collection  of  orthants  in  9in  such  that  D  forms  a  pseudomanifold, 
i.e.,  between  any  two  orthants  in  D  there  is  a  path,  in  D ,  of  “neighboring” 
orthants.  The  image  under  F  of  each  orthant  is  a  complementary  cone. 
If  the  complementary  cones  which  are  the  images  of  the  orthants  in  D  are 
all  positive  complementary  cones  is  it  then  the  case  that  the  restricted  map 
F  :  D  — *•  is  injective?  If  D  =  the  answer  is  “yes”  as  shown  in 

Murty  (1972).  If  D  is.  convex  we  can  reduce  the  problem  to  the  case  where 
D  is  3Jm  for  some  m  <  n  and  the  answer  is  again  “yes”  by  the  result  in 
Murty  (1972).  In  general  the  question  is  open.  It  should  be  noted  that  the 
LCP  structure  is  important.  If  we  were  to  require  the  function  F  on  D  to 
just  be  piecewise-linear,  with  the  pieces  of  linearity  being  the  orthants,  and 
the  determinants  of  the  matrices  defining  the  affine  functions  on  adjacent 
orthants  to  be  of  opposite  signs,  then  F  :  D  — *  3?n  would  not  necessarily  be 
injective.  As  an  example,  consider  D  =  SR3 ,  and  F  defined  as  follows  on 
the  different  orthants 


F(xi,x  2,x3) 


(xi, 12.13) 

(xi,  Xi  +  X2,  Xi  +  X3) 
'  ( Xi  -hx2,x2,x2  +X3) 
(xi  +  X3,X2  +  13,13) 
F(-x) 


if  xi,  X2,  X3  >  0 
if  xi  <  0,  X2  >  0,  x3  >  0 

if  xx  >  0,  X2  <  0,  X3  >  0 

if  Xi  >  0,  X2  >  0,  X3  <  0 

otherwise 


Then  F  is  not  injective  even  though  it  satisfies  all  the  other  restrictions 
mentioned. 
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